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ABSTRACT 

Studies of the formation of metal-free Population III stars usually focus primarily on 
the role played by H2 cooling, on account of its large chemical abundance relative to 
other possible molecular or ionic coolants. However, while H2 is generally the most im- 
portant coolant at low gas densities, it is not an effective coolant at high gas densities, 
owing to the low critical density at which it reaches local thermodynamic equilibrium 
(LTE) and to the large opacities that develop in its emission lines. It is therefore pos- 
sible that emission from other chemical species may play an important role in cooling 
high density primordial gas. 

A particularly interesting candidate is the H^" molecular ion. This ion has an LTE 
cooling rate that is roughly a billion times larger than that of H2, and unlike other 
primordial molecular ions such as H^ or HcH + , it is not easily removed from the gas 
by collisions with H or H2 . It is already known to be an important coolant in at least 
one astrophysical context - the upper atmospheres of gas giants - but its role in the 
cooling of primordial gas has received little previous study. 

In this paper, we investigate the potential importance of Hjj" cooling in primor- 
dial gas using a newly-developed H^ cooling function and the most detailed model of 
primordial chemistry published to date. We show that although is, in most cir- 
cumstances, the third most important coolant in dense primordial gas (after H2 and 
HD), it is nevertheless unimportant, as it contributes no more than a few percent of 
the total cooling. We also show that in gas irradiated by a sufficiently strong flux of 
cosmic rays or X-rays, H^" can become the dominant coolant in the gas, although the 
size of the flux required renders this scenario unlikely to occur. 



1 INTRODUCTION 



Over the past decade, we have made substantial progress in understanding how the very first stars in the universe formed. We 
know that in cosmological models based on cold dark matter (CDM), the first stars will form in small protogalaxies, with total 
masses of the order of 10 5 - 10 6 M^, an d that by a redshift z ~ 30 we expect to find at least one such star-forming system 
per comovmg Mpc 3 dYoshida et alJl2003h . We also know that although molecular hydrogen formation is inefficient, owing to 
the absence of dust, it is nevertheless the most abundant molecule in primordial gas, and is the main source of cooling at low 
densities, for temperatures between ~ 200 K and 8000 K. 



Furthermore, s imple semi-analytical estimates (e.g.. iTegmark et al.lll997l ). later confirmed by detailed simulations (e.g., 
Yoshida et ahlboOor ). demonstrate that H2 provides enough cooling in these small protogalaxies to allow the gas to collapse 
under the influence of gravity on a timescale comparable to its grav itational free-fall timescale, thereby a llowing star formation 
to occur. High-resolution, adaptive mesh simulations performed by Abel. Bryan, fc Norman ( 2000l . I2OO2I ) have taught us much 
about the dynamics of the gas in these first protogalaxies. They consider gas cooled only by H2, and find that gravitational 
fragmentation of the collapsing gas is inefficient and that therefore only a single, massive star will f orm during the collapse . 
This will then suppress further star formation through a variety of feedbacks (see the recent review by ICiardi fc Ferrarall2005T ) . 
Other simulations, making simplifying assumptions such as the adoption of spherical symmetry, have allowed us to examine 
the importance of various aspects of the physics of the gas which are currently difficult to treat in the high-resolution adaptive 
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mesh simulations (for instance, the dev e lopment of large optic al depths in the rotational and vibrational lines of H2 at high 
densities; see e.g. Omukai fc Nishi 19981 ; Ripamonti et al. 20021 ) . 

Much of the theoretical uncertainty that remains concerns the behaviour of dense gas during the later stages of gravita- 
tional collapse, and during the period of accretion that follows the formation of the first protostar. There is much that is still 
unknown here - for instance, a complete understanding of the mechanism by which the collapsing gas loses much of its initial 
angular momentum still eludes us - but in this paper we intend to focus on one relatively simple aspect: the identification of 
the dominant coolant(s) in the dense gas. 

As previously noted, H2 cooling has long been known to dominate at low densities, and it is frequently assumed that it also 
dominates at high densities. However, it is not at all obvious that this is actually the case. Two factors dramatically reduce the 
effectiveness of H2 as a coolant in high density gas. The first is the fact that the excited ro tational and vibrational levels of H2 
have small radiative transition probabilities, and hence long radiative lifetimes (r > 10 s; IWolniewicz. Simbotin. fc Dalgarno 
1998). This means that collisional de-excitation of excited H2 becomes competitive with radiative de-excitation at fairly low 



gas densities (n ~ 10 4 cm 3 ), and so as the number density n increases, the cooling rate of H2 quickly reaches its local 
thermodynamic equilibrium (LTE) value, given by 

Ah, ,iv] 



i,j>i 



AjiEj 



(1) 



where Aji is the transition probability for a transition from j — ► i, Eji is the energy of this transition, rij is the number 
density of H2 molecules in level j, computed assuming LTE, and where we sum over all bound levels i and over all bound 
levels j with energies greater than i. In the LTE limit, the cooling rate per H2 molecule is independent of the gas density, 
and is largely determined by the size of the transition probabilities. Since these are small, the LTE cooling rate is also small. 
Consequently, at high gas densities, a molecular species whose excited states have much shorter radiative lifetimes than those 
of H2 will provide far more cooling per molecule than H2. 

The second factor reducing the effectiveness of H2 cooling at very large n is the fact that the gas eventually becomes 
optically thick in the cores of the main rovibrational lines of H2. The effects of this cannot currently be treated fully in three- 
dimensional simulations, due to the high computational cost of solving the resulting r adiative transfer problem, but it has been 



dimensional simulations, due to the nign computational cost 01 solving tnc resulting radiative transfer problem, but it nas bcci 
modeled accurately in simple spherically symmetric, one-dimensional simulatio ns dOmukai fc Nishil [l998l ; Ripamonti et al 



2002l ; lRipamonti fc Abell20o4 ) and in an approximate fashion in three dimensions (jYoshida et al.ll2006l ) . These studies confirm 



that at densities n > 10 10 cm -3 , optical depth effects significantly suppress H2 cooling. 

Together, these factors combine to render H2 a fairly ineffective coolant in high density gas, despite the fact that at 
n > 10 10 cm -3 , several three-body H2 formation reactions 

H + H + H -> H 2 + H, (2) 
H + H + He -> H 2 + He, (3) 
H + H + H 2 -> H2 + H2, (4) 

rapidly convert almost all of the hydrogen in the gas to H2 ( Palla. Salpeter. fc Stahler 1983). It is therefore reasonable to ask 
whether cooling from any of the other molecular species present in the gas will become competitive with H2 cooling at these 
densities. 

One obvious possibility is deuterated hydrogen, HD. Its excited rotational and vibrational levels have radiative lifetimes 
that are about a factor of 100 shorter than those of H2, and so the HD cooling rate does not reach its LTE limit unt il 
n ~ 10 cm . It is also a far more effective coolant than H2 at low temperatures (T < 200 K; see e.g., iFlower et allboool ). 
This is due primarily to the fact that radiative transitions can occur between rotational levels with odd and even values of 
J, allowing cooling to occur through the J — 1 — > transition. The corresponding odd <-> even transitions in the case of H2 
represent conversions from ortho-H2 to para-H2 or vice versa, and are highly forbidden. Furthermore, at low temperatures 
the ratio of HD t o H2 can be significantly enhanced with respect to the cosmological D:H ratio by chemical fractionation (see 
e.g.. lGloverll2008l ). 

The role of HD cooling in early protogalaxies has been investigated by a number of authors. In the case of the earliest 
generation of protogalaxies, which form from very cold neutral gas that is never heated to more than a few thousand K during 
the course of the galaxy formation process, HD cooling appears to be unim portant, as the collapsed gas doe s not become cold 
enough for sufficient fractionation to occur to make HD cooling dominant ( Bromm. Coppi. fc Larson 20021) . 

The situation is rather different, however, in primordial gas cooling from an initially ionized state. In that case more 
H2 is formed, allowing the gas to cool to a lower temperature, at which point fractionation becom e s effective and HD 
cooling rapidly becomes dominan t (see e.g., Nakamura fc Umemural 20021 : Nagakura fc Omukai 2005 : Johnson fc Bromm 
20061 : IShchekinov fc VasilievlEooel ). However, the initial ionization required is much larger than that expected to be present 



in the earliest protogalaxies. 

Another molecule to ha ve attracted considerable attention is lithium hydride, LiH. This molecule has a very large dipole 
moment, /i = 5.888 debyes (jZemke fc Stwallevlll98oh . and consequently its excited levels have very short radiative lifetimes. 
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Therefore, despite the very low lithium abundance in primordial gas (xx,i = 4.3 x 10 10 , by number: see I Cvburtlliool ). it was 



thought for a time that LiH would dominate the cooling at very high densities (see e.g ., Lcpp & Shul] 



quantal calculations of th e rate of formation of LiH by radiative association (Da lgarno et al. 
199d : iBennett et alJliooih 



199' 



6: 



1984). However, accurate 



Gianturco fc Gori Giorgi 



Li + H -> LiH + 7, (5) 
have shown that the rate is much smaller than was initially assumed, while recent work by Defazio et al. ( 20051 ) has shown 



that the reaction 



LiH + H -> Li + H 2 



(6) 



has no activation energy and so will be an efficient destruction mechanism for LiH for as long as some atomic hydrogen 
remains in the gas. Consequently, the amount of lithium hydride present in the gas is predicted to be very s mall, even at very 
high densities, and so LiH cooling is no longer believed to be important ( Mizusawa. Omukai. fc Nishi 20051 ). 

ha ve att r acted little 
dl977l.ll978j). and 



In contrast to HD or LiH, the various molecular ions present in the gas, such as Hj~, Hj~ or HeH + 



attention. Some early work on H^" cooling in ionized primordial gas can be found in lSuchkov fc Sh chckinov 



its possible importance in hot, highly ionized conditions has recently been re-emphasized bv lYoshida et al.l (|2007l ). but aside 
from this, there has been little exploration of the role that cooling from these species might play in the evolution of primordial 
gas, presumably because the abundances of these species are expected to be small. It is this absence that the current paper 
attempts to rectify. 

We present here the results of a set of simulations of the chemical and thermal evolution of gravitationally collapsing 
primordial gas. These simulations use a very simple one-zone dynamical model for the gas, but couple this with a detailed 
chemical network and a comprehensive model of the various heating and cooling processes at work. Besides the coolants 
considered above, we include the effects of cooling from (in no particular order) H^", HD + , D^, Hj", H2D + , HD^, Dg~, HeH + , 
HeD + , HeJ, LiH + , LiD + , LiD and LiH J. We focus in particular on the possible role of H^". This ion has a very large number 
of excited rotational and vibrational levels that are energetically accessible at the temperatures of interest in primordial gas, 
and its vibrational levels have much shorter radiative lifetimes than those of H2 or HD. In LTE, its cooling rate per molecule 
is rou ghly 10 9 times larger than that of H2. It is known to be an importa nt coolant in planetar y atmospheres ( Miller et al. 
2000) and may also be an effective coolant in high density primordial gas ( Glover fc Savin 2006). Unlike ions such as Hf or 



HeH + , it does not react with H2, and is not easily destroyed by collisions with H, and so its abundance in high density gas 
should be large compared to the other molecular ions included in our model. 

The layout of this paper is as follows. In §2 we outline the numerical method used to simulate the thermal and chemical 
evolution of the protostellar gas. The chemical reactions included in the model are discussed in §3, and the thermal processes, 
in particular Hj" cooling, are discussed in §4. We present the results of our simulations in §5 and close with a brief discussion 
in 86. 



2 NUMERICAL METHOD 

We treat the thermal and chemical evolution of the gas using a one-zone model, in which the density is assumed to evolve as 



where ts = ^3iv/32Gp is the free-fall timescale of the gas and where 77 is an adjustable constant. In most of our simulations, 
we set 77 = 1, corresponding to free-fall collapse, but in a few runs we examine the effect of slowing down the collapse by 
setting rj < 1 (see i)5.7l for details). 

To follow the temperature evolution, we solve the energy equation 

de p dp 
dt = ~r? dt 



A + r > (8) 



where e is the internal energy density, p is the thermal pressure, A is the total cooling rate (which includes contributions from 
both radiative and chemical cooling, as outlined in |4| and F is the total heating rate. Since the temperature evolution is 
strongly coupled to the chemical evolution, we solve E quation [8] simultaneously with the chemical rate equations using the 
DVODE implicit ordinary differential equation solver ( Brown. Byrne, fc Hindmarsh 1989). To model the chemistry 



we use 



an extensive chemical network consisting of 392 reactions among 30 atomic and molecular species. Table [T] lists all 30 species 
considered. Tables IA1HA141 in Appendix A list the reactions included in this network, broken down by the type of process 
involved. These tables also give details of the rate coefficient or rate adopted for each reaction and the source of the data. In 
these tables and elsewhere in the paper, T is the gas temperature in K, T3 = T/1000 K, and T c is the temperature in units 
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Table 1. Initial fractional abundances in our reference calculation 



Species 


Initial abundance 


Species 


Initial abundance 


Species 


Initial abundance 


c~ 


2.2 X icr 4 


HD+ 




HeD+ 




H+ 


2.2 X 10" 4 


HD 


6.2 X lO" 11 


He+ 




H 


0.99978 


»t 




Li+ 


2.2 X lO" 10 


H" 




D 2 


1.6 x 10~ 15 


Li 


2.1 x lO" 10 






H 2 D+ 


0.0 


Li" 




H 2 


2.4 x 10~ 6 


HD+ 


0.0 


LiH+ 


0.0 


H 3 


0.0 


D 3 + 


0.0 


LiH 


0.0 


D+ 


5.7 x 10" 9 


He+ 


2.8 x 10~ 26 


LiD+ 




D 


2.6 x 10~ 5 


He 


8.3 x 10" 2 


LiD 


0.0 


D~ 




HeH+ 




LiH+ 





Notes: The quoted values are fractional abundances relative to the number density of hydrogen nuclei. 

Chemical species listed without initial abundances are assumed to be in chemical equilibrium, as described in the text. 



of eV. (Note that we assume, both here and throughout, that the kinetic temperature of the electrons is the same as that of 
the atoms and molecules). 

For many of these species, we followed the full time-dependent, non-equilibrium chemistry, but in some cases - the ions 
H - , D~, Li~, Hj, HD + , D^, HeH + , HeD + , HeJ, LiD + and LilTj" - chemical equilibrium is reached very rapidly, on a 
timescale of the order of t cq = 1/ (kdestn) ~ 10 9 nT 1 s, where fcdest ~ 10 -9 cm 3 s _1 is a characteristic destruction rate coefficient 
and n is the number density of hydrogen nuclei. At the gas densities considered in this work (1 < n < 3 x 10 13 cm -3 in 
the majority of our simulations), this chemical equilibrium timescale for these rapidly reacting ions is orders of magnitude 
smaller than the free-fall timescale of the gas, which is approximately 1.5 x lO^n -1 / 2 s _1 . It is therefore sufficient to use 
the equilibrium abundances for these ions. Some of the species remaining in our non-equilibrium model may also be close to 
chemical equilibrium during a large portion of the collapse, but are included in the non-equilibrium treatment because we 
cannot be sure that they always remain in equilibrium. Further details of our chemical network are given in Sj3] 

We assume elem ental abundan ces relative to hydrogen of 0.083, 2.6 x 10 -5 and 4.3 x 10 -10 for helium, deuterium and 
lithium, respectively I Cvburtll2004h . In £ ]5.5l we explore the effects of reducing the deuterium and/or the lithium abundance to 
zero, in order to asses the impact of the deuterium and lithium chemistry on the evolution of the gas. The initial abundances 
of the various molecular and ionic species in our standard model are summarized in Table [T] The values used for the initial 



H" 1 ^ H e + and H2 abundances, and the ratio of ionized to neutral lithium are taken from the calculations of IStancil et al 



(199. 



3) - specifically their model V. Deuterated species are assumed to have abundances that are a factor (2.6 x 10 5 ) JVd 



smaller than the abundances of the undeuterated equivalents, where No is the number of deuterium nuclei in the species in 
question. The electron abundance is computed assuming charge conservation. To assess our sensitivity to these initial values, 
we have also run several simulations with different initial abundances; the results of these runs are discussed in 3531 



For most simulations, we adopt an initial density rii = 1 cm -3 and an initial temperature T\ — 1000 K. The effects of 
altering n\ and Ti are examined in 35.51 All of the simulations are run until the density exceeds nt = 3 x 10 13 cm -3 . At higher 
densities, collision-induced emission from H2 quickly comes to dominate the cooling, and the minor species considered here 
are unlikely to be important coolants in this very high density regime. 

Finally, we ran all of our simulations starting at a redshift z = 20. However, the main influence of the redshift is to set a 
minimum temperature for the gas (since the gas cannot cool radiatively to below the CMB temperature, Tomb). As the gas 
temperature T 3> Tcmb at the densities of interest in this paper, we do not anticipate that changing 2 by a moderate amount 
will significantly affect our results. 



3 CHEMISTRY 



The chemical evolution of the gas is modeled using a chemical network consisting of 392 reactions amongst 30 neutral and 
ionized species. A list of all of the reactions included can be found in Tables IAT1IA14I This network, which to the best of our 
knowledge i s the l argest used to date f or the study of primordial gas chemistry, is based in part on previous compilat i ons b y 
Abel et all rtl997h.ISta.ncil et al.1 rtl998Tl.lGalli fc Pallal (|l99g| . l2002h . lLepp. Standi, fc Dalgarnol |2002h , IWang fc Stancill rt2002t ). 
and Walmslev. Flower, fc Pineau des Forets ( 20041 '). supplemented with additional reactions drawn directly from the chemical 
literature, as well as some whose rates have not (to our knowledge) been previously discussed. These latter are generally 
rates involving one or more deuterium nucl ei in place of hydrog en nuclei; and in estimating rates for these reactions, we have 
generally followed the same procedure as in lStancil et alj (1998): for a non-deuterated reaction with a reaction rate coefficient 
that has a power-law temperature dependence k oc T m , we have estimated a rate coefficient for the deuterated reaction by 
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multiplying this rate coefficient by a scaling factor (/xh/md)™, where /ih and /j,o are the reduced masses of the reactants in the 
non-deuterated and deuterated reactions, respectively. Some notable exceptions to this strategy are discussed in more detail 
in section i]3.1l below. 

For reactions where the presence of a deuteron increases the number of distinguishable outcomes and where no good 
information exists on the branching ratio of the reaction, we have assumed that each outcome is equally likely. An example of 
this is the dissociative attachment of HD with e~ (reactions AD6-AD7), which can produce either H and D~ or H" and D, in 
contrast to the dissociative attachment of H2 with e~ (reaction AD5) which can only produce H~ and H. For this particular 
example, this assumption gives branching ratios of 50% for reactions AD6 and AD7, respectively. 

In spite of the size of our chemical network, there remain a number of chemical processes that are not included. These 
are discussed in A3. 41 along with our justifications for omitting them. 



3.1 Discussion of selected reactions 

3.1.1 Photorecombination of H and He (reactions PR1 & PR3, Table \A2]) 

We assume that the ionizing photons produced by the recombination of H + to ground state H are immediately consumed by 
the ionization of atomic hydrogen (the on-the-spot approximation), and so we use the case B rate coefficient for hydrogen 



recombination (IFerland et al.lll992f ). We note that although the fractional abundance of atomic hydrogen becomes small at 
densities greater than 10 10 cm -3 , the number density of atomic hydrogen remains considerable, and so case B remains a good 
approximation . 

We also use the on-the-spot approximation to treat He + recombination, but in this case the net recombination rate is 
larger than the case B rate, as some of the photons produced by recombination directly into the n = 1 ground state are lost 
through photoionization of H rather than He. For He + recombinations directly into the ground state, occurring in primordial 
gas with a low fractional ionization and low molecular a bundance, approx imately 68% of the resulting photons are absorbed 
by H, with the remaining 32% being absorbed by He l|Osterbrocklll989h . Therefore, the effective He + recombination rate 
coefficient in these conditions is given by 

fcpR3 = 0.68fcpR3,rr,A + 0.32fcpR3,rr,B + fcpR3,di Cm 3 S~\ (9) 

where fcp R 3, r r,A an d fcpR3, rr ,B are the case A and case B rate coefficients, and fcpR3,di is the dielectronic recombination rate 
coefficient. This formula becomes incorrect once the H2 fraction of the gas becomes large, but as this occurs only at densities 
n > 10 10 cm -3 at which the He + fractional abundance is negligible, the error that is introduced by using this prescription for 
fcpR3 throughout the simulation is unimportant. 

It is also necessary to take account of the photoionization of H caused by the He + recombination emission. In addition 
to the contribution coming from He + recombination into the n = 1 ground state, there is an additional contribution made 
by photons produced during transitions from n = 2 to n = 1 in He; in other words, even pure case B recombination of 
He + produces H-ionizing photons. The proportion of case B recombinations that yield photons capable of ionizing hydrogen 
depends upon the relative populations of the n = 2 singlet and triplet states, and hence upon the electron density, but in th e 
low density limit, 96% of all recombinations to excited states produce photons that will ionize hydrogen 1 Qsterbrock 19891 ) . 



To model the effects of these photons, along with those produced by recombination direct to the n = 1 ground-state and by 
dielectronic recombination, we include in our chemical network a local H ionization rate per unit volume R p i, with a value 

Rpi = [0.68(fcp R 3,rr,A — fcpR3,rr,B) + 0.96fcpR3,rr,B + 2fc PR 3 idi ] n c n Ho + Cm" 3 S~ , 

= [0.68fcp R3 ,rr,A + 0.28fc P R3,rr,B + 2fc P R 3 ,di] «- c n Hc + Cm" 3 S~\ (10) 

where the three terms in brackets on the first line correspond to the contributions from recombination direct to the n = 1 
ground state, pure case B recombination (i.e., recombination to all states n ^ 2), and dielectronic recombination, respectively. 
(Note that every dielectric recombination produces two photons capable of ionizing hydrogen: one due to the radiative 
stabilization of the process , and one as the captured electron cascades to the Is level). If the electron density is large 



(n e > 10 3 cm~ 3 ; lciegjll987h . then more helium recombinations will result in two-photon tra nsitions from 2 S-l S, reducing 



the number of photons produced that are capable of ionizing hydrogen l|Qsterbrocklll989l ). However, the effect on R pl is 
relatively small, and in any case, we do not expect to encounter large abundances of He + in dense gas in the particular 
scenario that we are investigating. Therefore, adopting this simplified treatment at all n should be sufficient for our purposes. 



3.1.2 Dissociative recombination of (reactions DR4 & DR5, Table [Al 

For a long time, considerable disagreement has existed on the sub ject o f the H^" dissociative recombination rate. Measurements 
of the rate in merged beam experiments (e.g., ISundstrom et al.lll994T ) typically give values of the order of 10 -7 cm 3 s _1 for 
the rate coefficient at temperatures near room temperature, while measurements made in flowing afterglow experiments (e.g., 
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1 , an order of magnitud e smaller. At the same time, most 
theor etical calculations have indicated a smaller rate still, of the order of 10 -11 cm 3 s _1 (see lOrel. Schneider, fe Suzor-Weiner 



Smith fc Spanell 1993a, b) often give values of the order of 10 



2000, and references therein), which is in complete disagreement with the experimental measurements. 



However, it has recently become clear that three-body e ffects play a highly important role in the recombination of Hg 
in flowing afterglow experiments (see e.g., Glosfk et al. 20051 ). When proper allowance is made for these effects, the inferred 
tw o-body recombina tion rate is in good agreement with the results of the merged bea m experiments (see th e discu ssion 
in Glosfk et al. 2007 ). Moreover, recent theoretical calculations of the rate coefficient by Kokoouline fc Greene! ( 20031 ) that 



account for the effects of Jahn- Teller coupling between the electronic and vibrational degrees of freedom produce a result that 
is in good agreement with the experimental measurements, although disagreement at the level of a factor of two or so is still 
present at some energies. 

In our calculations, we theref ore take our value fo r the total Hjj" dissociative recombination rate coefficient from the recent 
ion storage ring measurements of lMcCall et al.l (|2004r i. To convert this total rate coefficient - the sum of the rate coefficients 
for reactions DR4 and DR5 - into a rate coefficient for each indivi dual reaction, we a dopt a branching ratio of 0.25 for reaction 

~ 1 19951 ). 



DR4 and 0.75 for reaction DR5, based on the measurements of Datz et al 



Strictly speaking, these values are only 
appropriate for temperatures T < 3000 K, but in practice the behaviour of the gas is not particularly sensitive to the values 
chosen. 



3.1.3 H3 formation by radiative association (reaction RA18, Table YASv 

The rate coefficient we quote for reaction RA18, the radiative association of H2 and H + to fo rm H^j~, was taken from the 
tudy of Gerlich fc Horninel 1 1992^ . and was the rate coefficient quoted by Galli fc Palla (1998) for this reaction. However, 

n 3 s _1 for this reaction. In 3531 we examine the effect of 



Stancil et al 



1 1998) prefer a much smaller rate coefficient of 10 



adopting this smaller rate coefficient. 



3.1.4 H2 formation by associative detachment of H (reaction AD1, Table 



The rate of this reaction is quite uncertain, and we have shown in prev ious work that this uncertainty can lead in some cases 
to a substantial uncertainty in the amount of H2 formed in the gas ( Glover. Savin, fc Jappsen 20061 ). However, we do not 
expect this uncertainty to significantly affect the results in this paper. At the high densities at which H3" cooling is potentially 
important, the dominant H2 formation pathway is three-body formation. This can produce a much larger molecular fraction 
than is possible via two-body reactions, and so uncertainty in the amount of H2 produced via the H - ion has a negligible 
impact on the evolution of the gas at high densities. In our simulations, we adopt a def ault value for t he ra te coefficient 



for reaction AD1 of fc AD i = 1.5 x 10" 9 T 3 " 01 cm 3 s~\ where T 3 = T/300 K, taken from lLaunav et all (|l99ll ). In 351)1 we 



demonstrate that our results are insensitive to this choice. 



We note also that the other reaction discussed in iGlover. Savin, fc Jappsen! (|2006l ). the mutual neutralization of H~ 
by H + (reaction MN1, Table IA8|) . is no longer a source of significant uncertainty in chemical models of primordial gas. 
Recent measurements of the cross-section for this reaction at astrophysic ally-relevant energ i es ma de by X. Urbain (private 
communication) yiel d values in good a greement with those obtained by Fussen fc Kubach ( 1986T ) and used as a basis for 
the rate coefficient of ICroft et al.l ( 199E). These measurem ents strongly suggest that the earlier measurements of the mutual 
neutrali zation cross-section made by Mpseleyet_aL 1 1970[ ) were somehow in error, and that rate coefficients based on them 



(see e.g. iDalgarno fc Lepplll987l : [Galli fc Pallal l 19981 ) ar e incorrect. The error in the rate coefficient for this reaction has thus 



been reduced from the order of magnitude discussed in Glover. Savin, fc Jappsen ( 20061 ) to an uncertainty of about 50% (X. 
Urbain, private communication). 



3.1.5 Collisional dissociation of H2 (reactions CD9-CD12, Table p47| ) 

The rate coefficients for the collisional dissociation of H2 by H (CD9), H2 (CD10), He (CD11) and e~ (CD12) are represented 
by functions of the form 

log h = ( n J n " ) log fci,LTE + ( : 1 —, J log fc, v =o, (11) 

\1+ n/n CI J \ 1 + n/n cl - J 

where n is the number density of hydrogen nuclei, ki is the collisional dissociation rate for collisions with species i, and fc v= o,i 
and fci/rE,i are the rate coefficients for this reaction in the limits in which all of the H2 molecules are in the vibrational 
ground-state (appropriate in low-density gas), or have their LTE level populations (appropriate for high density gas). The 
critical density, n cr , for H2 collisional dissociation in a gas containing a mix of H, H2, He and electrons is not well determined. 
For simplicity, we therefore assume that it is given by a weighted harmonic mean of the (better known) critical densities 
corresponding to reactions CD9, CD10, and CD11 considered individually, i.e., 
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l 

1 + XHc 



■I'll 



+ 



2.x 



II 2 



n C r,H «cr,H 2 ^cr,Hc 



(12) 



where ih, xh 2 , and XHe are the fractional abundances of H, H2 and He relative to the total number of hydrogen nuclei, we 
use the approximation that xn + 2xn 2 — 1> and where 

^ ^ _ To n n a 1 a l, — , n 007 f\~~ T<. \21 „ m — 3 



Wcr,H 2 



n cr ,Ho = 



dex [3.0 -0.416 log T 4 - 0.327 (log T 4 ) 2 ] cm" 3 , (13) 
dex [4.845 - 1.3 log T 4 + 1.62 (log T 4 ) 2 ] cm" 3 , (14) 
dex [5.0792 {l.O - 1.23 x 10" 5 (T - 2000)}] cm" 3 , (15) 

with T 4 = r/10000 K. The expression for n CT .H is f rom | Lepp fc Shulll 1 19831 ) . but has been decrea sed by an order of magn itude, 
as recommended by Martin. Schwar z i- ^JMandy ( 19961 ). The expression for n cr ,H 2 comes from Shapiro fc Kane (1987), and 
the expression for n cr ,Hc comes from lDove et al.l (|l987h . Note that this expression for the critical density assumes that in high 
density gas, n c < hh, so that electron excitation of H2 does not significantly affect the value of n cr . Other forms of averaging 
to obtain n cr are possible, of course, but we would not expect our results to be sensitive to our particular choice here, as any 
differences will only be seen for densities n ~ n CI , and in our simulations, gas at these densities is always far too cold for 
collisional dissociation of H2 to be important. 

To ensure that our adopted collisional dissociation rate coefficients and three-body H2 formation rate coefficients are 
consistent, we used the fact that in LTE, the equilibrium constant K obeys 

K = fc T Bl ACD9 = fcTB2/fccD10 = fcrB3/fccDll (16) 



and varies with temperature as IjFlower fc Harris! 120071 ) 



K = 1.05 x 10" 22 T" l515 exp 



^ 52000 ^ 



(17) 



to determine values for the rate coefficients of reacti ons CD9, CD10, and CP U in the LTE lim i t. How ever, we also ran some 
test simulations where we used rate coefficients from Lepp fc Shull ( 19831 ) and Shapiro fc Kane ( 1987 ) for reactions CD9 and 
CD10, regardless of the value of fcxBi or A;tb2. These simulations produced almost identical results, demonstrating that our 
results here are insensitive to our treatment of H2 collisional dissociation. 



.1.6 Collisional dissociation of HD and D2 (reactions CD13-CD20, Table 



For col lisions with electrons, accurate rate coefficients are available from lTrevisan fc Tennyson! (|2002ar ) and lTrevisan fc Tennyson 
I 2002b|). For collisions with H, H2, or He, however, we are unaware of a treatment in the literature. We have therefore assumed 
that the rate coefficients of these reactions in the v — and LTE limits are the same as for the corresponding H reactions 
(nos. CD9-CD11). For D2, we have also adopted the same value for the critical density, while for HD, we have increased it 
by a factor of 100 to account for the larger radiative transition probabilities. Note that although these rate coefficients are 
highly approximate, this probably does not introduce much uncertainty into the chemical model, as reactions 1X18 and 1X20 
(Table lAlO[l become effective at much lower temperatures and therefore dominate the destruction of HD and D2 in warm gas. 



3.1.7 Three-body H2 formation (reactions TB1-TB3, Table Wty 



Although unimportant at low densities, three-body reactions are the dominant source of H2 in high density primordial gas, 
and so these reactions represent an important part of our chemical network. Unfortunately, rate coefficients for these reactions 
are, in general, not known to a high degree of accuracy. For three-body collisions in which atomic hydrogen is the third body 
(react ion TBI), the situation is particularly bad. One commonly adopted rate coefficient is that of lPalla. Salpeter. fc Stahler 
(1983), who quote a rate coefficient 



fc T Bi,pss = 5.5 x 10" 2H T 



-1 6 - 

cm s 



(18) 



for this rea ction, based on experimental wo rk bv lJacobs. Giedt. fc Cohen! (| 1967T ) ■ Also in common usage is the rate coefficient 
adopted by I Abel. Bryan, fc Norman (2002), which is 



A;tbi,abn = 1.14 x 10" 31 T" - 38 cm 6 s" 1 T< 300K 

= 3.90 x 10" 30 T" 100 cm 6 s" 1 T > 300 K. (19) 

The low temper ature portion of this rate coef ficient is based on lOrell |l987h . w hile the high tempe r ature portion is an 
extrapolation by Abel. Bryan, fc Norm an (2002). This reaction is also discussed by Cohen fc Westberj 1 19831 ) in their large 
compilation and review of chemical kinetic data. They summarize a large number of different experimental measurements and 
argue that the precision of the data does not justify anything more elaborate than a constant rate coefficient 
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fc T Bi,cw = 8.8 x lQ~ Ai cm 6 s . 

Another possibility is found in ISchwenkel(|l990h , who gives calculated values at T = 3000K and T = 5000K of 1.4 x 10 



(20) 



and 8.2 x 10 33 cm 6 s 1 , resp ectively. These values are roughly an order of magnitude larger than those given by the 



Abel . Bryan, fc Norman! (120021 ) rate coefficient, and about 30% lower than the values given by the lPalla. Salpeter. fc Stahler 
( 1983. ) rate coefficient. 

More recently, Flower fc Harris! 1 2007 ) have argued in favor of a rate coefficient 



fcrm.FH = 1-44 x 10~^ D T 



-1.54 6 - 

cm s 



(21) 



which they derived from the rate coefficient of the inverse process (H2 collisional dissociation by atomic hydrogen, re- 

action CD9) by using th e prin ciple of detailed balance. This rate coefficient is approximately six time s larger than the 

Palla. Salpeter. fc Stahler, i| 19831 ) rate coefficient at T — 1000K, or approximately ninety times larger than the lAbel. Bryan, fc Norman! 
1 20021 ) rate coefficient. Unfortunately, the accuracy of a rate coefficient derived using detailed balance depends upon the ac- 



curacy with which the rate coefficient of the inverse process is known. In this case, that accuracy is poor, as the H2 collisional 
dissociation rate coefficient is not well constrained by exper iment at low temperatur e s (i.e. , T < 2000 K) owing to its small 
size at these temperatures. Flower fc Harris! 1 2007 ) used the Jacobs. Giedt. fc Cohen d 19671) fit to the collision al dissociation 



Martin. Schwarz. fc Mandvl dl996.). then a much 



rate coefficient, but if we instead use the calculated rate coefficient from 
smaller three-body H2 formation rate coefficient is obtained, which can be fit to within ~ 20% by (|Gloverll2008l) 



, _ „ .. r ,-31 T ,-0.464 6 - 

ktbi.gl = 7.7 x 10 1 cm s 



(22) 



There is thus an uncertainty of almost two o rders of magnitude in the rate coefficient for reaction TBI. In our simulations, 
we adopt the Abel. Bryan, fc Norman 1 20021 ) rate coefficient as our default value, but in EI5.6I we examine the effect of using 
a different rate coefficient. 

The rate coefficient for three-body H2 formation in collisions where H2 is the third body (react ion TB2) is known with 
greater precision, but nevertheless substantial uncertainty remains. iPalla. Salpeter. fc Stahler |l983l ) quote a rate coefficient 



A:tb2,pss = 6.9 x 10 30 T 10 cm 6 s 1 



(23) 



for this reaction, again taken from Jacobs. Giedt. fc Cohen (1967), while Cohen fc Westberg ( 19831 ) recommend instead 



fcTB2,cw = 2.8 x i(r 31 T~ ' 6 cm 6 s~\ 



(24) 



Flower fc Harris! 1 2007 ) assume that the ratio of fcTB2 to &tbi is the same as that measured by Jacobs. Giedt. fc Cohen ( 1967 ) 
i.e., one-eighth. Therefore, their rate coefficient for reaction TB2 is: 



7 -I o , n — 27rrn— 1.54 6 - 

ktb2,fh = 1.8 x 10 T cm S 



(25) 



The same assumption applied to the rate coefficient from Glover (2008) gives 



7 n enc w i a-32t-,-0.464 6 -1 

«tb2,gl = 9.625 x 10 T cm s . 



(26) 

Finally, calcul ations by Schwenke! 1 19881) using orbit al resonance theo ry find a rate coefficient that is about a factor of two 
lower than the Cohen fc Westberg! 1 19831 ) values, but Schwenke (1990) shows that one of the assumptions underlying his own 
orbital reson ance calculations is i nvalid , and provides revised values, obtained from a master equation approach, that agree 
well with the lCohen fc Westberg, (|l983h recommendation. 

These rate coefficients agree to within a factor of a few at T = 3000 K, co nsistent with the scatter in experimental 
determinations of the rate coefficient at this temperature ( Cohen fc Westberg 19831). but differ m ore by more than an order of 
magnitude at low temperatures. In our simulations, we use the lPalla. Salpeter. fc Stahlerl l|l983l ) rate coefficient as our default 
value, but we examine in i|5.6l the effect of altering /ctb2- 

We also included three-b ody formation of H2 via collisions with He (reaction TB3), using a rate coefficient from 
Walkauskas fc Kaufman! (119751 ). This reaction has not been included in previous treatments of the evolution of dense pri- 
mordial gas and so we wished to assess its effects. We found that reaction TB3 could be responsible for anywhere between 
0.1% and 10% of the total three-body H2 formation rate, depending on the temperature, the H2 abundance and the choice of 
rate coefficients for reactions TBI and TB2. Moreover, this estimate does not take into account the uncertainty in the rate 
of reaction TB3, which we have been unable to quantify, but which should probably be assumed to be comparable to the 
uncertainty in the other three-body rates. Thus, although it probably never dominates, reaction TB3 should be included if 
accurate modeling of H2 formation in dense gas is desired. 



3.1.8 Deuterated three-body reactions (TB4-TB9, TB11-TB13, TB17-TB31, & TB34-TB35, Table \Mj) 

In view of the large uncertainties present in the rate coefficients of many of the three-body reaction rates (particularly for 
reactions TBI and TB2, as discussed above), we consider that the most prudent course of action when estimating rates for the 
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deuterated forms of these reactions is simply to adopt the same values as for the non-deuterated reactions. Any uncertainty 
introduced by this assumpti on is likely dwarfed by the uncertainties arising from our poor knowledge of the non-deuterated 
reaction rates. We note that Flower fc Harris! 1 2007 ) follow a similar course of action in their study of three-body H2 and HD 
formation in primordial gas. 



3.1.9 Destruction of D2 by collision with H (reaction 1X20, Table BTfll) 

Our fit to the data collated by Mielke et al. ( 20031 ) for this reaction is accurate to within a few percent over the temperature 
range of the tabulated data, 200 ^ T ^ 2200 K. Outside of this range, our fit may be significantly inaccurate (although at 
low temperatures, the reaction rate is small enough that any inaccuracy is unlikely to be important). 



3.1.10 Photodissociation of H 2 and HD (reactions BP7 & BP 8, Table VKWl 

Table IAT21 lists the rates of these reactions in optically thin gas, given our assumed incident UV spectrum (see i]3.2[) . In 
optically thick gas, however, self-shielding of H2 by H2 and HD by HD can significantly reduce both of these rates, by factors 
/sh,H 2 an d fah. HD respectively. In s t atic, isothermal gas, these self-shielding factors can be calculated approximately using the 
prescription of iDraine fc Bertoldi (199a) together with an appropriate set of molecular data, provided that one knows the 
H2 and HD column de nsities. In gas which is in motion, with internal velocities comparable to or larger than the thermal 
velocity of the gas, the IDraine fc Bertoldi jl99rj ) treatment brea ks down, and one m ust use approaches that are either less 
accurate or more computationally expensive, as discussed in lGlover fc Jappsenl (|2007l ). However, in the one-zone calculations 
presented here, we know neither the H2 and HD column densities, nor the velocity structure of the gas, and so including even 
a highly approximate treatment of the effects of self-shielding is problematic. In our runs with a non-zero UV background, 
we therefore consider two limiting cases: one in which self-shielding is highly efficient and / s h.H 2 — /sh,HD = throughout the 
run, and one in which it is ineffective, and we remain in the optically thin limit throughout (i.e., / s h,H 2 — /sh.HD = 1). These 
two limiting cases bracket the true behaviour of the gas. 



3.1.11 Formation and destruction o/LiHj (reactions RA20, DR19, DR20, DR21 & CD26) 

To date, the LiHJ" ion has attracted little attention in the astrochemical literature. iKirbv fc Dalgarnol (|l978h considered the 
reaction chain 

Li++H 2 -> LiH++ 7 , (27) 
hmt + e~ -> LiH + H (28) 

(reactions RA20 and DR20 in our chemical model) as a possible source of LiH in the interstellar medium, but showed that 
even if the rate coefficient for the radiative association were assumed t o be very large (fcRA 20 ~ 10~~ 16 cm 3 s" 1 ), the resulting 
LiH abundance would be far too small to be observable. More recently, Standi et al. (1996) considered the LiH J ion in their 
comprehensive study of the lithium chemistry of the primordial intergalactic medium, but again reached the conclusion that 
its abundance would be very small, and so chose not to include it in their chemical model. However, to the best of our 
knowledge, there has been no previous investigation of the role that this ion may play in regulating the fractional ionization 
of very dense primordial gas. 



Previous work focussing on modelling the fractional ionization at high densities ijMaki fc Su sa 2004, 2007]) , in the context 
of the study of ambipolar diffusion in dense population III prestellar cores, has shown that once the free electron fraction 
falls below x ~ 10~ 10 , ionized lithium takes over from ionized hydrogen as the primary positive ion in the gas. It is there- 
fore important to ensure that all of the major loss routes for Li + are represented in the chemical model. In addition to 
photorecombination (reaction PR4), Li + can be removed from the gas by a number of reactions with atomic or molecular 

(29) 
(30) 
(31) 
(32) 
(33) 
(34) 

Most of these processes are highly endothermic, and so are not competitive with photorecombination even when x is small. How- 
ever, the two radiative association reactions are exothermic and deserve closer scrutiny. Radiative associati on with atomic hy- 
drogen (reaction RA10) has been included in a number of previous models of primordial gas chemistry (e.g. Stancil et al.lll99a ; 



hyd 


rogen: 




Li" 


^+H - 


Li + H+, 


Li" 


^+H - 


LiH++ 7 , 


Li+ 


+ H2 -> 


Li + H+, 


Li+ 


+ H 2 -* 


LiH + + H, 
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Galli fc Palla 1998h, and accurate q uantal calculations of the rate coefficient for this reaction are available ( Dalgarno et al,. 
199rj ; iGianturco fc Gori Gio rgi 1996 . However, at the densities of interest in the present case, efficient three-body formation 
of H2 ensures that the hydrogen is primarily molecular rather than atomic, and so renders radiative association with H2 (reac- 
tion RA20) the more important reaction. Unfortunately, we have been unable to locate any calcul ation of the rate coe fficient 



10" 



while Standi et al 



! 1996h 



quote 



of this reaction. Kirbv fc Dalearnol 1 1978h quote an upper limit of A:ra20 
an upper limit of &ra20 = 10~ 17 cm 3 s" 1 , but the true rate coefficient could be orders of magnitude smaller. In our reference 
model, we make the conservative assumption that the rate coefficient is of the same order of magnitude as that for reaction 
RA10, and hence adopt a value &ra20 = 10~ 22 cm 3 s" 1 . We investigate the effects of adopting a larger value in section i ]5,6,4 1 
The inclusion of reaction RA20 in our chemical network necessitates the inclusion of additional reactions: the dominant 
destruction mechanisms for LiHj. Unfortunately, there has been very little theoretical or experimental study of any reactions 
involving LiHj and so it is not even clear which processes dominate. The best studied destruction process is dissociative 
recombination (reactions DR19, DR20 and DR21): 

LiH^+e" -> Li + H 2 (35) 
-> LiH + H (36) 
-> Li + H + H. (37) 

Thomas et"ZI l|2006l ) have studied th s process experimentally using the CRYRING heavy ion storage ring, and have reported 
preliminary results regarding the branching ratio of the reaction, but have not yet reported any value for the total rate. C. 
Greene and collaborators are currently involved in a theoretical calculation of the total rate coefficient, but again have yet to 
publish any results. However, their preliminary findings suggest a total rate coefficient that is about 2.5-3 times larger than 
that for the dissociative recombination of H^j". (C. Greene, private communication). We therefore adopt a total rate coefficient 
2 x 10 -7 (3I0) ^ cm 3 s _1 for the dissociative recombination of LiHj , and use the values quoted by Thomas et al. (2006) 
for the branching ratios. 



Other reactions that could be important destruction mechanisms for LiHj include 



LiHJ + H 
LiH^ + H 
LiH J + H 2 
LiHj + He 



LiH + + H 2 
Li+ + H 2 + H 
Li+ + H 2 + H 2 
LiHe + +H 2 . 



(38) 
(39) 
(40) 
(41) 



None of these reactions appears to have previously been studied in the astrochemical literature, and so their rates are unknown. 
For simplicity, therefore, we include only a single representative example from this set of reactions, namely the collisional 
dissociation of LiH J by H2 (reaction CD26). As we are primarily interested in the role of LiH J within the highly molecular 
dense core, it is likely that this reaction will dominate, unless its rate coefficient is unusually small. In our reference model, 
we adopt a rate coefficient of fccD26 = 1.0 x 10~ 9 exp (— 22|^) cm 3 s _1 for this reaction; however, in £15.6.41 we examine the 
effect of adopting a smaller value. 



3.2 Photochemistry 

To compute rates for the photochemical reactions listed in Table IAT2I we assume that the gas is illuminated by an external 
background radiation field with the spectral shape of a 10 5 K black-body at energies hu < 13.6 eV, and which is zero at higher 
energies. This choice of spectrum is motivated by the fa ct that the brightest population III stars are expected to have high 
effective temperatures, T e s ~ 10 5 K (jCoiazzi et aljbooch . while the cutoff at the Lyman limit is intended to account for the 



effects of absorption by neutral hydrogen in the intergalactic medium. We quantify the strength of this background radiation 
field in terms of the flux at the Lyman limit, J(u a ) = 10~ 21 J21 erg s _1 cm -2 Hz -1 sr _1 . The rates listed in Table lA12l are 
computed assuming that J21 = 1.0, but scale linearly with J21 and so can easily be rescaled for other values of the background 
radiation field strength. 



3.3 Cosmic rays 

If cosmic rays are present, then they will directly ionize some species and indirectly photoionize and photodissociate others. 
Direct ionization is simple to treat, and the appropriate rates are listed in Table [A13l normalized by the cosmic ray ionization 
rate of atomic hydrogen, £h, which we treat as a free parameter. However, the indirect effects of the cosmic rays are harder 

to model accurately. 

The basic physics is straightforward, and was first discussed bv lPrasad fc Tarafdar (1983). They noted that the secondary 
electrons produced by cosmic ray ionizations are energetic enough to excite the electronic states of H2, and that the subsequent 
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radiative decay of these excited states would produce ultraviolet photons. In the Galactic context, the mean free path of these 



photo ns is small, and so the cosmic-ray induced photochemistry can be modelled as a purely local process (see e.g-. lGredel et al 



1989). 



In the population III star formation context in which we are interested, however, there are two main factors that complicate 
matters. First, the H2 fraction in the gas is small at densities n <C 10 10 cm -3 , and so most of the secondary electrons produced 
by the cosmic rays lose energy by exciting and ionizing atomic hydrogen, rather than molecular hydrogen. Second, an accurate 
treatment of the propagation of the photons produced by excited H and H2 is far more involved than in the Galactic ca se. The 



continuum opacity of metal- free gas is very small at most densities of interest ([Lenzuni. Chernoff. fc Salpeterlll99lr ). owing 
to the absence of dust absorption, and so the majority of the photons produced by cosmic-ray induced excitation of H2 have 
large mean free paths. On the other hand, Lyman-Q photons produced by the excitation of atomic hydrogen have small mean 
free paths, but scatter many, many times before escaping the gas (see e.g., Diikstra. Haiman. fc Spaansl 20061 ) . In neither case 



is it a good approximation to assume that all of the photons are absorbed locally in the gas, and so the simple local treatment 
developed for Galactic dark clouds no longer applies. 

If we consider only the effects of emission from H2, then to compute _Rx, the photoionization (or photodissociation) rate 
per unit volume of species X, one must use an equation of the form 

1 f°° f e(v x 'V~ T( "' x '' x) 

*x(x) = -/ o *,(v)j v (' |x ;_ x|2 dx'd, (42) 

where the volume V that we integrate over corresponds to the entirety of the protogalactic core, and where t(v, x',x) is the 
optical depth of the gas between point x and point x' at a frequency v. The photon emissivity e(v, x') is given in this case by 

e(u, x') = Ph 2 (^)Ch 2 wh 2 , (43) 

where Ph 2 i v ) dv is the probability that the cosmic ray ionization of H2 leads to the production of a photon with a frequency 
in the range v — > v + dv, and £h 2 is the cosmic ray ionization rate of H2. 

For gas at the center of a spherically symmetric protogalactic core, we can simplify Equation !42l to 

Rx= / e(v,r)e- T( "' r) dr dv. (44) 



Jo Jo 

where t(v, r) is the optical depth between the center of the halo and gas at a radius r, e(y, r) is the emissivity at r, and R 
is the core radius. However, even after making this simplification, calculation of Rx still requires more information than we 
have available in our one-zone calculation, namely the radial profiles of density, temperature and chemical abundances, which 
determine both r(v,r) and e(v,r). In their absence, we are forced to approximate. 

If we assume that the protogalactic core has a density structure with a 'core plus halo' form, i.e., 

n ( r ) = <i 1 1 w . (45) 
l_ n c (r c /r) a , r > r c 

then provided that a > 1, the integral in Equation [42] will be dominated by the contribution from the core of the density profile. 
If we further assume that the core is chemically homogeneous, and that r(v,r) ~ 0, then we can approximate Equation 1421 as 

poo 

Rx~ I crx(^)r c PH 2 (^)CH 2 n c ,H 2 dv, (46) 
Jo 

where n c ,H 2 is the number density of H2 within the core. Note that even if the point we are considering is not directly at the 
center of the core, Equation [46] remains a reasonable approximation to Rx, provided that we are considering a point within 
r c , and that r c <C R. Provided that our approximations hold, Equation 1461 allows us to reduce what is formally a non-local 
problem into one that can be treated as if it were local. 

To properly include the effects of hydrogen excitation, one would have to solve for the radiative transfer of the Lyman 
alpha photons within the collapsing protostellar core. However, as the outcome would be highly sensitive to the assumed 
density and velocity profiles of the gas, which are not available from our one-zone calculation, the wisdom of performing 
such a detailed calculation for each of our simulations that include cosmic rays is questionable; we run the risk of getting an 
answer that is completely determined by our assumptions, and that therefore is not robust. Instead, we have chosen a more 
conservative course, and have attempted to put limits on the effects of the Lyman alpha photons by considering two limiting 
cases: one in which they do not propagate significantly into the core of the protogalaxy, and do not contribute to Rx (which is 
thus given in this case by Equation 1461 above), and another in which the optical depth of the gas to the Lyman alpha photons 
is negligible, and Rx is given by a generalization of Equation 1461 



Rx~ ox(v)r c (PH 2 (^)CH 2 n c ,H 2 + Pn(v)Cnn c ^) dv, (47) 
Jo 

where Pn(v)dv is the probability that the cosmic ray ionization of H leads to the production of a photon with a frequency 
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v + Av, is the cosmic ray ionization rate of H, and n c ,H is the number density of atomic hydrogen in the 

assume, following lOmukai (2000) and 



in the range v 
core 

To evaluate Rx , it remains necessary to specify r c . In our calculations 
Omukai et all |2005h . that r c is given by the current Jeans length. 

In Table |A"l4l we list estimated values for crx.cff,H2 — Jq cr x(^ y )-PH2 iy) d/^ and 0"x,cff,H — Jq o"x(^).Ph(^) dv for both 
photoionization and photodissociation for a number of different chemical species. To compute these values, we assumed that 
all of the photons produced by excited hydrogen are emitted in the Lyman-a line, implying that Ph = 5(v a — v), where v a is 
the freq uency of Lyman-a and 5 is the Dirac delta function. For Pn 2 , we used estimated values based on the emission spectra 
given in lSternberg. Dalgarno. fc Leppl (|l987l ); note that these are likely accurate only to within a factor of a few. Given these 



values, the rates for the cosmic-ray induced photoionization and photodissociation of these species in our model cores can be 
calculated using Equation 1461 or l47l as appropriate. 



3.4 Neglected processes 

Although the chemical network presented in this paper is, to the best of our knowledge, the most comprehensive network used 
to date to simulate primordial gas chemistry, there remain a large number of possible reactions that we have not included. 
Below, we discuss which types of processes have been omitted, and why. 

(i) We do not include the formation of H2 or D2 by the radiative association of ground state atomic hydrogen or deuterium, 
respectively, on the grounds that the rate coefficients for these processes are negligible. 



(ii) We have not included reactions th at involve electronically excited atomic hydrogen (as considered in lLatter fc Black 



1991 or Rawli ngs. Drew, fc Barlowlll993l . for instance). We justify this omission by noting that the population of the n — 2 



electronic level of atomic hydrogen will be very small at the densities and temperatures considered in this work, on account 
of the large Einstein coefficient associated with the Lyman-a transition and the large energy separation of the n = 1 and 
n = 2 levels. For similar reasons, we have not included any reactions that require electronically excited deuterium, helium or 
lithium. 

(iii) We have restricted the range of chemical species considered to those with three or fewer atoms. In principle, the 
formation of larg er species is possible - for instance, Hg" can be formed from H^" by the radiative association reaction 



( Paul et al.Ml995T ) 



H+ + H 2 -» H+ + 7. (48) 

However, the chemical abundances of the three-atom species in our model are very small, and we expect the abundance of 
even larger species to be much smaller still. It therefore seems unlikely that they will play any significant role in the cooling 
or chemistry of the gas. 

(iv) We have omitted any photochemical reactions that require photons with energies greater than 13.6 eV, under the 
assumption that any such photons emitted by external sources of radiation will be absorbed in the intergalactic medium, or 
in the interstellar medium of the protogalaxy, before reaching the particular collapsing core under study. Moreover, since we 
consider only the initial collapse of the core, internal sources of radiation, such as a central protostar, fall outside of the scope 
of this work 

(v) We do not include processes that have negligible reaction rates at all temperatures treated in this work. This includes 
the prod uction of doubly-ionized heliuim He 2+ , or doubly or triply ionized lithium, Li 2+ and Li 3+ by collisional ionization 



(see e.g., Lepp. Stancil. fc Dalgarno! 2002 ). which are therefore omitted from the chemical model. 

(vi) We have ignored the effects of stimulated radiative association and stimulated radiative attachment, i.e., reactions of 
the form 

X + Y + 7b XY + 7 + 7b , (49) 

X + e~+7 b -> X~+7 + 7b, (50) 

where 7b represents a background photon. The influence of stimulated radiative association or attachment on the producti on 
of various molecules in primordial gas (Li H, HD, H~, Li - and HeH + ) has been investigated I Stancil fc Dalgarno! 1997a 



Il998l : Izygelman. Stancil. fc Dalgarno! 1993). generally for the case of a black-body radiation field. However, significant effects 



are found only for radiation temperatures T ra( j > 500 K, much larger than the CMB temperature at the redshifts of interest in 
this study. The background radiation fields considered in £|5.4l have the same shape below hv = 13.6 eV as a 10 K black-body, 
but have intensities that are orders of magnitude weaker than a true black-body radiation field with this temperature, and so 
are also unimportant in this context. Therefore, it is clear that the influence of these simulated processes will be negligible, 
(vii) We have omitted collisional processes such as 

H2 + H2 ^H + H + H + H, (51) 
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or 

H + H -> H+ +H+ + e" +e", (52) 

that have energy thresholds corresponding to temperatures significantly higher than those considered in this work, as we do 
not expect these processes to play an important role in low temperature gas. 

(viii) We do not include dissociative charge transfer reactions involving or its isotopologues, e.g. 

H + H+ -> H + H + H+, (53) 
H + HD+ -> H + D + H+, (54) 

because at the temperature of interest in this study, the cross-sections for these proces ses are far smaller than those for the 
equivalent non-dissociative charge transfer reactions (IKrst ic 2002; Krstic fc Janevll2003l ) 

H + Hj -> H 2 + H+, (55) 
H + HD+ -f HD + H+. (56) 

(ix) We have not included charge transfer from He + to Li 

He + +Li^ He + Li+, (57) 
or its inverse 

He + Li+ -> He+ + Li. (58) 

The first of these reactions is unimportant in comparison to charge transfer from H + owing to the low He + abundances we 
find in our simulations. The second reaction is negligible at T < 10000 K due to its large endothermicity. 

(x) We have omitted all collisional dissociation reactions caused by minor molecules or ions, e.g. HD, Li, LiH, etc. Collisional 
dissociation reactions involving HD, such as 

H 2 + HD -> H + H + HD, (59) 
will be unimportant compared to the analogous reactions involving H2, while reactions of the form 

XY + Li -> X + Y + Li, (60) 

will be unimportant due to the very small abundance of lithium relative to hydrogen. 

(xi) We have ignored a large number of possible three-body processes: specifically, every process that involves any species 
other than H, H2 or He as the third body. At the densities at which three-body reactions become significant, the abundances 
of these three species are orders of magnitude higher than the abundances of any other species, and so it is easy to justify the 
omission of these minor contributions. 

(xii) We do not include transfer reactions involving two ions of the same charge, e.g. 

H+ + H+^H++H+ (61) 

as the mutual Coulomb repulsion of the ions renders these reactions ineffective at the temperatures considered in this work. 

(xiii) We do not include the double ionization of H2 or He by cosmic rays, i.e. 

H 2 + C.R. -» H + +H + +e _ +e- (62) 
Hc + C.R. -» He ++ +e-+e", (63) 



as th e fraction of cosmic ray ionization events leading to these outcomes is expected to be very small IjGlassgold fc Langer 
19731 ). 



(xiv) We assume that in mutual neutralization reactions involving H J or one of its isotopologu e s, com plete breakup of the 
ion is unlikely to occur; this is in line with e.g. the detailed chemical network of Le Teuff et al. ( 2000f ) . which includes the 
processes 

Et + H~ -> H 2 + H 2 (64) 
and 

Et + H~ -> H 2 + H + H, (65) 
but not 

Et + E" ->H + H + H + H. (66) 

(xv) A number of possible reactions involving LiH J have been omitted, as have any reactions involving LiHD + or LiDj. 
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Table 2. Processes included in our thermal model 



Species 


Process 


Collision partner (s) 


Refs. 


Cooling: 








H 2 


Rovibrational lines 


H, H 2 , He, H+, e" 


1, 2 


H 2 


Collision-induced emission 


H 2 


2 


HD 


rovibrational lines 


H 


3 


LiH 


rovibrational lines 


H, H 2 


1 




rovibrational lines 


H, H 2 , He, e- 


5 


H 


resonance lines 


c~ 


6 


CMB photons 


Compton scattering 


c~ 


6 


H+, HD+, D+ 


Rovibrational lines 


H, e" 


7 


Minor species 


Rovibrational lines 


H, H 2 , He, e" 


8 


H+ 


Recombination 


c~ 


9 


He+ 


Recombination 


e — 


10 


H 


Collisional ionization 


c~ 


11 


H 2 


Collisional dissociation 


H, H 2 , He, e- 


12 


H+ 


Charge transfer (reaction CT19) 


H 2 


13 



Heating: 

H 2 Formation 12 

H 2 Photodissociation 14 

H 2 Ultraviolet pumping 15 

Cosmic rays Ionization/excitation 16 



Note: See the ap propriate subsections in j|4"lfo r det ails of how we have d ecided whic h collision partners to include. 

References: 1 - IWrathmall fc Flowerl d2007h: 2 -IRipamonti fc Abeil d2004h : 3 - iLipovka. Nunez-Lopez fc Avila-Reesej j2005h : 4 - 

pel h004): 6 - iBlackl lll98ll) , ICenl Jl992h ; 7 - See g5l 8 - Se e 



iLcicieiiLca. _l — v vi enjiniieiii oc i iuwui \ z/uu i a ~ iviuciiiiuiiii oo nuui 

iGalli fc Pallal jl998h: 5 -iNeale, Miller, fc Tennvsonl dl996l' ). lOka fc Epc' 



9 -iFerland et al] 1 11992 
iBlack fc Dalgarnol (fl977 



Hummer fc Storevl dl998h: 11 -[Janev et al 
Burton. Hollenbach, fc Tielenj dl990T) : 16 - 



Il987h: 12 - See M a nd 33*91 13 - ISavin et all l|2004h : 
Goldsmith fc Laneerl (l978) 
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As we have already discussed in i\3. 1.111 we know very little regarding the values of the rate coefficients for many of the 
most important formation and destruction processes for LiH J, rendering its abundance highly uncertain. In view of this large 
uncertainty, there is little to be gained by adding in additional, equally uncertain but less important processes involving LiH J , 
or by considering the chemistry of the deuterated forms of the molecular ion. 



4 THERMAL PROCESSES 

Our model of the thermal behaviour of the gas includes the effects of heating and cooling from a large number of different 
radiative and chemical processes. A full list of the processes included is given in Table [2j while a more detailed discussion is 
given below. 



4.1 H 2 cooling 



In our treatment of H 2 rovibrational cooling at low densiti es, we include the effects of collisions between H 2 and H, H 2 , He, H + 
and e~, using fitting formulae from Glover fc Abel (2008). At high densities, we use the standard LT E cooling function (see 
e.g., Hollenbach fc Mckeelll979l ). We assumed the usual value of 3:1 for the H 2 ortho-para ratio, which [Glover fc Abell (|2008h 
have demonstrated is a good approximation for all temperatures in the range in which H 2 cooling is important. Although the 
revised treatment of H 2 cooling presented by Glover fc Abel ( 20081 ) can make a significant difference to the thermal evolution 
of the gas in some circumstances - notably, in gas with a substantial fractional ionization - we do not expect it to have a 
significant impact on our current results, as at the densities of greatest interest in this paper, H 2 is well within the LTE 
cooling regime. 

An importan t source of inaccuracy at high densities is the treatment of the opacity of the H 2 emission lines. In our 
models, we follow IRipamonti fc Abell |2004h and model optically thick H 2 cooling with the expression 



A 



H 2 , thick 



A H2 ,thin x min [l, (n/no) ] 



(67) 



where Ah 2 .thick and AH 2 ,thin are the optically thick and optically thin cooling rates, respectively, n is the number density of 
hydrogen nuclei, and where no = 8 x 10 9 cm -3 and {3 = 0.45. This approximation works well for modelling gas at the center 
of a collapsing c ore, but is less accur ate when used to treat H 2 cooling in the surrounding envelope (N. Yoshida, private 
communication). Yoshida et al. (2006) present a more accurate approach based on the computation of escape probabilities 
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Table 3. Numerical coefficients used in our analytical fit to the H^j~ cooling rate 





LTE 
20 < T < 400 K 


LTE 

400 < T < 10000 K 


n -> 
20 < T < 10000 K 




-1.6583133 x 10 4 


9.5033824 x 10 3 


-7.9192725 


ai 


5.0808831 x 10 4 


-1.7832745 x 10 4 


-43.505799 


02 


-5.9475456 x 10 4 


1.2847118 x 10 4 


41.100652 


<*3 


2.8459331 x 10 4 


-3.9079919 x 10 3 


-17.327161 


0*4 


1.9988968 x 10 3 


-2.8286326 x 10 1 


3.3895649 


as 


-8.6370305 x 10 3 


3.7394515 x 10 2 


-0.24931287 


ae 


4.0429912 x 10 3 


-1.1130317 x 10 2 


0.0 


«7 


-8.2863818 x 10 2 


1.4187579 x 10 1 


0.0 


18 


6.5975582 x 10 1 


-6.9969136 x 10" 1 


0.0 



for each individual H2 line using the Sobolev approximation. However, this treatment requires dynamical information, in the 
form of the local velocity gradient, that is not availabl e in an y meaningful form in our simple models, and so for our current 
study we must be content with the iRipamonti fc Abell (|2004h approach. 

At very high densities (n > 10 14 cm -3 ), cooling from H2 becomes dominated by collision-induced emission (CIE). When 
an H2 molecule collides with a hydrogen or helium atom, or another H2 molecule, the particles involved briefly act as a 
'supermolecule' with a non-zero electric dipole, which has a high probability of emitting a photon. Because the collision time 
is very short, the resulting collision-induced emission lines are very broad, an d typically merge in t o a c ontinuum. In our model, 
we model the CIE cooling rate with a power- law approximation taken from IRipamonti fc Abell (|2004t ). valid for gas in which 
xh 2 > 0.5: 

Acie = 4.578 x lfr 49 T 4 n H2 n erg s" 1 cm" 3 . (68) 

Although cooling from collision-induced emission contributes only 12% of the total cooling at the density reached at the end 
of our simulation, ri{ = 3 X 10 13 cm" 3 , we have verified in test runs that at higher densities it very rapidly becomes the 
dominant form of cooling, justifying our decision to end our simulations at this point. 



4.2 H3 cooling 



In the LTE limit, it is straightforward to calculate the H3 cooling rate using the data presented by Neale. Miller, fc Tennyson 
1 1996h . In our simulations, we use tabulated values computed directly from the Neale. Miller, fc Tennyson! line list. However, 
for the convenience of readers, we also provide an analytical fit of the form 



log 



A 



HJ.LTE 



= ^a l (logT) 1 , 



(69) 



where A, 



is the cooling rate per unit volume in the LTE limit, and where the values of the a t coefficients are listed 



in Table [3] This fit is accurate to within 25% for temperatures in the ra nge 20 < T < 400 K and to within a few percent for 
temperatures in the range 400 < T < 10000 K. At high temperatures, the Neale. Miller, fc Tennyson ( 19961 ) line list is known 
to be incomplete, and so for T > 3000 K, we may systematically underestimate the cooling due to Hg~. However, as we shall 
see in fjSl we never find gas at these temperatures in the regime where H^" is potentially important, and so this incompleteness 
will not affect our results. 

At densities where H3" is not in LTE, the calculation of the H^ cooling rate presents more of a problem . A commonly 
used approximation for dealing with the cooling from molecular species (see e.g., Hollenbach fc Mckee 19791 ) is to compute 
the cooling at a density n using the expression 

A = Alte (70) 
l + (n cr /n)' [m> 

where Alte is the LTE cooling rate per unit volume, and where the critical density n CI is given by n cr /n = Alte/A„^o, 
where A n ^o is the cooling rate per unit volume in the n — > limit. For n <C n cr and n 2> n CI , this expression is highly 
accurate, while for n ~ n cr it does a reasonable job of capturing the basic behaviour, at a far smaller computational cost than 
a full level population calculation would require. We adopt this approximation in our treatment of Hjj" cooling, reducing the 
problem of calculating A H + to one of calculating A H + n ^ Q . Here, however, we hit a problem. To compute A H + n ^ , we must 
evaluate 



A H+,n->0 = n H+X! C ' (U - B(U 

3 



(71) 
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where Coj is the rate of collisional excitation from the HjJ" ground stat^], here denoted as level 0, to an excited level j, and 
Eoj is the energy difference between level and level j. The collisional excitation rate for transitions from — > j is simply 



Coj = qOj^nn + Qo.7,H 2 flH 2 + Q'Oj.He'lHe + 1oj,e 



(72) 



where qoj,H, Qoj,u 2 , Qoj,ne arid qoj, e - are the collisional excitation rate coefficients for collisions with H, H2, He and e~, 
respectively; tih, «-h 2 , RHe and n e - are the corresponding particle number densities; and where we have ignored the effect 
of collisions with protons (which are unimportant in the case of a positively charged ion such as H3") or with minor ionic or 
molecular species such as H~ or HD. The difficulty in computing the collisional terms, and hence the low density limit of the 



iit co oling rate, arises because most of the required collisional excitation rate coefficients are unknown. iFaure fc Tennyson 



[ 2003T I give rate coefficients for the collisional excitation of a number of low-lying rotational states by collisions with electrons, 
but in the high density, low ionization conditions of interest in this study, collisions with electrons are unimportant, and 
analogous datasets for collisions with H, He or H2 are not available. 



To deal with this problem, we have used an approach introduced bv lOka fc Eppl (|2004j). They computed rate coefficients 



for rotational transitions in Hg~ caused by collisions with H2 by making use of the principle of detailed balance and by assuming 
that the collisional transitions are completely random (i.e., that they obey no selection rules). These assumptions led them 
to suggest rate coefficients of the form 



Ej — Ei 
2kT 



(73) 



for transitions between an initial level i and final level j, where Qi and gj, are the statistical weights of levels i and j, 
respectively, Ei and Ej are the corresponding level energies, and Kij is a normalizing factor given by 



K i:i = C 



1+ E 



9n 



1/2 



exp 



E„ 



(1/2)(E S +Ei) 



2kT 



(74) 



where C i s the total collision rate, which is independent of i and j, and where the summation does not include levels i or j. 
Although 1 Oka fc Eppl (J200J) consider only pure rotational transitions, the same scheme can be used to treat ro- vibrational 
transitions. 

In our treatment, we assume that the Oka & Epp scheme can be used to treat collisions with atomic hydrogen and 
helium as well as H2, and hence are able to determine the temperature dependence of the set of collisional excitation rate 
coefficients for each collider (goj,H, <Zoj,H 2 and goj,He); we ignore collisions with electrons, on the grounds of the very small 
electron abundance that exists at the densities where Hj cooling is potentially important. Using these collisional excitation 
rate coefficients, we can then construct Coj via Equation 1721 from which A H + n ^ Q follows via Equation [FT] The overall 

normalization of the cooling rate remains uncertain, as it depends on the H3" number density, and on the total collision rates 
with each of H, H2 and He, which we can write as Ch, Ch 2 and Ch 6 - If we define the total collision rate C to be the sum of 
these three unknowns 



C — Ch + Ch 2 + CHe, 

and write the low density H3" cooling rate as 



A 



H+,n^0 



H 3 ' 



*o"h+> 



(75) 



(76) 



then it is easy to show that the combination i H + n ^ /C is completely determined. We have computed an analytical fit to this 
quantity, using a fit of the form of Equation 1691 The fitting coefficients are listed in Table [3] This fit is accurate to within 1% 
over the temperature range 20 < T < I0000K. To convert from L H + n _ >0 /C* to L H + n ^ , we must fix the size of our remaining 
free parameter, the total collision rate C. In most of our simulations, we assume that C is given by 
9 



C 



2.2 x I(T m ?ih + 1.9 x I0~ 9 n H , + 8.1 x IO~ 10 n Hc s^ 1 



(77) 



which is the sum of the Langevin rates for collisions betwee n Hj~ and H, H2 and He, respectively. These Langevin rates were 
computed using polarizabilities for He and H2 taken from Huiszoon fc Brielsl 1 19931 ); the exact value was used for H. The 
true value of C is unlikely to be very much larger than this, but could be considerably smaller, and so in £|5 . 2 1 we examine 
the sensitivity of our results to our choice of value for C. Finally, we note that as both A H + LTE and A H + n ^ Q are directly 

proportional to n„+ , the Hi" critical density is independent of n„+ . Thus, once C is specified, n cr can be trivially computed 

3 3 
using our numerical fits. 



1 The ( J, K) = (1, 1) rotational level of the vibrational ground state; occupation of the (J, K) = (0, 0) level is forbidden by the Pauli 
exclusion principle 
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We assume that the Hg emission remains optically thin throughout our simulations. For a subsonic collapse in which the 
effect of large-scale velocity gradients are unimportant compared to the local Dopple r broadening of the emis sion lines, the 
optical depth at line center corresponding to a given emission line can be written as IjMihalas fc Mihalaslll984h 

„2 



f/j 



-Aji- 



N, 



(78) 



where gj and gt are the statistical weights of levels i and j, Vij is the frequency of the transition from level j to level 
i, Aij is the corresponding spontaneous radiative transition rate, Ni is the column density of absorbers in level i, and 
A^d = (yij /c)(2kT/m) 1 ^ 2 is the Doppler width of the line, where m is the mass of the Hg" ion. For simplicity, we have neglected 
the effects of stim ulated emission. Illustrative value s for Vij and Aij for a strong vibrational transition are Uij ~ 8.46 x 10 13 Hz 



and Aij ~ 94 s 1 ( Neale. Miller, fc Tennyson 19961 ). and so for this transition 



lO" 15 ^. 



(79) 



where we have assumed a gas temperature of 1000 K. Therefore, in this example, the line becomes optically thick only once 
Ni > 10 15 cm -2 . As many of the Hg~ emission lines are considerably weaker than this example, and as the column density of 
Fit ions in any particular level i can clearly be no larger than the total Ht column density, N„+ , it is safe to conclude from 
this analysis that optical depth effects are unlikely to significantly affect the cooling rate until N H + > 10 15 cm -2 . 

We now investigate whether we expect our models to reach this column density in . Our one-zone dynamical model 
does not contain any information about the overall structure of the collapsing core and so does n ot predict AT H + directly. 
However, based on the results of more detailed numerical simulations ( Abel. Bryan, fc Norman 20021 : Yoshida et all 20061 ). we 
assume that the protostellar core has a density profile that is well approximated by a power law n(r) oc r~ 2 ' 2 , and that it is 
collapsing subsonically. With this assumed density profile, the column density of hydrogen nuclei along a radial ray from a 
point r to the edge of the core is given by 



n(r')dr' , 

\ ' core / 



(80) 
(81) 



where n(r) is the number density of hydrogen nuclei at r and r CO rc is the radius of the core. As we shall see in £0 the H^ 
abundance in the collapsing gas typically varies only slightly with density below some threshold density nthr and then declines 
sharply for n > nthr- The value of n t hr depends on factors such as the cosmic ray ionization rate and the speed of the collapse, 
but even in the most extreme models (e.g., run CR5; see i)5.3[) n t hr = 10 11 cm -3 , while in general it is much smaller. Therefore, 
almost all of the contribution to N„+ comes from gas at densities n < nthr, and hence at radii r > rthr, where r t hr is the 

3 

radius such that n(rthr) = n t hi- If we assume thar rthi <C r CO rc, or in other words that the density distribution of the collapsing 
core extends to densities n <C nthr, then Equation [81] tells us that the column density of hydrogen nuclei between r t hr and 
the edge of the core is approximately 



N to t(r t hv) ^ -rthrnthr 
6 



(82) 



Denoting the fractional abundance of HjJ" at rthr as s H +(r t hr), and assuming that it remains constant for r > rthr, we can 
therefore write the H^ column density between r t h r and the edge of the core as 



N n +(rtuA 



-x„+(r thl )r thl n thl 

6 3 



(83) 



For a typical protog alactic core, simulations have s hown that a density n t h r = 10 11 cm 3 corresponds to a radius 



rthr 



10 1 



15 cm (see e.g., Abel. Bryan, fc Norman 20021 ; Yoshida et al. 2006^ . Furthermore, the results presented in H5.3I 



demonstrate that in run CR5, 



(rthr 



4.4 x 10" 11 . We therefore obtain N +(r thl ) ~ 4 x 10 15 cm -2 . At r > r thr , the 



3 3 

Hg" column density is smaller, but at r < rthr, it does not grow significantly larger. In this particular case, the core may be 
marginally optically thick, albeit only in the strongest lines. However, in most of our models, n t hi is significantly smaller, and 
x„+ (r t hr) is orders of magnitude smaller. In these runs, it is clear that the core remains optically thin. 



4.3 HD cooling 

To model HD cooling, we use the cooling function of Lipovka. Nunez-Lopez, fc Avila-Reese 1 2005h . Although formally valid 
only in the temperature range 100 < T < 2xl0 4 K, we have compared its behaviour at lower tempe r ature s with an explicit cal- 
culation of the cooling rate made usi ng radiative de-excitation rates from Abgrall. Roueff, fc Viala (1 19821) and collisional rates 

extrapolated from those computed by Wrathmall. Gusdorf. fc Flower ( 2007h . We find that the Lipovka. Nunez-Lopez, fc Avila-Reesel 
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( 2005h rate remains reasonably accurate down to temperatures as low as 50 K, with error s no greater than 20%, and that even at 
T = 30K it remains accurate to within a factor of two. At temperatures T ^> 100K, the Lipovka. Nunez-Lopez, fc Avila-Reesei 



2005) cooling rate slightly underestimates the effects of HD cooling compared to the newer calculations of lWrathmall. Gusdorf. fc Flower 



2007]), presumably owing to the more accurate vibrational excitation rates used in the latter study, but the differences are 



never great er than about 50%, and in any case occur in the temperature regime in which H2 cooling dominates. The break- 
down of the Lipovka. Nunez-Lopez, fc Avila-Reese! 1 20051 ) fit at very high temperatures (T > 20000 K) is unimportant, as the 
gas in our models never reaches this temperature. 

To correctly model the effects of HD cooling at low temperatures, it is necessary to take the effects of the cosmic microwave 
background into account. We do this approximately, by using a modified HD cooling rate, Ah D , defined as 



Ahd = Ahd(T) — Ahd(Tcmb) 



(84) 



where Ahd^) and Ahd(Tcmb) are the unmodified HD cooling rates at the gas temperature T and the CMB temperature 
Tcmb, respectively. 

The quoted range of densities for which the Lipovka. Nunez-Lopez, fc Avila-Reese 1 20051 ) cooling function is valid is 
1 < n < fO 8 cm -3 . To extend the range of the cooling function to densities n < f cm -3 , we assume that at these densities the 
HD cooling rate scales proportionately to the number density of HD times the number density of colliders (i.e. as n 2 ), and 
hence that 



A H D(n = n) = (n') 2 A H D(n = 1) 



(85) 



for n' ^ 1 cm -3 , where Ahd(h-) is the HD cooling rate per unit volume (with units erg cm -3 s _1 ) at gas number density n. 
To extend the cooling function to high densities, n > fO 8 cm -3 , we assume that the HD molecule is in LTE and hence that 
the HD cooling rate per unit volume is independent of the number density of colliders and scales linearly with n; or in other 
words, that Ahd/^hd is independent of n. In view of the fact that I <C ii C i,hd <§C 10 8 cm~ 3 , where i! ci .hd is the HD critical 
density, both of these assumptions appear well justified. 



The Lipovka et al. (2005) cooling function only includes the effects of collisions between HD and H. However. iFlower et al 



(2000) 'rave shown that the influence of the H2/H ratio on the HD cooling rate is very small, and so the Lipovka et al. (2005) 
cooling function should remain reasonably accurate even after the molecular fraction becomes large. Moreover, collisions 
between HD and other species (electrons, H + , etc.) can be neglected compared to collisions with H on account of the much 
larger abundance of the latter at the gas densities of interest. 

Finally, we assume that the HD rovibrational lines remain optically thin throughout all of our runs. In practice, this is 
probably not the case: the strongest HD lines will become saturated once the HD column density exceeds JVhd = 10 22 cm -2 , 
and an analysis similar to that performed for H^ in the previous section suggests that this will occur in our model cores once 
n > 3 x 10 13 cm -3 . However, HD is only a minor coolant at these high densities, and so we can safely neglect optical depth 
effects on HD cooling without significantly affecting the thermal evolution of the gas in our models. 



4.4 LiH cooling 

To treat cooling from LiH, we use the cooling function given in Galli fc Palla (1998). This is the low density limit of the 
LiH cooling rate and so is strictly valid only for gas densities significantly below the LiH critical density, n C r,LiH- However, 
the transition probabilities for the rotational and vibrational transitions of LiH are very lar ge, on account of th e molecule's 



large dipole momen t. This means that the LiH critical density is large, n cr ,LiH — 10 12 cm 3 ( Lepp fc Shulll 1984 ). and so the 



Galli fc Pallal |l99Sl ) cooling function is a reasonable choice over most of the range of densities covered by our simulations. At 
very high densities, we would expect LiH to begin to reach LTE, and our approximation to break down; at these densities, our 
continued use of the Galli fc Palla ( 19981 ) cooling function means that we will overestimate the effectiveness of LiH cooling. 
Despite this, we find LiH cooling to be ineffective at all densities (see i)5. II below), suggesting that if we were to use a more 
accurate treatment of LiH cooling at high densities it would not significantly alter our conclusions. 

In principle, we should adjust the LiH cooling rate to account for the effects of the CMB, just as we do for the HD cooling 
rate. However, as LiH cooling proves to be unimportant at all densities, this correction is also unimportant and its omission 
does not significantly affect the thermal evolution of the gas. 

We assume that the LiD cooling rate is the same as the LiH cooling rate. While this is a crude approximation, in practice 
the LiD abundance is so small that its contribution is always negligible and thus the choice of LiD cooling rate is unimportant. 



4.5 H+, HD+ and D+ cooling 



To treat cooling from these molecular ions, we use the same approach as in Glover fc Abell ( 20081 ) . At low densities, most 



Hj cooling occurs due to collisions wi th free electrons and neutr al hydrogen atoms; collisions with He and H2 excite Hj at 



comparable rates to collisions with H (jRoberge fc Dalgarnolfli&j ), but are unimportant due to the low abundances of these 
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species relati ve to atomic hydrogen. To m odel the cooling due to collisions with electrons, we use the vi brational excitation rate 
coeffic ients of Sarpal fc Tennyson 1 19931 ). while for collisions with neutral hydrogen, we use a fit to the Suchkov fc Shchekinov 
1 197fih rate coefficient provided to us by D. Galli (private communication); note that this is a factor of ten smaller than the 
rate coefficient given in iGalli fc Pallal (1998), owing to a normalization error in the latter paper. 

At high densities, the vibrational levels of will be in LTE. In this regime, the cooling rate is given approximately by 



A 



r+ TTF , = 2.0 x 10~ 19 T 01 expi 



3125\ 



This fit is from Glover fc Abell ( 20081 ) and includes contributions from all vibrational states v 8 (higher vibrational states are 
not expected to cont ribute significantly at the temperatu res of interest in this work) . It wa s computed using level energies from 
Karr fc Hilicol ( 20061 ) and radiative transition rates from Posen. Dalgarno. fc Peek ( 19831 ). The effects of rotational excitation 
were not included, but are unlikely to change this expression by a large amount, owing to the very small Einstein coefficients 
associated with pure rotational transitions in H^". 



At intermediate densities, we assume that the HJ vibrational cooling rate is given approximately by the function 



V, 



1 + n cr / r 



(87) 



where n CI /n — A H + LTE /A H + n _^ is the critical density, and where A H + n _ >0 is the cooling rate in the low density 



limit. This is given by 



A H+,n-»0- L H+.c n o 



±l n , rl xl„ 



(88) 



where L„+ 



and L, 



are the cooling rates per H^" ion per unit collider density for collisions with electrons and atomic 



hydrogen, respectively, taken from Sarpal fc Tennyson ( 19931 ) and Suchkov fc Shchekinov 1 19781 ) as noted above. 

To model cooling from vibrational transitions in HD + , we assume, in the absence of better information, that the low 
density cooling rate is the same as that used for H J . However, since HD + has much larger radiative transition rates than H J , 
the LTE cooling rate for HD + is much larger than that for H^~. Accordingly, we use the following functional fit for the HD + 
LTE cooling rate: 



1 exp 



A H d+,lte = 1.09 x 10 
at temperatures T ^ 1000 K and 
A 



HD+ ,LTE 



5.07 x 10 _12 T 014 exp | 



2750 \ 

—jT I n HD + 



2750 \ 

—jT J %D+ 



(89) 



(90) 



at T > 1000 K. These fits are from Glover fc Abell 1 20081) and were calculated u sing HD + level energies from Karr fc Hilico 
1 20061 ) and radiative transition rates from Peek. Hashemi- Attar, fc Beckel ( 19791 ). For densities between the low density and 
LTE limits, we again use a function of the form of Equation [87] to compute the cooling rate. 

Finally, to model cooling, we simply assume that the same rates apply as for Hj cooling. In practice, the very small 
size of the typical abundance renders this process completely unimportant. 



4.6 Other radiative coolants 

In addition to the coolants discussed above, we also include a treatment of cooling from a number of other minor molecular 
ions that are present in the gas. Specifically, we include the effects of cooling from H2D + , HD^~, Djj~, HeH + , HeD + , HeJ , 
LiH + , LiD + and LiHjj". Since appropriate collisional data are not available for most of these species, we treat their contribution 
to the cooling rate in an extremely simple fashion. We assume that the contribution to the cooling rate made by a species i 
with number density m can be written as 

A* = kT dcnA m, (91) 

where n c is the number density of a collider c, C; c is the rate coefficient for inelastic collisions between i and c, and where we 
sum over all possible colliders. For collisions with H, H2 or He, we assume that C; c is given by the Langevin rate, while for 
collisions with electrons we conservativ ely assume that Cj c = 10 ~ 6 cm 3 s" 1 , which is comparable to the total rate coefficients 
found for other molecules, such as H^" jFaure fc Tennvsoni r2003). Collisions with all other species can be and are neglected. 



In constructing this approximation we have assumed that each collision with i transfers an amount of energy kT, all of 
which is subsequently radiated. In practice, this procedure is likely to significantly overestimate the cooling provided by i, 
for several reasons. For one thing, it is not clear that the mean amount of energy transferred in a collision will always be 
~ kT, since collisions that transfer AE <<iC kT are possible while collisions that transfer AE 3> kT are highly unlikely. More 
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importantly, this procedure neglects effects such as the collisional de-excitation of excited levels that will significantly limit 
cooling at high densities. However, these simplifications are unlikely to significantly affect our results, since even when we 
use these overestimates for the cooling rates, we find that the contribution of these minor species to the total cooling rate is 
negligible (see fj5|. 

As well as these minor coolants, we also include two forms of cooling that are of great importance in hot, ionized gas. 
The first is cooling from electron impact exci tation o f atom ic hydrogen (Lyman-a cooling) . We treat this using a rate from 
Cer] (|l992h . which is itself based on a rate in iBlackl (| 198 IT) . However, we note that for temperatures T < 8000 K, Lyman-a 



cooling is completely negligible, and so it does not significantly affect the outcome of the simulations presented in this paper. 

The se cond proces s is cooling due to the Compton scattering of CMB photons by free electrons. This is also treated using 
a rate from Cen (.1992), but again plays very little role in the thermal evolution of the gas, since it is important primarily at 
low densities (n < 1 cm" 3 ), even when the gas is initially highly ionized. 

Finally, we note that we do not include the effects of cooling from D2. As D2 is a homonuclear molecule, it suffers from 
the same drawbacks that H2 does with regard to low temperature cooling. However, as the results in £|5 - 1 1 demonstrate, it 
generally has a chemical abundance that is many orders of magnitude smaller than H2. Thus, in contrast to HD, it appears 
highly unlikely that D2 cooling is ever significant. 



4.7 Radiative heating 



We include two forms of radiative heating that can be significant if a strong ultraviolet background is pres ent. The first is the 
photo dissociation of H2. We calculate the H2 photodissociation rate as discussed in [J3] and then, following iBlack fc Dalgarno 
(1977]), we assume that each photodissociation deposits 0.4 eV of heat into the gas. Note that although the photodissociation 
of other ionic and molecular species (e.g., H^~, HD) will also heat the gas, they are unimportant when compared to H2 
photodissociation owing to the low abundances of the other species relative to H2. 

The second form of radiative heating included in our thermal model arises due to the population of excited vibrational 
states of H2 produced by radiative pumping by the UV field. At high densities, this leads to heating of the gas, as most of 
the excited molecules u ndergo collisional de-exc itation. We adopt a radiative pumping rate that is 8.5 times larger than the 
photodissociation rate ( Dmine_& i 13erto i ldi Il996 ) , and assume that each excitation transfers an averag e of 2(1 + n CI /n)~ 1 eV 
to the gas ( Burton. Hollenbach. fc Tielenslll99C ), where n cr is the H2 critical density, calculated as discussed in §3] above. 



4.8 Cosmic ray heating 

Following Goldsmith fc Langer 1 19781 ). we assume that each primary ionization deposits 20 eV of energy into the gas, giving 
us a heating rate 

c 



: 3.2 x 10~ 



(lO-^s- 1 ) 



-1 -3 
n ergs s cm , 



.10-^' S" 

where ( = J^. ( t and we sum over all species listed in Table lA"l3l 



(92) 



4.9 Chemical heating and cooling 

Any exothermic chemical reaction will potentially heat the gas, while any endothermic reaction will cool it. In practice, 
however, the effect of most reactions on the gas temperature is small, and only in a few cases do we need to take chemical 
heating or cooling into account. 

In highly ionized gas, cooling due to the recombination of hydrogen and helium can be a significant effect, particularly 
at temperatures which are too low for Lyman-a cooling to be effective. However, recombination cooling becomes ineffective 
once the fractional ionization of the gas falls below x ~ 0.01, and it therefore plays no role at the high gas densities of interest 
in this study. 

Other forms of chemical cooling included in our model occur due to the collisional dissociation of H2 (reactions CD9, 
CD10, CD11 and CD12), the destruction of H2 by charge transfer with H + (reaction CT2), and the collisional ionization of 
hydrogen and helium (reactions CI1 & CI2). Cooling from these processes may be of some importance at very early times 
in simulations starting at high temperatures (T > 10 4 K), but in general the gas temperature is too low for these sources of 
cooling to be significant. 

As far as chemical heating is concerned, the most significant process is H2 formation heating. When H2 is formed by 
reaction AD1 or reaction CT1, it preferentially forms in an excited vibrational state, with an energy comparable to the 
exothermicity of the reaction (3.73 eV for reaction AD1, 1.83 eV for reaction CT1). In low density gas, this energy is simply 
radiated away, but for n > 10 4 cm -3 , most is instead converted into thermal energy by collisional de-excitation of the newly 
formed H2. H2 formation via H~ or Hjj" therefore acts as a minor heat source in gas with n > 10 4 cm -3 . 
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Three-body formation of H2 also heats the gas, since the third body in the collision generally carries away additional 
energy equal to the binding energy of the new H2 molecule, 4.48 eV. In our reference simulation, this is the dominant heat 
source for densities 5 x 10 10 < n < 2 x 10 12 cm -3 . 



5 RESULTS 

Our simple one-zone dynamical model of gravitationally collapsing primordial gas contains a number of free parameters. To 
fully explore the role of Hjj" cooling and its sensitivity to these free parameters, it is necessary to perform a large number 
of calculations. However, discussion of the results of all of these calculations to the same level of detail would not only be 
extremely tedious, but would also run the risk of obscuring our main results. Therefore, we proceed by first discussing in 
detail in i )5. li the results of a single calculation - our reference model, hereafter denoted as computational run REF - before 
highlighting in the subsequent sections the differences in outcome (if any) that result from alterations in our free parameters. 
Full details of all of the runs discussed here can be found in Table [4] 



5.1 The role of H3 cooling 

We begin our study by investigating the outcome of our reference calculation, run REF, whose parameters are indicated in 
Table [4] In Figure [1] we show how the fractional abundances of 28 of our 30 chemical species vary with density during the 
course of the collapse. For clarity, we have divided these species into four sets on the basis of the elements that they contain, 
and illustrate the evolution of each set separately in Figures [TJi-[l]l. The two species that are not plotted - He + and HeJ - 
have abundances that remain negligibly small throughout the calculation. 

Figure [T] demonstrates that the evolution of the H J abundance passes through four distinct phases. In the first phase, 

at n < 10 s cm -3 , the ratio of Ht to H2 remains approximately constant: i h +/j;h, ~ 10 -9 . This is a consequence of the 

3 

balance between two main processes: the formation of Hj by the radiative association of H2 and H + (reaction RA18) and its 
destruction by dissociative recombination (reactions DR4 and DR5). If we assume that these reactions dominate the formation 
and destruction of H^~ , then the corresponding equilibrium abundance of H3" is given by 

_ fcR A i8n H+ a:H 2 , , 

3 (fc D R4 + KDR5jn c 

which reduces to 

&RA18 1n -9 /„,% 

3 (fe D R4 + KOR5) 

if x c = a; H +, which is a very good approximation at these densities, as Figure [TJi demonstrates. 

The second phase in the evolution of the H^" abundance begins at a density of around 10 s cm -3 , when there is a sudden 
decrease in the H^" abundance. This decline is caused by the fact that at these densities, dissociative recombination is no 
longer the only significant destruction mechanism. Increasingly, H3" is also destroyed by reaction TR17: 

H+ + H^H++H 2 . (95) 

Although this reaction is endothermic, the temperature of the gas at these densities (T>800K; see Figure[2]) is high enough to 
make this mechanism significant in comparison to reactions RA18, DR4 and DR5, owing to the very low fractional ionization 
of the gas. 

The third phase occurs at n ~ 10 10 cm -3 as the decline in the H^ abundance is briefly halted by an increase in the Hj 
formation rate. This is caused by the increase in the H2 abundance at these densities, which itself is driven by the onset of 
efficient three-body formation of H2. 

Finally, at a density n > 10 11 cm -3 , the H3" abundance decreases once more, owing to the rapid loss of the few remaining 
free H + ions from the gas, and the consequent disruption of the major H^ formation mechanisms. The reaction responsible 
for this loss of H + ions is the same as the reaction driving the formation of Hjj~, namely RA18. If this were the only process 
operating, then it would convert all of the H + in the gas into Hjj~ in a time t conv , given approximately by 

tconv ^ "J • (^^) 

KRAI8WH2 

Comparing this timescale with the free-fall timescale, ta ~ 1.4 x 10 15 n~ 1//2 s, and taking Arais = 10~ 16 cm 3 s _1 , we find that 
tconv < tg if n > bO/x^. For xh 2 ~ 10~ 3 , this gives a critical density n CO nv ~5x 10 7 cm -3 . Therefore, in the absence of any 
other effects, conversion of H + to Hjj~ should occur rapidly once the number density exceeds n C onv- In practice, however, a 
second effect intervenes. The steady increase in the gas temperature at these densities soon results in reaction TR17 becoming 
a major destruction mechanism for H J , as noted above. Destruction of H3" by reaction TR17 produces H J ions, most of which 
are then destroyed by reaction CT3 
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Table 4. List of simulations run 



Run 


ni (cm 3 ) 


I] (K) x H+ 


■1-21 


Ch (s ) 


C/C Ic f 


'/ 


Notes 


REF 


1.0 


1000 


2.2 x 10~ 4 


0.0 


0.0 


1.0 


1.0 




CI 


1.0 


1000 


2.2 x 10" 4 


0.0 


0.0 


0.1 


1.0 




C2 


1.0 


1000 


2.2 x 10~ 4 


0.0 


0.0 


10.0 


1.0 




CR1 


1.0 


1000 


2.2 x 10~ 4 


0.0 


ID" 20 


1.0 


1.0 




CR2 


1.0 


1000 


2.2 x 10~ 4 


0.0 


10 -19 


1.0 


1.0 




CR3 


1.0 


1000 


2.2 x 10" 4 


0.0 


10 -18 


1.0 


1.0 




CR4 


1.0 


1000 


2.2 x 10" 4 


0.0 


io- 17 


1.0 


1.0 




CR5 


1.0 


1000 


2.2 x 10" 4 


0.0 


io- 10 


1.0 


1.0 




CR6 


1.0 


1000 


2.2 x 10~ 4 


0.0 


10 -16 


1.0 


1.0 


No Ht cooling 


CR7 


1.0 


1000 


2.2 x 10" 4 


0.0 


IO" 20 


1.0 


1.0 


No PT mechanism 


CR8 


1.0 


1000 


2.2 x 10~ 4 


0.0 


ltr 18 


1.0 


1.0 


No PT mechanism 


CR9 


1.0 


1000 


2.2 x 10~ 4 


0.0 


io- 10 


1.0 


1.0 


No PT mechanism 


CR10 


1.0 


1000 


2.2 x 10" 4 


0.0 


IO" 20 


1.0 


1.0 


'Maximal' PT mechanism 


CR11 


1.0 


1000 


2.2 x 10" 4 


0.0 


io- 18 


1.0 


1.0 


'Maximal' PT mechanism 


CR12 


1.0 


1000 


2.2 x 10" 4 


0.0 


io- 10 


1.0 


1.0 


'Maximal' PT mechanism 


UV1 


1.0 


1000 


2.2 x 10~ 4 


IO" 4 


0.0 


1.0 


1.0 


Optically thin 


UV2 


1.0 


1000 


2.2 x IO" 4 


IO -2 


0.0 


1.0 


1.0 


Optically thin 


UV3 


1.0 


1000 


2.2 x IO" 4 


1.0 


0.0 


1.0 


1.0 


Optically thin 


UV4 


1.0 


1000 


2.2 x IO" 4 


IO" 4 


0.0 


1.0 


1.0 


loll M ... ^ = Tel-, tTT~) ^ o 


UV5 


1.0 


1000 


2.2 x 10" 4 


IO -2 


0.0 


1.0 


1.0 


tah H „ — fcVi TTF> — 


UV6 


1.0 


1000 


2.2 x 10" 4 


1.0 


0.0 


1.0 


1.0 


fah H ,-. — Tsh TTn — 


Nl 


0.03 


1000 


2.2 x IO" 4 


0.0 


0.0 


1.0 


1.0 




N2 


30 


1000 


2.2 x 10 -4 


0.0 


0.0 


1.0 


1.0 




Tl 


1.0 


100 


2.2 X 10" 4 


0.0 


0.0 


1.0 


1.0 




T2 


1.0 


1000C 


2.2 X 10~ 4 


0.0 


0.0 


1.0 


1.0 




XI 


1.0 


1000 


10~ 6 


0.0 


0.0 


1.0 


1.0 




X2 


1.0 


1000 


10- 2 


0.0 


0.0 


1.0 


1.0 




X3 


1.0 


1000 


1.0 


0.0 


0.0 


1.0 


1.0 




ELI 


1.0 


1000 


2.2 x 10~ 4 


0.0 


0.0 


1.0 


1.0 


No D 


EL2 


1.0 


1000 


2.2 x IO" 4 


0.0 


0.0 


1.0 


1.0 


No Li 


EL3 


1.0 


1000 


o o * j in — 4 

2.2 X 10 


0.0 


0.0 


1.0 


1.0 


No D or Li 


RA 


1.0 


1000 


2.2 x 10~ 4 


0.0 


0.0 


1.0 


1.0 


^RA18 from ref. 1 


AR1 


1.0 


1000 


2.2 x 10~ 4 


0.0 


0.0 


1.0 


1.0 


See ij5.6.3| 


AR2 


1.0 


1000 


2.2 x 10~ 4 


0.0 


0.0 


1.0 


1.0 


See 95.6.31 


3B1 


1.0 


1000 


2.2 x IO" 4 


0.0 


0.0 


1.0 


1.0 


/ctbi from ref. 2 


3B2 


1.0 


1000 


2.2 x 10~ 4 


0.0 


0.0 


1.0 


1.0 


fcxBl from ref. 3 


3B3 


1.0 


1000 


2.2 x 10~ 4 


0.0 


0.0 


1.0 


1.0 


&TB2 f rom re f- 4 


3B4 


1.0 


1000 


2.2 x IO" 4 


0.0 


0.0 


1.0 


1.0 


&TB2 f rom re f- 3 


LP1 


1.0 


1000 


2.2 x IO" 4 


0.0 


0.0 


1.0 


1.0 


See f|5.6.4| 


LP 2 


1.0 


1000 


2.2 x 10~ 4 


0.0 


0.0 


1.0 


1.0 


See ^5.6. 4| 


LP3 


1.0 


1000 


2.2 x IO" 4 


0.0 


0.0 


1.0 


1.0 


See ij5.6.4| 


DYN1 


1.0 


1000 


2.2 x IO" 4 


0.0 


0.0 


1.0 


0.(5 




DYN2 


1.0 


1000 


2.2 x IO" 4 


0.0 


0.0 


1.0 


0.3 




DYN3 


1.0 


1000 


2.2 x IO" 4 


0.0 


0.0 


1.0 


0.1 
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HA + H — > H2 + H 



(97) 



producing H + ions. Therefore, most of the H + ions that are removed from the gas by reaction RA18 are replaced by this chain 
of reactions. A net loss of H + from the gas occurs only if the Hjj" ion produced by reaction RA18 is destroyed by dissociative 
recombination, rather than by reaction TR17. The proportion of the destroyed by dissociative recombination is given by 



(fcDR.4 + fcDRE 



(&DR4 + fc>R5)«.c + fcTR17TOI ' 



(98) 



where we have assumed that dissociative recombination (reactions DR4 and DR5) and reaction TR17 are the only significant 
processes destroying H^". Consequently, the net rate at which H + ions are removed from the gas is a factor /dr. slower than 
was assumed in our calculation of t conv above, and hence the actual timescale on which the majority of the H + ions are 
removed is given by 
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Figure 1. (a) Chemical evolution of the gas in our reference calculation. Fractional abundances are plotted for H2 (upper solid line), 
H^" (lower solid line), e _ (upper dashed line), H+ (lower dashed line), H (upper dash-dotted line), H^" (lower dash-dotted line) and 
H — (dotted line), (b) As (a), but showing the fractional abundances of He (solid line), HeH + (dashed line) and HeD + (dotted line). 
The abundances of He + and HeJ remained negligibly small throughout the simulation and are not plotted here, (c) As (a), but for 
the fractional abundances of HD (upper solid line), D (upper dashed line), D2 (upper dash-dotted line), D+ (central solid line), H2D+ 
(upper dotted line), D — (lower dash-dotted line), HD^~ (lower solid line), HD" 1 " (central dashed line), (lower dotted line) and 
(lower dashed line), (d) As (a), but for the fractional abundances of Li (upper solid line), Li+ (upper dashed line), LiH (upper dotted 
line), Li~ (upper dash-dotted line), LiH^~ (lower dotted line), LiH + (lower solid line), LiD (lower dash-dotted line) and LiD + (lower 
dashed line). 



floss — j. tconv ■ 
JDR 

Now, fi oss depends on the electron density through /dr, and so as long as H + is the dominant source of free electrons, decreasing 
its abundance increases fioss, thereby preventing rapid removal of the H + ions from the gas. However, once the H + abundance 
falls below ~ 10 , it is singly ionized lithium, Li + , that becomes the dominant positive ion. At this point, further decreases 
in x K + have very little effect on t\ OBS . Taking x e - = 10 -11 and T = 1000 K, and assuming that (fcoR.4 + &DRs)^e <S ^tr^^h, 
we find that /dr — 10 _4 Xy , and hence t\ OBB ~ 10 20 xn/nu 2 S. Thus, at the point at which Li + first becomes the dominant 
positive ion, which occurs around n ~ 10 9 cm -3 in our reference simulation, fi oss 3> tfj. However, the wholesale conversion of 
H to H2 by three-body reactions that begins to set in at around this density rapidly decreases fi oss , and at n ~ 10 11 cm -3 it 
becomes shorter than the free-fall timescale of the gas. At this point, most of the remaining H + ions are lost from the gas, 
following which H^~ formation largely ceases. Since the destruction of H3" by reactions DR4, DR5 and TR17 is unaffected by 
the fall-off in the H + abundance, the end result is a very rapid fall-off in the H^" abundance. 

To determine whether the small amount of HjJ~ that forms in the gas is enough to significantly affect its thermal evolution, 
we have compared the H;j" cooling rate per unit volume to the total cooling rate per unit volume. The results are plotted 
in Figure [3] We see that even at its moment of peak effectiveness, which occurs at n ~ 10 8 cm -3 , H^~ contributes no more 
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Log n (cm -3 ) 

Figure 2. Temperature evolution as a function of gas density in our reference calculation, run REF. 




Log n (cm 3 ) 

Figure 3. Ratio of the Ht cooling rate to the total cooling rate, plotted as a function of density, for run REF. 



than about 3% of the total cooling rate. At lower densities, the Hj ions, which are not yet in LTE, contribute less of the 
cooling because they undergo fewer collisions. At higher densities, on the other hand, the effectiveness of Hg" is reduced by 
the significant decrease in its chemical abundance, even though each individual ion contributes more cooling than at 
n = 10 8 cm" 3 . 

If is ineffective, then what about other potential co olants, such as LiH, Hj~ and i ts deu terated counterparts, or other 
ions such as HeH + or H2D + ? As far as LiH is concerned, iMizusawa. Omukai. fc Nishil (|2005l ) have already shown that far 
too little forms for it ever to be a significant coolant, a result which we confirm (c.f., our Fig. (TJ1 with their Fig. lb). The 
contributions made by the other possible coolants are assessed in Figures |4^, and|4p, where for convenience we plot only the 
sum of the contributions from two sets of species. For the ions in one of these sets (Hj, HD + and Dj; Fig. |4^), we have 
cooling functions that should be at least reasonably accurate; for those in the other set (H2D + , HD^ , D^", HeH + , HeD + , He^~, 
LiH + , LiD + and LiH J; Figure HJd), we use the highly approximate treatment described in £14.61 It is clear from the figures 
that none of these species contribute significantly to the total cooling rate, which is unsurprising given their extremely small 
abundances throughout the range of densities examined here (see also Figure [1} • 



5.2 Sensitivity to the choice of Hj collision rate 

As we saw in the previous sub-section, one of the factors preventing H^ from becoming a dominant coolant in our reference 
calculation is the fact that its abundance begins to decrease, owing to the increasing importance of destruction by collisions 
with hydrogen atoms, before its cooling rate has reached its LTE limit. If the Hg~ ion were to reach LTE earlier than we have 
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Figure 4. (a) As Figure[3] but for cooling from H^" and its deuterated counterparts. Note the difference in horizontal and vertical scales 
from Figure |3] (b) As (a), but for the sum of the contributions of the minor coolants discussed in £14.61 
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Figure 5. Ratio of the Hg cooling rate to the total cooling rate, plotted as a function of density, for runs in which the total Hg collision 
rate C was varied. Results are plotted for runs REF (solid line), CI (dashed line) and C2 (dash-dotted line), corresponding to C = C IC { , 
0.1C re f and 10C re f, respectively. 

assumed - in other words, if the low density limit of its cooling rate were to be larger - then cooling would have more 
effect. We have therefore explored the effect of altering C, the total H^" collisional excitation rate coefficient that is the single 
free parameter in our treatment of cooling. The value of C in our reference model, hereafter C rc f , is given by Equation 1 771 
Increasing it or decreasing it compared to this value has the effect of, respectively, increasing or decreasing the low-density 

cooling rate; or, equivalently, decreasing or increasing the critical density at which reaches LTE. 

In Figure O we show the effect that increasing or decreasing C by a factor of ten has on the contribution made by H3" to 
the total cooling rate. As one would expect, the result is rather dramatic. In particular, it is clear that if C = 10 C re f, then 

cooling does contribute significantly to the total cooling rate around densities n ~ 10 8 cm -3 . However, such a large value 
for C seems unrealistic, given that in our expression for C rc f, we are already assuming that collisions occur at the Langevin 
rate. Collisional excitation by electrons could in ideal circumstances give one a large value for C, but it is clear from Figure [T] 
that the electron abundance is orders of magnitude too low in the present case for collisions with electrons to be important. 
Moreover, even if C were as large as 10 C re f , the extra cooling provided by the ions would have only a small effect on the 
temperature evolution, as Figure [6] demonstrates. 

On the other hand, if C is smaller than we have assumed, then Hjj~ cooling has even less effect. Therefore, despite the 
uncertainties in our treatment of cooling at low densities, our main result - that cooling is, in general, unimportant 
- seems robust. 



26 S. C. 0. Glover & D. W. Savin 




rt 5 10 6789 

Log n (cm -3 ) Log n (cm -3 ) 



Figure 6. (a) Temperature evolution as a function of gas density for runs in which the total Hg collision rate C was varied. Results are 
plotted for runs REF (solid line), CI (dashed line) and C2 (dash-dotted line). Note that the solid and dashed lines are not distinguishable 
in this plot, (b) As (a), but focussing on a smaller range of densities and temperatures, to better show the difference between the runs. 
The results of runs REF and CI remain barely distinguishable. 



5.3 Influence of cosmic rays 

In i )5.1l we saw that a major reason for the dramatic decrease in the H^j" abundance at gas densities n > 10 s cm -3 is the loss 
of the remaining H + ions from the gas. In the absence of any H + ions, can no longer be produced directly by reaction 
RA18 (formerly, the dominant production mechanism), while its production from reactions involving H^" or HeH + is also 
disrupted, as the main formation routes for these species also depend upon the availability of H + . Clearly, therefore, one way 
of significantly enhancing the abundance of at high densities would be to provide a source of H + ions at high densities. 
Alternatively, the abundance could also be enhanced if there were a suitable source of ions in dense gas. 

One mechanism capable of producing both H + and H[l~ ions in dense gas is the partial ionization of the gas by an 
external flux of cosmic rays. In view of the many uncertainties and unknowns regarding the c omposition, energy spec trum 
and energy density of the cosm ic rays produced by the earliest supernovae, summarized in Sta cy fc Bromml l|2007l ) and 



Jasche. Ciardi. fc Ensslinl (|2007l ), we use a highly simplified treatment of their effects. We assume that all of the uncertainties 



can be folded into a single free parameter, £h, the cosmic ray ionization rate of atomic hydrogen, and that the cosmic ray 
ionization rates of other atoms and molecules have the same scaling with £h as they are commonly assumed to have in the 
local ISM. Secondary effects resulting from the Prasad- Tarafdar mechanism are treated as outlined in i )3.3l note that we 
assume in this first set of models that Lyman-a photons make a negligible contribution to the secondary photochemical rates. 

In Figure [7Ji, we show how the temperature evolution of the gas changes as we increase £h- We show in the figure results 
from five runs that included cosmic rays: CR1, CR2, CR3, CR4 and CR5, with (h = 1CT 20 , 10" 19 , 10" 18 , KT 17 and 10" 16 s _1 , 
respectively. We also plot the temperature evolution in our reference run REF, for the purposes of comparison. We see that 
as we increase the cosmic ray ionization rate, the gas gets colder. In particular, for £h ^ 1CP 18 s _1 , the gas is able to cool 
to T < 100 K, indicative of the fact that in these runs, HD becomes the dominant low-temperature coolant. This is a simple 
consequence of the ionization produced by the cosmic rays: the additional free electro ns allow more H2 to be produced than 
in our reference run (see Figure [7p), and so the gas can cool to lower temperatures. IStacv fc Bromml ((2007) find a similar 
effect in their recent study of the effects of cosmic rays on primordial star formation, and also show that the effect of the 
cosmic rays on the temperature evolution becomes significant once £h ^ 10 -19 s _1 ; we find a slightly larger critical value here, 
possibly due to the differences in our treatment of H2 cooling. 

In Figure [7jj, we show how the increase in £h affects the contribution that cooling makes to the total cooling rate. 
We see that as the ionization rate increases, H3" cooling becomes steadily less effective at low densities, particularly for 
(h ^ 10~ 18 s -1 , owing to the growing importance of HD cooling at these densities. Above n ~ 10 s cm -3 , however, the 
contribution from H3" cooling increases with increasing £h- As Figure (7JI illustrates, this is a consequence of a significant 
increase in the high density abundance in these runs compared to run REF, which is an expected consequence of the 
greater availability of H + and at high densities in the runs with non-zero (h- 

Figure 0: also demonstrates that if £h 3* 10~ 18 s _1 , then is responsible for > 10% of the total cooling over several 
orders of magnitude in gas density. Moreover, if £h *5 10~ 17 s _1 , then there is a brief period in which it provides > 50% of 
the total cooling. Thus, for cosmic ray ionization rates of this order of magnitude, Hjj~ cooling is clearly significant and Hjj" 
may even be the dominant source of cooling at densities n ~ 10 10 -10 11 cm -3 . 
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Figure 7. (a) Temperature evolution as a function of gas density for runs in which the cosmic ray ionization rate £h was varied. Results 
are plotted for runs REF (upper solid line), CR1 (upper dashed line), CR2 (dash-dotted line), CR3 (dotted line), CR4 (lower solid 
line) and CR5 (lower dashed line), corresponding to £ H = 0.0, 10" 20 , 10~ 19 , 10~ 18 , 10~ 17 and 10" 16 s _1 , respectively, (b) As (a), but 
showing the evolution of the H2 abundance with density in the same set of runs. Note that in this plot the lower solid and dashed lines 
correspond to runs REF and CR1, respectively, while the upper solid and dashed lines correspond to runs CR4 and CR5, respectively, 
(c) As (b), but showing the ratio of the Hg~ cooling rate to the total cooling rate. The lower solid and dashed lines on the right hand 
side of the plot correspond to runs REF and CR1, respectively, the dash-dotted and dotted lines to runs CR2 and CR3, respectively, 
and the upper solid and dashed lines on the right hand side of the plot to runs CR4 and CR5, respectively, (d) As (c), but showing the 
evolution of the Hjj~ abundance with density in the same set of runs. 



How plausible is it that the cosmic ray ionization rate in primordial gas will be as large as 10 17 s x ? This value is 



comparable with the standard estimates for the cosmic ray ionization rate in dense gas in the local ISM (see e.g., Bergin et al 



1999: van der Tak fe van Dishoeckllioool ). and lower than recent estimates of the rate in diffuse gas (see e.g., iMcCall et al 



2003), and so this value is not prima facie unreasonable. However, our requirement that the cosmic rays penetrate to very 
high gas densities means that they must be highly energetic. If we use the same simple model for our collapsing protostellar 



core as in £14.21 then at n 



the core radius is r ~ 10 cm, and the col umn density of the core is N 



To penetrate to this depth, the cosmic rays must have energies of at least 100 MeV (I Stacy &i Bro mm 2007), which means that 
given reasonable assumptions regarding the shape of the cosmic ray energy spectr um, the main co n tribut ion to the cosmic 
ray ionization rate will come from cosmic rays with roughly this energy. Following IStacv fe Bromml (|20070 . we can estimate 
the required energy density in 100 MeV cosmic rays as 



Ucn, 



erg cm 



(100) 



If we fu rther assume that most high-re dshift cosmic rays are pr oduced in the supernova remnants left by pair-instability 
supernovae (jHeger fe Wooslevl I2002T ). then IStacv fe Bromml (J2007J) show that the total cosmic ray energy density Ucn. is 
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related to the cosmological star-forming rate per unit comoving volume, by 

-3 (Pcr\ ( E S n \ /l + z\ 3 / 2 



Ucr(z 



= 2 x 10 



-15 



erg cm 



fp 



2 x 10- 3 M" 1 



V 0.1 J V 10 52 er. 



21 



2 x 10-2 M yr- 1 Mpc" 



(101) 



where pgr is the fraction of the supernova explosion energy, -Esn, that is used to accelerate cosmic rays, and /pisn is the 
number of pair-instability supernovae per solar mass of stars formed. Given reasonable values for pea, ^sn and /pisn, this 
relationship implies that to produce a cosmic ray energy density of order 10 -13 erg cm -3 , we require a star formation rate 
per unit volume of ord er 1 yr" 1 Mpc~ 3 , two to three orders of magnitude larger than current estimates of the population 



III star formation rate ( Yoshida et al. 20031 : Bromm fc Loebl 2006h . Moreover, this estimate assumes that essentially all of the 
cosmic ray energy density is in 100 MeV cosmic rays; if we allow for a significant fraction of cosmic rays with smaller energies, 
as are required to produce the ionization in low density gas in our simplified model, then the required cosmic star formation 
rate increases still further. 

From this argument, we can conclude that any extraga lactic cosmic ray backg round will be too small to produce the 
effect that we desire. How about local sources of cosmic rays? IStacv fe Bromml l|2007l ) show that much higher energy densities 



can be produced close to individual supernova remnants, but to get an energy density of 10 -13 erg cm -3 one would have to 
be within ~ 10 pc of the rem nant, near enough tha t the gas would have been strongly processe d by the ultraviolet radiation 
of the supernova progenitor IjGlover fc Brand! 1200 ll ; I Whalen. Abel, fc Norma n 2004; Susa 2007). Consequently, this scenario 
for producing a high cosmic ray ionization rate also does not appear promising. 

Furthermore, even if we assume that it is possible to maintain a large at high densities, and that Hjj" cooling does 
briefly become dominant, it is possible to show that its effects on the temperature evolution of the gas remain small. In 
Figure [8] we compare the temperature evolution in runs CR5 and CR6. In both runs, we have set £h = 10 -16 s _1 , but in 
run CR6 we have artificially disabled cooling. We see that at densities 10 s < n < 10 12 cm -3 the temperature in run CR5 
is smaller than the temperature in run CR6, as expected. However, the difference is relatively small, and the temperature 
evolution is qualitatively similar in both cases. 

Finally, we have examined the importance of cosmic-ray induced photoionization and photodissocation (the Prasad- 
Tarafdar mechanism, discussed in section [3~3"Tl by peforming several additional simulations. In runs CR7, CR8 and CR9, we 
took £h = 10~ 20 , 10~ 18 and 10 _16 s _1 , respectively, but neglected the effects of the Prasad- Tarafdar mechanism completely. In 
runs CR10, CR11 and CR12, we adopted the same cosmic ray ionization rates, but maximized the effects of the Prasad- Tarafdar 
mechanism by assuming that the Lyman-a photons produced by secondary excitations of atomic hydrogen could propagate 
freely into the core of the protogalaxy, and could contribute to the total secondary photoionization and photodissociation 
rates there (cf. our standard treatment, where we assume that the Lyman-a photons are unable to penetrate into the core). 

In Figure [9^, we compare the temperature evolution of the gas in these six runs with the evolution in runs CR1, CR3, 
and CR5, which have the same values of £h, but which include the effects of cosmic-ray induced photoprocesses. In Figure [9p, 
we show a similar comparison of the ratio of Hj cooling to total cooling in these runs. 
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Figure 9. (a) Sensitivity of the temperature evolution to our treatment of the Prasad- Tarafdar mechanism. Solid lines correspond to 
runs using our default treatment, dashed lines to runs that neglect its effect entirely, and dash-dotted lines to runs that maximize its 
effects by including the effects of Lyman-a photons as if the gas were optically thin to them. The lower, middle and upper sets of curves 
correspond to runs with £jj = 10 — 20 , 10 — 18 and 10 — 16 s — 1 , respectively. Note that many of the lines in this figure are indistinguishable, 
(b) As (a), but showing the ratio of the H^j~ cooling rate to the total cooling rate in the same set of runs. The lower, middle and upper 
solid lines on the right-hand side of this plot correspond to runs CR1, CR3 and CR5, respectively, the lower, middle and upper dashed 
lines on the same side correspond to runs CR7, CR8 and CR9, and while the lower, middle and upper dash-dotted lines on that side 
correspond to runs CR10, CR11 and CR12. (Note that the results of runs CR5 and CR12 are indistinguishable in the plot, as are the 
results of runs CR1, CR7 and CR10). 



Below n = 10 s cm -3 , the Prasad- Tarafdar mechanism has no discernable effect on the evolution of the gas. At higher 
densities, however, its effect is to suppress cooling. More specifically, cosmic-ray induced photodissociation of (reaction 
CP3) reduces the amount of formed via reaction TR3, leading to a reduction in the abundance at these densities, and 
hence an overall reduction in the effectiveness of Hg" cooling. Nevertheless, the effect of the enhanced cooling in runs including 
its effects is slight, as can be seen from Figure [9^,. 

Figure|9]also allows us to assess the impact of the inaccuracy in our treatment of the Lyman-a photons photons produced 
by secondary excitations of atomic hydrogen. Our results indicate that when £h is small, the Prasad- Tarafdar mechanism 
has little effect on the amount of cooling that occurs, and so any inaccuracies in our treatment of it are unimportant. In 
high-^H runs, the Prasad- Tarafdar mechanism is more important, but we find that we obtain very similar results with and 
without the inclusion of the Lyman-a photons, and so any inaccuracy in their treatment is again unimportant. 



5.4 Influence of a radiation background 

In most of our calculations, we have assumed that any external sources of radiation have a negligible effect on the evolution 
of our collapsing protogalactic gas. At the epoch corresponding to the formation of the very first stars, this assumption is 
well-justified: the cosmic microwave background does not significantly affect the gas chemistry at redshifts z < 100 and no 
other sources of radiation yet exist. Once population III star formation has begun, however, the situation changes. Radiation 
from massive population III stars or from their remnants can affe ct pri mordia l gas t hrough a variety of mechanisms, as 
discussed in detail in the recent reviews by ICiardi fc Ferrara HI)!) and ICiardil (|2008r ). In this section, we explore how an 
external radiation field can influence the thermal evolution of primordial gas and affect the role of H J cooling. 



5.4-1 Ultraviolet radiation 



One of the most important forms of radiative feedback in the high-redshift Universe is the build-up of a soft ultraviolet 
background at photon energies hv < 13.6 eV. Photons from this background that are absorbed in the Lyman and Werner 
band transitions of H2 can cause photodissociation, and since these photons can propagate to large cosmological distances 
through the intergalactic medium, the strength of the background and the size of the associated photodissociation rate can 
both become considerable. This soft UV background is therefore expected to have a significan t effect on th e evolution of 
prim ordi al gas within protogalax ies IjHaiman. Rees. fc Loebl 1 19971 : lHaiman. Abel, fc Rccs 2000J; but see also I Wise fc Abel 
2007 and O'Shea fc Norman [2008 1 for evidence that the effects of the Lyman- Werner background may be less important than 
previously supposed). 
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Figure 10. (a) Temperature evolution as a function of gas density in runs REF (solid line), UV1 (dashed line), UV2 (dash-dotted line) 
and UV3 (dotted line). The strength of the ultraviolet background in these runs was J21 = 0.0, 10 — 4 , 10 — 2 and 1.0, respectively, and 
the gas was assumed to remain optically thin throughout its evolution, (b) As (a), but for runs REF, UV4, UV5 and UV6. Runs UV4, 
UV5 and UV6 had the same UV background field strengths as runs UV1, UV2 and UV3, respectively, but in this case we assumed that 
the gas was optically thick in the Lyman- Werner lines of H2 and HD. The results of the four runs are indistinguishable in the plot. 



To investigate the impact of such an ultraviolet background on our results, we have run several models with non-zero 
backgrounds: runs UV1, UV2, UV3, UV4, UV5 and UV6. As previously noted in £13.21 we assume that the spectral shape of 
the background is that of a diluted 10 K black-body, with a sharp cutoff at 13.6 eV. The only free parameter is then the 
normalization of this spectrum. We choose to normalize it by specifying its strength at the Lyman limit. In runs UV1, UV2 
and UV3, the field strength is J21 = 10 -4 , 10 -2 , and 1.0, respectively, where J21 is the flux at the Lyman limit in units of 
10 -21 erg s _1 cm -2 Hz -1 sr _1 . In these three runs, we assume that self-shielding by H2 and HD is not effective, and that the 
gas remains optically thin to the external radiation field throughout the simulation. In runs UV4, UV5 and UV6, the field 
strength is the same as in runs UV1, UV2 and UV3, respectively, but we assume that H2 and HD self-shielding is so effective 
that the H2 and HD photodissociation rates are negligible. The true behaviour of the gas lies between these two limiting cases. 

In Figure [T0k . we show how the gas temperature evolves in optically thin runs UV1, UV2 and UV3, as well as in run 
REF for comparison. The corresponding behaviour in runs UV4, UV5 and UV6 is shown in Figure [TOb . In the optically thin 
runs, the effect of the ultraviolet background is to increase the temperature of the gas at early times; quite dramatically so in 
the case of run UV3, where the minimum temperature reached by the gas is T ~ 900 K, compared to only T ~ 200 K in our 
reference calculation. This temperature increase is an obvious consequence of the photodissociation of H2 in low density gas, 
as can be seen clearly by comparing these results with those from the runs in which H2 photodissociation was assumed to be 
negligible. 

As far as H^" cooling is concerned, Figure [TT] demonstrates that it remains ineffective in both sets of runs. In the 
optically thin runs, H2 dissociation at early times reduces the effectiveness of H2 cooling, but also significantly reduces the Hjj" 
abundance. The net effect is to reduce the amount of cooling coming from H J to below the level that it has in our reference 
run. At densities n > 10 9 cm -3 , on the other hand, the effect of the UV background is to enhance cooling from H3 . This 
occurs because the H + abundance does not decline so quickly in the runs with the UV background, owing to the higher gas 
temperature, and so more free protons are available for making Hg at very high densities. This effect boosts the contribution 
of H^J" to the cooling rate in run UV3 by about a factor of four compared to our reference run. Despite this, however, the 
contribution from H^" remains unimportant. 

In runs UV4, UV5 and UV6, the behaviour of the contribution from H3 is somewhat different. The H2 abundance in these 
runs evolves in almost the same manner as in run REF, as does the temperature. Therefore at early times, there is no difference 
in the Hj contribution. At n > 10 7 cm -3 , however, a difference does become apparent between run REF and runs UV5 and 
UV6, with the H3" contribution falling off faster the more the strength of the UV background is increased. This behaviour 
is again a result of a change in the behaviour of the H + abundance at high densities. In this case, the H + abundance falls 
off faster at high density when the UV field strength is increased. This occurs because the ultraviolet background maintains 
a higher Li + fraction in the gas than in our reference calculation. Because the Li + abundance is larger, it contributes more 
free electrons to the gas, and so the fraction of Hjj" ions that are destroyed by dissociative recombination in the high density 
regime becomes larger. Consequently, fewer of the H + ions destroyed by reaction RA18 are recycled by reaction TR17, and so 
the H + and H^ abundances fall off more rapidly than in our reference model. A similar effect is not seen in the optically thin 
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Figure 11. (a) Ratio of the H^" cooling rate to the total cooling rate in runs REF (solid line), UV1 (dashed line), UV2 (dash-dotted 
line) and UV3 (dotted line), (b) As (a), but for runs REF (solid line), UV4 (dashed line), UV5 (dash-dotted line) and UV6 (dotted line). 



runs because it is more than offset by the effects of the higher gas temperature, which increases the rate of reaction TR17, 
thereby decreasing the fraction of H^" ions that are destroyed by dissociative recombination. 



5.4.2 X-rays 



X-rays are another form of radiation that can affec t the evolution of primordial gas. By providing an additiona l source of ioniza- 



tion in dense gas, they can promote H 2 formation IjHaiman. Rees. fc Loeblll997l ; lHaiman. Abel, fc Ree s 2000; G lover fc Brand 



20031 ; iMachacek. Bryan, fc Abellbooj ). in much the same manner as the cosmic rays considered in section [531 However, just as 
in the case of cosmic rays, for X-rays to materially affect the importance of cooling, they must be able to penetrate deeply 
into the collapsing protostellar core, to densities n ~ 10 10 cm -3 or more, corresponding to column densities N ~ 10 25 cm -2 
or more. As a core with this column density is opaque even to hard X-ray photons, the photon flux required to produce a 
significant photoionization rate is considerable. 

Given reasonable assumptions regarding the shape of any high-redshift hard X-ray backgrou nd, the dominant contribution 



to the photoionization rate at n = 10 10 cm 3 comes from X-ray photons with energies E ~ 3 keV ( Yan. Sadeghpour. fc Dalgarnol 



1998). The hydrogen ionization cross-section at this energy is approximately 10 -24 cm 2 , and so for a column density N = 
10 25 cm -2 , the gas has an optical depth r ~ 10 for a 3 keV photon. Each of the photons t hat is absorbed is responsibl e 
for roughly 100 ionizations, once the effects of secondary ionization are taken into account ( Dalgarno. Yan. fc Liu 19991 ). 
Therefore, a flux Fx of 3 keV photons incident on the cloud exterior produces an ionization rate Rx given approximately by 



Rx 



100 (10" 
1.1 x 10" 



l F x Au) ■ 
'Fxs- 1 , 



where the second line follows if we assume that Av ~ 1 keV /h. Inverting this expression, we obtain 
Rx 



9 x 10" 



10- 



photons s 1 cm 2 Hz 



This corresponds to an X-ray background field strength at E — 3 keV of 
Rx 



4.2 x 10" 



U0- 17 a- 1 ) 



-1 -2 TT -1 

erg s cm Hz , 



(102) 



(103) 



(104) 



which is orders of m agnitude larger than any plausible range of values for the high-redshift X-ray background (see e.g., 
Glover fc Brandll2003r ). We can therefore rule out an extragalactic X-ray background as the source of the required hard X-ray 
photons. 

Furthermore, although hi gher X-ray fluxes can be maintained close to strong X-ray sources such as miniquasars ( Haiman. Abel, fc Rees! 
2000l : iKuhlen fc Madaul 12005). even n this case it is diff icult to produce a s ignific ant ionization rate. For example, the lumi- 
nosity at 3 keV of the model miniquasars considered by Kuhlen fc Madau ( 20051) , assuming their hardest spectral model, is 
Lx ^ 10 22 ergs -1 Hz" 1 . To see a flux Ix = 4.2 x 10 17 erg s 1 cm 2 Hz 1 from this miniquasar, one must therefore be within 
a distance 



32 S. C. 



0. Glover & D. W. Savin 



o 
o 




5 10 rt 4 6 8 10 12 

Log n (cm -3 ) Log n (cm -3 ) 

Figure 12. (a) Temperature evolution of the gas for various different initial densities. Results are plotted for runs REF, Nl and N2, 
which have initial densities n; = lcm -3 (solid line), rij = 0.03cm -3 (dashed line) and n; = 30cm -3 (dash-dotted line), respectively, (b) 
As (a), but showing the ratio of the cooling rate to the total cooling rate for the three different models. 




of it. Gas this close to the miniquasar would have been strongly processed by the ultraviolet radiation of its progenitor, and 
is not a promising place to expect to find further star formation. 

We therefore consider it likely that the hard X-ray flux seen by most collapsing protostellar cores will be far too small to 
significantly affect the production of at high densities. Moreover, even if somehow a sufficiently large flux was produced, 
we would expect its effects to be very similar to those of the cosmic rays considered in H5.3I Accordingly, we do not consider 
it necessary or time-efficient to examine the effects of a hard X-ray flux in any greater detail. 



5.5 Sensitivity to initial conditions 

5.5.1 Altering the initial density and temperature 

In order to verify that our main results are not sensitive to the initial temperature or density assumed in our models, we 
have performed several calculations with different initial densities or temperatures. In runs Nl and N2, we set rij = 0.03 and 
30 cm -3 , respectively, while keeping all of the other input parameters fixed. The effect that this has on the thermal state of 
the gas is illustrated in Figure 112b . where we compare the temperature evolution in runs Nl and N2 with the evolution in 
our reference calculation, run REF. It is clear fro m the Figure that the temperature evolution o f the three runs is strongly 
convergent, in line with previous findings (see e.g., Palla. Salpeter. fc Stahler 1983 ; Omukai 200(J). Consequently, it comes as 
no surprise to find that the contribution that H;J~ makes to the cooling in the three runs is not greatly affected by the choice 
of ni, as shown in Figure [12b . 

In runs Tl and T2, we set 2] = 100 and 10000 K, respectively, and performed a similar comparison, which is illustrated 
in Figure [TS] Again we find that the results of runs Tl, T2 and our reference run REF converge well, although the differences 
in this case are slightly larger than those that occur when ni is varied. 



5.5.2 Altering the initial fractional ionization 

We have also explored the effect of altering the initial fractional ionization of the gas. In runs XI, X2 and X3, we set the 
initial H + abundance to 10 -6 , 10 -2 or 1.0, respectively. We also rescaled the initial D + abundances by a similar amount. 
However, the initial He + abundance was not altered. As Figure [T4k demonstrates, altering the initial fractional ionization in 
this way has a dramatic effect on the temperature evolution of the gas. In run XI, the low abundances of free electrons and of 
H + delay the formation of H2 and limit the amount that can form. The gas therefore undergoes a period of adiabatic heating 
that lasts for much longer than in our reference calculation. Furthermore, once enough H2 has formed to cool the gas, the gas 
temperature remains significantly higher than in run REF. However, the two runs eventually converge at n ~ 10 10 cm -3 , as at 
this density three-body processes dominate the formation of H2, and so the H2 abundance, and hence the thermal evolution 
of the gas, are no longer sensitive to the fractional ionization. 

In runs X2 and X3, on the other hand, the enhanced initial ionization has the effect of promoting the formation of H2, 
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Figure 13. (a) Temperature evolution of the gas for various different initial temperatures. The results plotted are from runs REF, Tl 
and T2, which had initial temperatures of T- = 1000 K (solid line), T ; = 100 K (dashed line) and T ; = 10000 K (dash-dotted line), 
respectively, (b) As (a), but showing the ratio of the H^" cooling rate to the total cooling rate for the three different models. 



and cooling the gas more than in our reference run. Moreover, the cooling provided by this extra H2 is sufficient to lower the 
gas temperature to a level at which chemical fractionation between HD and H2 becomes highly effective. The resulting boost 
in the HD abundance allows it to dominate the gas cooling and to cool the gas to a lower temp erature than could be reached 
by H2 cooling alone. Note, however, that in contrast to the results of pre vious studies (see e.g., Nakamura fc Umemuralfeood ; 
Nagakura fe OmukailEx)! ; Ijohnson fc Brommlliooj ; lYoshida et al.ll2007l ). the gas does not reach the CMB temperature This 



is a consequence of the relatively rapid rate of collapse assumed here (c.f. ^5.7p and of the re vised treatment of H2 cooling 
used in this work, which tends to render H2 cooling less effective, as explored in more detail in Glover fc Abel 1 20081 ) . In any 
case, the period of HD dominance lasts for only a short time, as illustrated in Figure [14b . As the HD level populations near 
their LTE values, HD cooling becomes less effective and the gas starts to warm. As it warms, chemical fractionation becomes 
less effective and the HD:H2 ratio declines. Once the gas has warmed to T ~ 200 K, which occurs at a density between 10 5 
and 10 6 cm -3 , the amount of HD remaining in the gas is no longer sufficient to maintain HD as the dominant coolant; H2 
becomes dominant once more, with HD thereafter relegated to a minor role. 

In Figure IT4"b . we show how the ratio of the H3" cooling rate to the total cooling rate varies in these runs. We see that if 
the initial fractional ionization is lowered, the contribution of Hj to the cooling is lowered at n < 10 9 cm~ 3 and increased at 
n > 10 9 cm -3 relative to our reference calculation. The lowered importance of H^ cooling at low densities is a result of the 
higher gas temperature: the greater temperature sensitivity of the H2 cooling rate compared to the H3" cooling rate makes 
the former more effective in comparison to the latter as the temperature is raised. At high densities, the gas temperatures 
converge, and the difference between the runs has a different cause: the lower H2 abundance (discussed at the beginning of 
this subsection) increases the time required to convert all of the H + to H^. For this reason, the H + abundance in n ~ 10 9 - 
10 10 cm -3 gas in the low ionization run is higher than in the reference run, with the result that the H^" abundance and the 
Hjj~ contribution to the cooling are also marginally higher. 

If the initial fractional ionization is increased, an interesting effect occurs. The contribution of H^~ to the total cooling 
rate is considerably suppressed at densities n < 10 6 cm -3 and n > 10 8 ' 5 cm -3 compared to the contribution in our reference 



run, but is slightly enhanced at densities 10 6 < n < 10 8 ' 5 cm 3 . The reduction in the effectiveness of Hj at low densities in 

n6 ™-3 



these runs is a result of the previously noted strong enhancement of the HD cooling rate at n < 10 cm , as can clearly be 
seen by comparing Figures [Tib and I14fc . On the other hand, the suppression of H3" cooling at n > 10 8 ' 5 cm -3 results from 
the rapid loss of H + in the gas driven by reaction RA18, which occurs more rapidly than in the reference run owing to the 
greater H2 abundance in these runs. Between these two density regimes, there is a small range of densities in which the H J :H2 
ratio remains relatively large, and where the gas temperature is ~ 300-500 K. At these temperatures, the ratio of the Hjj" to 
H2 cooling rates is larger than at T = 1000 K, owing to the greater temperature dependence of the H2 cooling rate, but the 
temperature is not low enough for the HD abundance to be significantly enhanced by chemical fractionation. These conditions 
are therefore close to ideal for Hj cooling, and the fact that even in this case the Hjj" contributes no more than a few percent 
of the total cooling helps to strengthen our conclusion that it is generally of little or no importance. 




T3 




Log n (cm 3 ) 

Figure 14. (a) Temperature evolution of the gas in runs REF (solid line), XI (dashed line), X2 (dash-dotted line), and X3 (dotted line), 
which had initial H+ abundances of £ H + = 2.2 X 10 -4 , 10 — 6 , 10 — 2 and 1.0, respectively, (b) As (a), but showing the ratio of the HD 
cooling rate to the total cooling rate in these four runs, (c) As (b), but for the ratio of the H^" cooling rate to the total cooling rate. 



5.5.3 Changing the elemental composition 

Finally, we have examined the effect of changing the elemental composition of the gas by removing all of the deuterium (run 
ELI), lithium (run EL2) or both (run EL3). Although not physically realistic, these runs do provide a convenient way to 
examine the roles that the deuterium and lithium play in the overall thermal evolution of the gas. In Figure [T5l we show how 
the gas temperature evolves in run ELI in comparison to our reference calculation, run REF. We see that the omission of 
deuterium has a noticeable effect on the temperature evolution at densities 10 3 < n < 10 6 cm -3 , and a very slight effect on the 
temperature at densities n > 10 12 cm -3 . Examination of the contribution of HD cooling to the total cooling rate in run REF 
(plotted as the solid line in Figure [14b ) demonstrates that at these densities, HD contributes significantly to the total cooling 
rate; indeed, at its peak at n ~ 10 4 cm -3 , it contributes almost a third of the total cooling. At first sight, this result is rather 
surprising, as it is often assumed that HD cooling is unimportant in primordial gas unless the gas has a large initial fractional 
ionization, as in runs X2 or X3 discussed above. However, it appears that the conventional wisdom is wrong on this point ; 



our results here are consistent wit h those of previous studies that have included HD (see e.g.. lBromm. CoppL fc Larson 



point; 

booi 



Mizusawa. Omukai. fc Nishi 2005), and show that although HD is never the dominant coolant, it does contribute enough to 



the total cooling rate at densities n ~ 10 4 -10 5 cm -3 to warrant inclusion in future models of population III star formation. 

In run EL2, the temperature evolution is essentially the same as in run REF, while in run EL3, the evolution is the same 
as in run ELI, indicating that lithium does not play a significant role in the cooling of the gas, in agreement with the results 
presented in i)5.1l note that runs EL2 and EL3 are not plotted in Figure 1151 as they would not be distinguishable from the 
existing lines. In Figure [TBI we investigate the size of the contribution that H^~ makes to the total cooling rate in runs ELI, 
EL2 and EL3; we also plot the result from run REF for purposes of comparison. The most obvious point to note is that in 
the runs without lithium, namely EL2 and EL3, the Hjj" contribution no longer falls off sharply at high densities, although it 
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Figure 15. Temperature evolution as a function of density in runs REF (solid line), which had the standard cosmological deuterium 
abundance, and ELI (dashed line), in which the deuterium abundance was set to zero. 
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Figure 16. Contribution of Ht cooling to the total cooling rate in runs with no deuterium (ELI; dashed), no lithium (EL2; dash-dotted) 
and no deuterium or lithium (EL3; dotted), along with the results of our reference run REF (solid) for comparison. Note that the dashed 
and solid lines are barely distinguishable from each other; similarly, the dotted and dash-dotted lines are not distinguishable in this plot. 



remains too small to be significant. This is easy to understand, given our previous discussion of the chemistry of the gas at 
high densities and low fractional ionizations (see SJ57TJ. As previously noted, in the absence of lithium, the net rate of removal 
of H + ions from the gas via reaction RA18 decreases as the fractional ionization decreases, since an increasing fraction of the 
H3" created by reaction RA18 is converted back to H2 and H + by reactions TR17 and CT3, instead of being destroyed by 
reactions DR4 and DR5. Therefore, the H + removal timescale, ti oss , remains longer than the free-fall timescale throughout 
the simulation, and the H + abundance falls off gradually at high densities; the rapid fall-off that occurs in our reference run 
once iioss < iff does not take place. Consequently, the corresponding rapid fall-off in the Hjj~ abundance also does not occur, 
as the Hjj~ formation rate never becomes negligible. 

Although the presence or absence of lithium does not affect the conclusions of our current study, as in either case cooling 
from Hjj~ is negligible, it is clear that if one is interested in determining the fracti onal ionization of th e gas accurately at very 
high densities, it is vital to include Li and Li + in the chemical model (see also iMaki fc Susall2004l . who come to a similar 
conclusion) . 
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Figure 17. Ratio of the Hg cooling rate to the total cooling rate in runs REF ( solid line) and RA (dash ed li ne). These runs use d 
values for Airais that differed by four orders of magnitude, and that were taken from lGerlich fc Horning! |l992) and lStancil et al 1 l| 19981 ). 
respectively. 



5.6 Sensitivity to uncertainties in the chemical rate coefficients 

5.6.1 Reaction RA 18 

As we have already discussed in Sj3] large uncertainties exist in the rate coefficients of a number of the processes included in 
our chemical model. Of particular relevance to this paper is the huge uncertainty that appears to exist in the value of the rate 
coefficient for reaction RA18, the formation of Hjj" by the radiative association of H2 with H + . In Figure | 17| we compar e the 
contribution of Hjj" cooling in our reference calculation (solid line), in which we adopt the large Gerlich fc Horningi (,1992) rate 



coefficient for react ion RA18, with the contribution of H^ cooling in run RA, a similar calculation that adopts the smaller 
Standi et al. (1998) rate coefficient (dashed line). We see that Hjj" is less effective in the latter case, and that its effectiveness 



also peaks at a later point in the simulation. 

However, the reduction in the contribution is less than one might expect given the very large difference in ^rais 
between the two runs. The reason for this is that although the rate of H3" formation by radiative association is strongly 
suppressed, other formation mechanisms remain unaffected. The HjJ" formed in run RA is produced primarily by the 
familiar reaction 

H+ + H 2 ^H+ + H, (106) 
with the necessary Hj coming mainly from reaction RA3, namely 

H + (107) 

The persistence of a significant H^{" contribution at later times in run RA than in run REF is a clear consequence of the fact 
that the rate at which H + ions are removed from the gas by reaction RA18 is smaller in the former run than in the latter. 
As a result, the familiar rapid fall-off in the H + abundance that occurs once the H + removal timescale, tioss, becomes smaller 
than the free-fall time (see £15. ip takes place at a later time in run RA than in run REF, and hence the corresponding fall-off 
in the Hjj~ abundance also occurs at a later point in the evolution of the gas. 



5.6.2 Reactions TBI and TB2 



Large uncertainties also exist in the rate coefficients for three-body H2 formation (reactions TBI and TB2). As illustrated in 
Figure [181 these uncertainties significantly affect the temperature evolution of the gas at densities n > 10 8 cm -3 , particularly 
the uncertainty in the rate of reaction TBI. The use of larger values for the three-body reaction rates leads to faster production 
of H2 at high densities, and hence a greater H2 cooling rate. This has the effect of slowing the rise in the gas temp erature at 
these densities, which may affect the ability of the gas to fragment at late times (see lClark. Glover, fc Kless en 2008). However, 
in the present context, the effect of the faster H 2 formation rates is to make cooling by H^" even less effective at late times 
than in our reference calculation, as demonstrated in Figure 1191 where we examine the effect that varying the rate of both 
reactions has on the contribution that Hjj" makes to the total cooling rate. 
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Figure 18. (a) Temperature evolution as a function of density, plotted for runs REF (solid line), 3B1 (dashed line) and 3B2 (dash- 
dotted line). We use dif f erent values for the rate of reaction TBI in these three runs. In run RE F, we use the rate coefficient from 
lAbel Bryan, fc Normanl (120021) . in run 3B1 the rate coefficient from lPalla. Salpeter, fc Stahlerl l ll983l ) and in run 3B2 the rate coefficient 
from lFlower fc Harris! fl2007l) . (b) As (a), but for runs REF (solid line), 3B3 (dashed line) and 3B4 (dash-dott ed line). In these runs, the 
rate of reaction TB2 is varie d. In run REF, we use the rate coefficient from l Palla, Salpeter, fc Stahlerl lll983l ). while runs 3B3 and 3B4 
use the rate coefficients from lColi cn & Wcstbcrg (1983]) and lFlower fc Harri si l|2007h . respectively. 




5.6.3 Reaction AD1 



The uncertainty in the rate of reaction AD1 discussed in Glover. Savin, fc Jappsen 1 20061 ) also affects the temperature evolution 
of the gas. We have examined two cases, runs AR1 and AR2. In run AR1, we use a rate coefficient /cadi = 0.65 x 10 -9 cm 3 s _1 
for reaction AD1, taken from I Glover. Savin, fc Jappsenl |2006l ), which is a plausible lower limit on the rate coefficient. In 
run AR2, on the other hand, w e use a rate coefficient &adi = 5.0 x 10 -9 cm 3 s _1 for reaction AD1, again taken from 
Glover. Savin, fc Jappsenl l|2006l ). which is a plausible upper limit. 

In Figure [20k . we show how the temperature of the gas evolves in these two runs, as well as in run REF for comparison. 
It is clear from the figure that the uncertainty in /cadi has only a slight impact on the temperature evolution of the gas. In 
Figure 120b we show a similar plot of the ratio of the cooling rate to the total cooling rate. Again, the rate coefficient 
uncertainty has only a small effect. This res ult is in line with previous work s howing that these are unimportant when starting 
from cold, low ionization initial conditions I Glover. Savin, fc Jappsen 20061 ) . If, instead, we star t with hot, ionized gas, then 
the effect of the uncertainties on the temperature evolution is much greater I Glover fc Abel! 2008). However, even in this case, 



the largest effects are seen at densities n< 10 4 cm 3 , far below the densities at which could conceivably become important, 
and so our basic conclusion regarding the unimportance of cooling remains unaffected. 
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Figure 20. (a) Temperature evolution as a function of density in runs REF (solid line), ARl (dashed line) and AR2 (dash-dotted line). 
In run ARl, the value for the rate coefficient of reaction AD1 is chosen so as to minimize H2 production, while in run AR2, H2 production 
is maximized, as discussed in more detail in the text, (b) As (a), but showing how the ratio of the Hg~ cooling rate to the total cooling 
rate varies in these three models. 



5.6.4 Reactions RA20 and CD 26 

Finally, we have investigated the effects of varying the rates of two of the key reactions involved in the LiH J chemistry: LiH^ 
formation by radiative association of Li + and H2 (reaction RA20) and the collisional dissociation of LiHj by H2 (reaction 
CD26). As we have already discussed in section [3.1.111 the rate coefficients for both of these reactions are unknown. In 
our reference model, we adopted values of &;ra20 = 1CP 22 cm 3 s" 1 and &CD26 = 1-0 x 10~ 9 exp (—2^2 J cm 3 s -1 for the 
rate coefficients; i.e., a small value for reaction RA20 and a large value for reaction CD26. These choices serve to minimize 
the role played by LiH^" in the chemical evolution of the gas, and hence we considered them to be the most conservative 
options in the circumstances. In run LP1, we adopted instead a much larger value for the rate coefficient of reaction RA20, 
feRA20 = 10 -17 cm 3 s _1 , but kept the same value for fccD26 as in our reference model. In run LP2, we used our reference value 
for A:ra2o, but adopted a much smaller value for the rate coefficient of reaction CD26: fccD26 = 1-0 X 10~ 13 exp f— cm 3 s _1 . 
Finally, in run LP3, we altered both rate coefficients, using the larger value for A;ra20 and the smaller for &CD26- 

In Figure [21] we show how the contribution of Hjj" cooling to the total cooling rate varies with density in runs LP1, LP2 
and LP3, along with run REF for comparison. It is clear that in runs LP1 and LP2, the behaviour is essentially the same as 
in our reference run. This can be understood if we consider the timescale on which Li + ions are destroyed by the reaction 
sequence 

Li+ + H 2 -> LiUt + 7, (108) 
followed by 

LiUt + e" -> products. (109) 
This sequence of reactions removes Li + on a timescale 

tloBB = " T- (HO) 

KRA20WH2 JDR 

where here /dr is given by 

_ (fcpR19 + fc D R2Q + fc DR .2l)w e - (111) 

(&DR19 + &DR20 + fcDR2l)n c - + fccD26^H 2 ' 

and represents the fraction of LiHj ions destroyed by dissociative recombination, rather than by collisional dissociation. In 
our reference run, at a gas density n = 10 10 cm -3 , the temperature T ~ 1000 K and the electron abundance x e - ~ 5 x 1Q -11 , 
and hence /dr — 1.4 x 10~ 7 and ti oss ~ 1.4 x 10 19 s, many orders of magnitude longer than the dynamical timescale. Thus, in 
our reference run, the LiH^" chemistry has almost no effect on the Li + abundance. In run LP1, A:ra20 is a factor of 10 5 larger 
than in our reference run, and in run LP2, fccD26 is a factor of 10 4 smaller, and so in both runs, ti OBB is significantly reduced. 
However, it still remains far greater than the free-fall timescale, which is ~ 10 s at this density. Thus, in these runs, the 
LiHj chemistry still has almost no effect. 

In run LP3, however, where we both increase &ra20 and decrease A?cd26, ^loss is reduced by a factor of 10 9 , making it 
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Figure 21. (a) Evolution of the ratio of the Hjj~ cooling rate to the total cooling rate as a function of density in runs REF (solid line), 
LP1 (dashed line), LP2 (dash-dotted line) and LP3 (dotted line); note that the first three of these lines are not distinguishable in the 
plot. In run LP1, the rate coefficient for reaction RA20 was increased by a large factor compared to our reference value, while in run 
LP2, the rate coefficient for reaction CD26 was decreased by a large factor. In run LP3, both changes were made. 



iioss ~ 1.4 x 10 10 s at n = 10 cm -3 , of the same order of magnitude as the free-fall collapse time. Moreover, as t\ OBB scales 
with density as ti oss oc while the free-fall time scales as ts oc n -1 / 2 , ti oss becomes smaller than iff at densities not very 
much greater than 10 10 cm -3 . Therefore, in this run, the LiH^ chemistry does have a noticeable effect on the Li + abundance, 
reducing it by a factor of roughly fifty in comparison to the reference run by the end of the simulation. This reduction in the 
Li + abundance reduces the number of free electrons available for destroying H^, and so limits the rate at which its abundance 
declines at very high densities, much as in runs performed without any lithium (c.f. section f5.5.3|l . Nevertheless, it is clear 
from Figure [21] that this change in the lithium chemistry does not change our basic results: H J cooling remains ineffective, 
albeit somewhat less ineffective at high densities than in our reference run. 



5.7 Sensitivity to the details of the dynamical model 

A major limitation of our current study is the highly simplified dynamical treatment that we use in our one-zone model. In 
all of the calculations that we have presented so far, we have assumed that the gas is collapsing grav itationally at the free-fall 



rate. However, it is well k nown from more detailed three-dimensional hydrodynamical models (e.g.. lAbel. Bryan, fc Norman 



20021 ; lYoshida et al.|[200f3) that in realistic primordial clouds, the collapse speed is significantly slower than the free-fall rate 



owing to the non-negligible gas pressure. As a result, the gas takes longer to evolve than assumed here, and the impact 
of compressional heating is also somewhat smaller. A crude way of taking this into account in a one-zone calculation is to 
artifically slow down the collapse of the gas. In other words, instead of assuming that the density evolves according to the 
standard free-fall relationship 

^ = £- (112) 

dt tff 

where ts is the free-fall time, we can instead assume that 

with r\ < 1. The effect of this change is to lengthen the time taken for the gas to collapse to any given density by a factor 

(iA?)- 

In Figure 122b we show the effect that slowing the collapse in this fashion has on the temperature evolution of the gas 
by comparing the results of three runs, DYN1, DYN2 and DYN3, with 77 = 0.6, r\ = 0.3 and r\ = 0.1, respectively, with the 
results of our reference run REF. We see that reducing r\ leads to a reduction in the temperature of the gas throughout the 
run, a simple consequence of the reduction in the compressional heating rate. Interestingly, in the r\ = 0.3 and r\ = 0.1 models, 
the reduced heating allows the gas to cool to temperatures low enough for chemical fractionation to strongly enhance the 
HD fraction, allowing HD cooling to further cool the gas down to temperatures close to Tomb- The fact that this effect is 



not seen in more realistic hydrodynamical models (e.g.. iBromm. Coppi. fe Larson 2002) suggests that in these models we are 



overestimating the extent to which gas pressure slows the collapse, at least at early times. 
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Figure 22. (a) Temperature evolution as a function of density in three runs in which the collapse parameter r\ was varied. Results are 
plotted for runs REF (solid line), DYN1 (dashed line), DYN2 (dash-dotted line) and DYN3 (dotted line), which had r? = 1.0, 0.6, 0.3 
and 0.1, respectively, (b) As (a), but showing the ratio of the cooling rate to the total cooling rate in the same three runs, (c) As 
(a), but showing the evolution of the Hg~ abundance in the three runs. 



In Figure 122b . we show how reducing 77 affects the contribution that H f cooling makes to the total cooling rate. In the 
77 = 0.6 run, H^~ cooling is slightly more effective than in the reference run at densities 10 5 <n<10 8 cm -3 . This is a consequence 
of the slightly lower temperature of the gas at these densities in run DYN1 compared to run REF, which decreases the H2 
cooling rate more than the cooling rate. At lower densities, cooling in run DYN1 is slightly less effective than in 
run REF, as the lower gas temperature makes HD cooling more effective at these densities in the former run than in the 
latter. At higher densities, Hg~ cooling becomes far less effective in run DYN1 than in run REF owing to a more rapid falloff 
in the abundance, as illustrated in Figure 122b . As in earlier runs, the reason for this rapid falloff is that the timescale 
for the removal of H + from the gas by conversion to H^" followed by destruction of the H^~ by dissociative recombination, 
becomes shorter than the dynamical timescale of the gas at these densities. In most of the runs that we have studied, the high 
gas temperature at these densities allows reaction TR17 to interfere with this process by converting most of the to Hj, 
following which reaction CT3 restores the original proton to the gas. As we saw in i]5.1l the effect of this is to lengthen the 
time required to remove all of the H + , which delays the precipitous falloff until n ~ 10 11 cm -3 . However, in run DYN1, the 
lower gas temperature means that reaction TR17 is less effective, and so the delay is much shorter. Moreover, the dynamical 
timescale itself is longer. Consequently, the rapid falloff occurs at a lower density. 

A similar kind of behaviour is seen in runs DYN2 and DYN3. However, in these runs the much lower temperature at low 
densities renders HD cooling dominant for significantly longer, while the lower temperature of the gas at high densities, plus 
the longer dynamical timescale, allow the rapid falloff in the H^" abundance to occur sooner. Consequently, H^" cooling never 
becomes significant in these runs. 
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6 CONCLUSIONS 

We have examined the contribution that cooling makes to the total cooling rate of gravitationally collapsing primordial 
gas in a wide range of different models using a newly-developed cooling function along with the most detailed model of 
primordial gas chemistry published to date. Our results demonstrate that in general H^" cooling is not important, although it 
comes close to being so at densities n = 10 7 -10 9 cm -3 , contributing at its peak a few percent of the total amount of cooling. 
We come to this conclusion despite making several assumptions (regarding the collapse rate of the gas, the formation rate of 
Hj{" by radiative association, and the collisional excitation rate of its excited vibrational states) that favour HjJ" cooling, and 
thus we have confidence that our conclusion is robust. 

As H^" comes so close to being important, it is instructive to examine why it ultimately fails to dominate. This can be 
ascribed to a combination of two main effects. First, at high densities (n > 10 s cm -3 ), the gas temperature becomes high 
enough to make the endothermic reaction TR17 

+ H + H 2 (114) 

the most important destruction mechanism, which significantly suppresses the Hjj~ abundance at these densities. Second, 
the formation rate of is strongly suppressed by the rapid removal of H + from the gas at densities n > 10 11 cm -3 . At these 
densities, the fractional ionization of the gas is so low that the main loss route for the H + is conversion to H^, followed by Hjj~ 
dissociative recombination, and once the timescale for H + removal via this combination of reactions becomes short compared 
to the free-fall time, the H + abundance decreases by orders of magnitude within a short space of time, effectively switching off 
the formation of H J . Moreover, destruction of H3" by dissociative recombination remains effective despite the fall-off in x H + 
thanks to the contribution of electrons from ionized lithium, Li + , which for n>3x 10 8 cm -3 is the most abundant positive 
ion in the gas. 

In our study, the only situation in which we found Hjj" cooling to be important is if the gas is illuminated by a strong 
flux of cosmic rays or X-rays. If the incident flux is strong enough to produce an ionization rate >10 -18 s _1 at densities 
n ^ 10 10 cm -3 , then the high density H^j" abundance can be significantly increased, and H^~ can even become the dominant 
coolant. However, the necessary flux of cosmic rays or hard X-rays is difficult to produce in the high-redshift Universe. As 
the estimates in §§ 15.31 and 15.4.21 demonstrate, the required flux is orders of magnitude greater than the size of any plausible 
extragalactic background, and will only be achieved within gas that is very close to a local source (i.e., within 5-10 pc). 
However, gas that is this close to a supernova remnant or miniquasar will have been strongly affected by radiative feedback 
from the progenitor star and so is not a promising place to expect to find ongoing population III star formation. 

Furthermore, even if the ionization rate is high enough to make Hjj~ an important or dominant coolant at high densities, 
the effect of H^ cooling on the thermal evolution of the gas remains relatively small; the difference it makes to the temperature 
evolution at n > 10 s cm -3 is smaller than the error introduced by the uncertainty in the three-body H2 formation rate (reaction 
TBI). 

Finally, our model has also allowed us to explore the effects of cooling from the other minor ionic and molecular species 
present in the gas (e.g., H^~, H2D + , LiH, etc.). Despite making rather optimistic assumptions regarding the cooling from these 
species, we find that they are orders of magnitude less effective than Hjj~ at cooling high density gas, and hence are never 
significant. 
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Table Al. Chemical processes: collisional ionization (CI) 



No. Reaction Rate coefficient (cm 3 s 1 ) Ref. 

CI1 H + e" -> H+ + e~ + e" fe C ll = exp[-3.271396786 x 10 1 1 

+ 1.35365560 x 10 1 lnT e 

- 5.73932875 x 10°(lnT o ) 2 
+ 1.56315498 x 10°(lnT c ) 3 

- 2.87705600 x HrHlnTe) 4 
+ 3.48255977 x 10- 2 (lnT e ) 5 

- 2.63197617 x 10- 3 (lnT e ) 6 
+ 1.11954395 x KT 4 (lnT c ) 7 

- 2.03914985 x 10- 6 (lnT e ) 8 ] 

CI2 D + e- -» D+ + e" + e" k C n = fccii 1 
CI3 He + e" -> He+ + e" + e" fc CI3 = exp[-4.409864886 X 10 1 1 

+ 2.391596563 x 10 1 lnT c 

- 1.07532302 x lO^hi^) 2 
+ 3.05803875 x 10°(lnT o ) 3 

- 5.68511890 x lO-^lnTe) 4 
+ 6.79539123 x 10- 2 (lnT c ) 5 

- 5.00905610 x KT 3 (lnT e ) 6 
+ 2.06723616 x 10" 4 (lnT c ) 7 
-3.64916141 x 10- 6 (lnT c ) 8 ] 

CI4 Li + e" -> Li+ + e" + e" fc CI4 = 3.11 x 10- 8 T° 163 exp (-S^SS) 2 

Notes: T is the gas temperature in K, T3 = T/300 K, and T e is the gas temperature in eV. 



fNotes: 1 is tnc gas temperatu re m K, 1 3 = 1 /30 U K, an 
References: 1 - jjanev et~aH <|l987t l: 2 - IVoronovl (|l997ft 



Table A2. Chemical processes: photorecombination (PR) 



No. 



Reaction 



Rate coefficient (cm 3 s 1 ) 



Notes 



Ref. 



PR1 

PR2 
PR3 



H+ +e" 
D+ +e" 
He+ +e 



H + 7 
•D + 7 
-» He + 7 



14 ( 31^11 ) 1 ' 500 [ 1 . 0+ (iifsa) u ' 4ur ] 



fcpRi = 2.753 x 10 

fcpR2 = fcpRl 
fcpR3,rr,A = lO" 1 ^" ' 5 [12.72 - 1.615 log T 

- 0.3162(logT) 2 + 0.0493(logT) 3 ] 

fcpR3,rr,B = lO- 11 ^ - 5 [11.19 - 1.676 log T 

- 0.2852(log T) 2 + 0.04433(log T) 3 ] 
fc P R3,di = T- 15 [5.966 x lO" 4 exp (^f 

+ 1.613 x 10- 4 exp(- 555200 ' 
- 2.223 x 10" 5 ' 



455600 " 



PR4 Li+ + e - — > Li + 7 fcpR4,, 



1.036 x 10" 



exp I 



T > 
. 107.7 ; 



T 

-898200 

T 
-0.5 



)] 



1 -°+( l.l77 T xl0v ) t 



1.0- 

-1.388 



.177X10 7 

fc P R4,di = T- 1 - 5 [2.941 X 10~ 5 exp ( r 

+ 6.068 xlO- 5 — ^- 702400 



exp j, j 

7.753 x 10- 7 exp(^fl2a)] 



Case A 
Case B 



Notes: T is the gas temperature in K. Note that the recently revi sed valu e s for P R1 and for the 
radiative recombination portions of PR3 and PR4 presented by [Badnelj ( 120060) do not differ 
from the older rate coefficients quoted here by more than a couple of percent at the temperatures 

of interest in this study. 

References: 1 -iFerland et al.l dl992h; 2 -iHummer fc Storevl jl998T) ; 3 - iBadnelll ll2006at) ; 4 - 
IVerner fc Ferfandl l ll996l) : 5 - iBautista fc Badnelll d2007h 

APPENDIX A: CHEMICAL NETWORK 

In Tables IA1HA141 we list the chemical reactions included in our model of primordial gas, along with the rate coefficients 
adopted and the references from which these rate coefficients were taken. Some of these reactions are discussed in more detail 
in £13.11 In these tables, T is the gas temperature in K, T3 = T/300 K, and T c is the gas temperature in units of eV. 
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Table A3. Chemical processes: dissociative recombination (DR) 



No. 


Reaction 




Rate coemcient (cm° s 


Notes 


Ref. 


DR1 


H+ + e~ -» 


H + H 


fc DR1 = 1.0 X 10~ 8 

= 1.32 x io- 6 T-° 76 


T ^ 617 K 
T > 617 K 


1 


DR2 


HD+ + e- - 


-> H + D 


fc DR2 = 7.2 x icr 8 T-°- 5 




2 


DR3 


D+ +e- -* 


D + D 


fc DR3 = 3.4 x 10- 9 T 3 -°' 4 




3 


DR4 


H+ +e" -» 


H 2 + H 


fc D R4 = 2.34 x 10- 8 T 3 " ' 52 




4 


DR5 


H+ + e- -> 


H + H + H 


fc DR5 = 4.36 x 10- 8 T 3 " - 52 




4 


DR6 


H 2 D+ +e" 


-> H + H + D 


fcDR6 = 4.38 x 1G- 8 T 3 - - 5 




5 


DR7 


H 2 D+ +e" 


-»H 2 + D 


fc DR7 = 4.2 x 1G- 9 T 3 " - 5 




5 


DR8 


H 2 D+ +e" 


-» H + HD 


fc DR8 = 1.2 x 10- 8 T-°' 5 




5 


DR9 


HD+ + e- - 


-» D + D + H 


fc DR9 = 4.38 x 1CT 8 T 3 " - 5 




6 


DR10 


HD+ + e" - 


-» D 2 +H 


fc DR10 = 4.2 x Kr 9 T-°- 5 




6 


DR11 


HD+ + e" - 


-> HD + D 


fc D Rii = 1.2 x 10- 8 T-°- 5 




6 


DR12 


D+ +e~ -> 


D 2 + D 


fc D Ri 2 = 5.4 x 10- 9 T-°- 5 




7 


DR13 


D+ +e- -> 


D + D + D 


fc D Ri3 = 2.16 x io- 8 r 3 -°- 5 




7 


DR14 


HeH+ + e" 


-» Hc + H 


fc D Ri4 = 3.0 x 10- 8 T-°- 47 




8 


DR15 


HeD+ + e- 


-» He + D 


fc D Ri 5 = 3.0 x 10- 8 T-°- 47 




9 


DR16 


He+ + e~ - 


-> He + He 


fcDRie = 6.1 x 10- u T 3 -°- 9 




10 


DR17 


LiH+ + e- 


— > Li + H 


fc D Ri7 = 3.8 x 10- 7 T-°- 47 




11 


DR18 


LiD+ + e" 


-> Li + D 


fcDRis = 3.8 x 10- 7 T-°- 47 




12 


DR19 


LiH+ + e- 


-> Li + H 2 


fc D Ri9 = 1-6 x 10- 7 T-°- 5 




13 


DR20 


um + e - 


— > LiH + H 


fe D R20 = 2.0 x 10- 8 T-°- 5 




13 


DR21 


LiHJ" + e- 


-> Li + H + H 


fc D R 2 i = 2.0 x 10- 8 T-°- 5 




13 



Notes: T is the gas te mperature in K, a nd T3 
References: 1 ISchne ider et all ( 1994 

Walm slev, Flower, fc Pineau des Foreti 
1996h : 6 [Roberts. Herbst. fc Milieu 
19971 ): 8 - iGubermanl lll994h; 9~ ~ 



ICarata. Orel, fc Suzor-Weinerl lll999h : 
reaction: 13 — [Thomas et al.l l|i 




r/300 K. 



— w^u^ _ v^r ■■ « ^ ^ ^ ^1, , 

- iMcCall et al.l ll2004h; 5 - lLarsson et al. 


on Larsson et al.l 


199a): 7 — lLarsson et al. 


(Il998h. based on 


GubermarJ dl994h; 10 - 



(2006 



Stancil etafl jl998 

11 - iKrohn et alj feoplh : 12 - same as corresponding H 
C. Greene (private communication) 



S. C. 0. Glover & D. W. Savin 



Table A4. Chemical processes: charge transfer (CT) 



No. 



Reaction 



Rate coefficient (cm 3 s 1 ) 



Notes 



Ref. 



CT1 

CT2 
CT3 
CT4 
CT5 
CT6 
CT7 
CT8 
CT9 
CTfO 

CT11 
CTf2 
CTf3 
CTf4 
CTf5 
CTf6 
CTf7 
CTf8 
CTf9 



H + D+ -> D + H+ 



H + D" -► D + H" 
H + H+ ->■ H 2 + H+ 
H + HD+ — ► HD + H+ 
H + D+ -> D 2 + H+ 



H 
H 
H 
H 



He + H+ + 7 
He + He + H+ 
■+ LiH + H+ 
LiD+ -> LiD + H+ 



-He+ - 
•He+ - 
- LiH+ 



D + H+ -^H + D+ 

D + H" -^H + D" 
D + H+ -» H 2 + D+ 
D + HD+ -» HD + D+ 
D + D+^D 2 +D+ 
D + He+ -> He + D+ + 7 



fccTi = 2.06 x 10" 10 T - 396 cxp (-f ) 

+ 2.03 x 10" 9 T" - 332 
fc CT2 = 6.4 x fO- 9 T 3 - 41 
fc CT3 = 6.4 x 10" 10 
fc CT4 = 6.4 x 10" 10 
fc CT5 = 6.4 x 10" 10 

feCT6 

k CT7 = 1.0 x 10" 9 
fc CT8 = f .0 x 10" 11 cxp | 
fc CT9 = 1.0 x 10" 11 cxp | 

fcCTIO 



1.25 x 10- 15 T 3 °- 25 



67900 

T 
67900 



D + He^ 



He + He + D+ 



D + LiH+ -> LiH + D+ 
D + LiD+ -> LiD + D+ 
H 2 +H+ ^H + H+ 



fcCTll 
fcCTl 2 
fcCT13 
fcCT14 
fcCT15 
fcCT16 
fc CT17 

fcCT18 
fcCT19 



2.0 x 10" 10 T - 402 exp 
- 3.31 x io- 17 T 148 
6.4 x 10" 9 T 3 - 41 



6.4 x 10 



-10 



6.4 x 10" 10 

6.4 x 10" 10 

1.1 x 10" 15 T°- 25 

7.5 x 10 
1.0 x 10 



-10 






1-11 


exp | 


' 67900 


> T 


1-11 


exp | 


' 67900 


T 



} 



[-3.3232183 x 10" 7 
+ 3.3735382 x 10" 7 lnT 

- 1.4491368 x 10" 7 (lnT) 2 
+ 3.4172805 x 10" 8 (lnT) 3 

- 4.7813720 x 10" 9 (lnT) 4 
+ 3.9731542 x 10" 10 (lnT) 5 

- 1.8171411 x 10" n (lnT) 6 
+ 3.5311932 x 10" 13 (lnT) 7 ] 
xcxp( - 212 T 37 - 15 ) 



CT20 


H 2 + D+ -> 


D + H+ 


fcCT20 




^CT19 






CT21 


H 2 + He+ - 


* He + H+ 


fcCT21 




7.2 X 10~ 


15 




CT22 


H 2 + He+ - 


+ Hc + H + H+ 


&CT22 




3.7 x 10- 


14 cxp ( 


35\ 

t) 


CT23 


H 2 + Li+ -» 


Li + Hj 


fcCT23 




3.0 x 10" 


lOrp-1.5 
J 3 


cxp [- 


CT24 


HD + H+ - 


■H + HD+ 


&CT24 




&CT19 






CT25 


HD + D+ - 


• D + HD+ 






&CT19 






CT26 


HD + Hc+ - 


■+ He + HD+ 


&CT26 




7.2 x 10" 


15 




CT27 


HD + Hc+ - 


-> He + H+ + D 


&CT27 




1.85 x 10 


- 14 exp 




CT28 


HD + Hc+ - 


■+ He + H + D+ 


fcCT28 




1.85 x 10 


-14 

exp 


81 


CT29 


HD + Li+ ->• Li + HD+ 


fcCT 2 9 




&CT23 






CT30 


D 2 + H+ -> 


H + D+ 


fcCT30 




&CT19 






CT31 


D 2 + D+ -> 


D + D+ 


fccT31 




^CT19 






CT32 


D 2 + He+ - 


* He + D+ 


^CT32 




2.5 X 10" 


14 




CT33 


D 2 + He+ - 


* He + D+ + D 


fcCT33 




1.1 x 10" 


13^-0.24 
J 3 


CT34 


D 2 + Li+ -» 


Li + D+ 


fcCT34 




&CT23 






CT35 


He + H+ -> 


H + He+ 


fcCT35 




1.26 x 10 


- 9 T" - 75 exp(- 












4.0 x 10" 


37j>4.74 




CT36 


He + D+ -> 


D + Hc+ 


fcCT36 




&CT35 







T ^ 10000 K 
T > 10000 K 



2 
3 
4 
4 
5 
6 
7 
4 
1 

2 
4 
4 
4 
8 
9 
4 
4 
10 



4 
11 
11 
12 
4 
4 
4 
13 
13 
4 
4 
4 
14 
14 
4 
15 



Is H% cooling important? 



Table A4 - continued 



No. 



Reaction 



Rate coefficient (cm 3 s 1 ) 



Notes Ref. 



CT37 
CT38 
CT39 
CT40 
CT41 
CT42 
CT43 
CT44 
CT45 
CT46 
CT47 



Li + 
Li + 
Li + 
Li + 
Li + 
Li + 
Li + 
LiH 
LiH 
LiD 
LiD 



H+ - 

H+ - 

D+ - 

D+ - 

H 2 + - 
HD+ 

D 2 + " 
+ H+ 
+ D+ 
+ H+ 
+ D+ 



H- 
H- 
D + 
D + 
H 2 



- Li+ 

- Li+ - 
-Li+ 
-Li+ - 

-Li+ 
-> HD + Li+ 
D 2 + Li+ 
-»H+ LiH+ 
-> D + LiH+ 
-> H + LiD+ 
-> D + LiD+ 



7 



7 



fcCT38 
fcCT39 
fcCT40 
fcCT41 
^CT42 
fcCT43 
fcCT44 
&CT45 
&CT46 
fccT47 



: 2.5 X 10 

: 1.7 X 10 

: 8.0 X 10 

: 1.1 X 10 

: 3.0 X 10 

: kcT41 
: fcCT41 

: 2.0 x 10" 15 
: 2.0 x 10" 15 
: 2.0 x 10~ 15 
: 2.0 x 10~ 15 



1 exp y l91n 
13y-0.051 



22 T|- 8 exp 
13^-0.051 

IOt^-1.5 



CX P ( 282000 ) 



T 
1800, 



exp (- 



r) 



16 
17 
18 
19 
20 

4 

1 

21 
21 
21 
21 



Notes: T is the gas temperat ure i n K, and T3 = T/300 K 



References: 1 - Savin (20021) ; 2 -iDalearno fe McDowelll dl956h . scaled by D re duced mass; 3 



1989); 



iKarpas. Anicich. fe Huntress! (Tl979f) : 4 — sam e as corresp onding H rea c tion; 5 
Estimate by IStancil et ail (Il998l) . based on IStancil et all dl993h 



1996); 8 



IZvgelman et al 
7 [§ 



Savin et al 



11 



6 - Estimate by IStancil et ail (1998), based on IStancil et ail dl993f ): 7 - IStancil et al . 
Zygelman et al.| |l989l) , scaled by D red uced mass; 9 — As ref. 6, but scaled by D reduced 
Barlowf(ll984h : 12 - Estimate, based o n low-energy e xtrapola- 
(11997 13 - t otal rat e coefficient fromlBarlowl <1984h. branch- 
(|l989h ; 14 - IWalmsIev, Flower, fc Pineau des Foretsl 
16 - iKimura. Dutta fc Shimakural lll994l) ; 17 



ll2004h~ 



mass; 10 

tion of cross-section in Wutteetal^ 
ing ra tios from iPineau des Forets et al 

1120041); 15 - IKimura et ail dl993h; .... 

IStancil fc Zygelman! Jl996h; 18 - IKimura, Dutta. fc Shimakural dl994l) , scaled by D reduced 
mass: 19 - IStancil fc Zygelmanl dl996l). scaled by D reduced mass; 20 - From detailed balance 



applied to inverse reaction; 21 - iBodo et al. I feOOll) 



S. C. 0. Glover & D. W. Savin 

Table A5. Chemical processes: radiative attachment and radiative association (RA) 



No. 



Reaction 



Rate coefficient (cm 3 s 1 ) 



Notes 



Ref. 



RA1 



fcrtAi 



: dcx[ 



RA2 
RA3 

RA4 
RA5 



RA6 
RA7 
RA8 
RA9 
RA10 

RA11 
RA12 
RA13 
RA14 
RA15 
RA16 
RA17 
RA18 
RA19 
RA20 
RA21 
RA22 
RA23 
RA24 
RA25 
RA26 
RA27 
RA28 
RA29 
RA30 
RA31 



D 
H 

H 
H 



H 
H 
H 
H 
H 



H 



+ 



■7 
-7 



D+ -+ HD+ + 7 
D -> HD + 7 



"2 

HD+ 

He+ 
Li+ - 



D 
D 

H 2 4 
H 2 -t 
Li+ 
HD - 
HD - 
D 2 4 
D 2 -4 



He + H+ 
He + D+ - 
He + Hc+ 
Li + e - — 
Li + H+ - 
Li + D+ - 
Li + H-> 



Hj+7 
-> H 2 D+ +7 

HD+ + 7 
-» HeH+ + 7 
► LiH+ + 7 



D + H+ — i 
D + D+ - 
D + H+^ 
D + HD+ 
D + D+ - 
Hc+ - 
■Li+ - 
4-H+- 



-7 



HD 4 

D++7 
H 2 D+ +7 
-» HD^ 



d; 



2 T7 
3 t7 
HeD+ + 7 
LiD 4 
Hj 



-7 



D+ - 
-H 2 - 
-H+ 

-D+ 
H+ - 
D+ - 



-3+7 
H 2 D+ +7 
. LiH+ + 7 
* H 2 D+ +7 
> HD^ 



HD 



+ 



f7 
-7 



D3++7 



-+ HeH+ + 7 
-> HeD+ + 7 
-+ He+ + 7 
■ Li - + 7 
-> LiH+ + 7 
-» LiD+ + 7 
LiH + 7 



-17.845 + 0.762 log T 
+ 0.1523(log T) 2 

- 0.03274(logT) 3 ] 

= dex[-16.420 + 0.1998(log T) 2 

- 5.447 x 10~ 3 (logT) 4 
+ 4.0415 x 10- 5 (logT) 6 ] 

fc R A 2 = ^RAl 

fc RA3 = dex[-19.38 - 1.523 log T 

+ 1.118(logT) 2 - 0.1269(logT) 3 ] 

fc RA4 = 3.9 x 10- 19 T 3 18 exp (f?) 

fc RA5 = 10~ 25 [2.80202 - 6.63697 In T 

+ 4.75619(lnT) 2 - 1.39325(lnT) 3 
+ 0.178259(lnT) 4 - 0.00817097(ln T) 5 ] 
= 10- 25 exp [507.207 -370.889 In T 
+ 104.854(lnT) 2 - 14.4192(lnT) 3 
+ 0.971469(lnT) 4 - 0.0258076(lnT) 5 ] 

fc RA6 = 1.5 x 10- 17 T 3 18 exp I 

fe RA 7 = 1.2 x 10~ 17 T 3 1 ' 8 exp I 

fc R A8 = 1-1 x 10- 17 Tj- 8 exp I 

fe RA9 = 4.16 x 10- 16 T 3 " - 37 exp 1 

fc R .Aio = dcx [-22.4 + 0.999 log T 
-0.351(log T) 2 ] 

fcRAXl = 3.9 X 10" 

feRA12 = 1.9 X 10- 

fc R M3 = 7-0 x lO-^'i^ejci, ( 

fc R A14 = 5.2 x 10" 

fc R A15 = 4.3 x 10" 



' 201 

. T / 
'20^ 

* T I 
'20_S 

v T , 



' ( 87600 ) 



- 19 T 3 1 - 8 exp 

- 19 T 3 L8 ex P 
-18711. 

- 18 2*- 8 exp 



-16n 



RA32 Li + D -> LiD + 7 



fc R A16 = 5.0 x 10 
fe R A17 = 1.5 x 10 
fc RA 18 = 1.0 X 10- 
fc RA 19 = 1.0 X 10" 
fcRA20 = 1.0 X 10- 
feRA21 = 1.0 X 10- 
fcRA22 = 1.0 X 10" 

fc RA23 = 1-0 x 10- 

feRA24 = 1.0 X 10- 
fe R A 2 5 = 8.0 x 10- 
fc R A26 = 1.0 X 10" 
fe R A27 = 4.76 X 10 
fc R A28 = 6.1 X 10" 
fe R A 2 9 = 4.8 x 10- 
fc R A30 = 6.4 x 10- 
feRASl = 10- 20 {3.22 

+ 6.3 x 10- 
fc R A32 = 5.5 x 10- 20 : 



-0.37 



'20> 
T , 
'»> 
T , 
'2CP 
T . 
'20> 
T . 
'201 
. T , 



exp 
'exp ( 



- 22 Q 

-16 
-16 
-22 
-16 
-16 
-16 
-16 

-20 T -0.24 exp 
-19 T -0.24 exp 

- 2 »T 3 L82 exp 
17 T 0.58 exp 
14 J.-0.49 
14J1-0.49 



( 87600 ) 
7000 / 



T 
" 4000 
T 
4000 



29 N 

T , 



( 17 2 Oo) 



+ [0.0657(T/1000)- : 
^T/IOOO) - 837 ]" 1 } 



T < 6000 K 



T > 6000 K 



10 < T 4. 200 K 



T > 200 K 



3 
3 
6 
6 
6 
6 
9 

10 
11 
12 
11 
11 
11 
11 
13 
6 

11 
15 
16 
6 
17 



exp I 



T 

' 3300 / 



Notes: T is the g as temperature in K, and T 3 = T/300 K . 

Ref erences: 1 - IWishartl lll979l) : 2 - iRamaker fe Peekl lll976h : 3 - IRamaker fe Peek! il976ri 
and Frommhold & Pickett al978f) ■ scaled by D reduced mass; 4 - iDickinsonl l l2005l) ; 5 - 
iDalgarno fc Mc Dowelll jl956T); 6 — Same a s cor responding H reaction , but scaled by D reduced 
mass; 7 - |Kraem er 1 Spirko 1 fc Jufeld jl995l) ; 8 - IDalgarno et ah I lll996f) :lGianturco &: Gori G iorgi 
Jl996h ; 9 - IStancil et al.l (Il996l'), scaled by D reduced mass; 10 - iGerlich fc Horning! l ll992l) ; 11 - 
esti mate, based on iGerh chfc Horningl l ll992j): highly uncertain ; 12 — estimate - see also H3.1.11I 



13 
16 



. Spirko. fc Kraemerl jlQQSl'); 14 -IStancil et alj jl993h : 15 - iRamsbottom et alj Jl994h : 
Dalgarno et al.l lll996h : 17 iBennett et al.l ( l2003h 



Is H% cooling important? 

Table A6. Chemical processes: associative detachment, dissociative attachment and associative 
ionization (AD) 



No. Reaction Rate coefficient (cm 3 s 1 ) Notes Ref. 



AD1 


H + 


H" -* 


H 2 + e" 


fcADl = 


= 1.5 x 10- 9 T 3 " 01 






1 


AD 2 


D + 


H~ — » 


HD + e~ 


&AD2 = 


= 1.5 x 10- 9 T 3 " 01 






2 


AD3 


H + 


D" -> 


HD 4- e — 


^AD3 - 


= 1.5 x 10- 9 T 3 " 01 






2 


AD4 


D + 


D" -> 


D 2 + e - 




= 1.6 x 10- 9 T" 01 






2 


AD 5 


H 2 - 


f e - -> 


H + H- 


^AD5 - 


= 2.7xlO- 8 T- 1 - 27 exp( 43 ° 00 ) 




3 


AD6 


HD 


+ e" - 


+ H + D- 


^AD6 = 


= 1.35 x 10 _9 T -1 - 27 exp (- 


43000 ^ 
T J 




4 


AD7 


HD 


+ e" - 


+ D + H~ 


&AD7 - 


= 1.35 x 10 _9 T -1 - 27 exp (- 


43000 \ 
T J 




4 


AD8 


D 2 - 


f e" -> 


D + D _ 


&AD8 - 


= 6.7 x 10 _11 T _1 - 27 exp (- 


43000 S 
T ) 




1 


AD9 


H+ 


+ H- - 


- H+ + e" 


^AD9 = 


= 6.9 x io~ 9 T-°' 35 
= 9.6 x i -7 T -0-90 




T < 8000 K 
T > 8000 K 


5 


AD 10 


H+ 


+ D" - 


-+ HD+ + e" 


&AD10 


= 1.1 x 10- 9 T-°- 4 






2 


AD11 


D+ 


+ H- - 


-» HD+ + e~ 


feADll 


= 1.1 x 10- 9 T-°- 4 






2 


AD12 


D+ 


+ D- - 


- D 2 + e ~ 


^AD12 


= 1.3 x 10- 9 T 3 "°- 4 






2 


AD 13 


H* 


+ H" - 


-H 3 + + e- 


&AD13 


= 2.7x 10- 10 T 3 - a485 e X p| 


' T \ 




6 


S 31200/ 


! 


AD 14 


H 2 + 


+ D- - 


-> H 2 D+ +e" 


&AD14 


= 2.24 x 10- 10 T 3 -°- 49 exp 1 


> T i 




6 


I 43600 , 




AD15 


HD+ +H" 


-> H 2 D+ +e" 


&AD15 


= 2.9 x 10- 10 T-°' 485 c X p( 


f T \ 




2 


, 31200 1 




AD16 


HD 


^+D" 


-> HD+ +c- 


&AD16 


= 3.7x 10- 10 T-°' 485 e X p( 


> T \ 




2 


,31200/ 




AD17 


D 2 + 


+ H~ - 


-» HD+ + e" 


A; AD 17 


= 3.0 x 10- 10 T- ' 485 exp( 


' T \ 




2 


,31200 / 




AD18 


D 2 + 


+ D~ - 


- D+ + e- 


&AD18 


= 3.9 x 10- 10 T- a485 exp( 


' T \ 




2 


, 31200/ 




AD19 


Li 4 


•H- -» 


LiH + e" 


&AD19 


= 4.0 x 10" 10 






7 


AD20 


Li 4 


•D- -t 


LiD + e~ 


&AD20 


= 4.0 x 10" 10 






2 


AD21 


Li- 


+ H -» 


LiH + e" 


&AD21 


= 4.0 x 10- 10 






7 


AD22 


Lr 


+ D -> 


LiD + e" 


&AD22 


= 4.0 x lO" 10 






2 



Notes: T is the g as temperatur e in K, and T3 = T/300 K. 

References: 1 - lLauna v et al. l d 1 99 lh : 2 — S ame as corresponding H rea ction, but scaled by 
D red uced mass; 3 - ISchulz fc Asundil (1196711: 4 - IXu fc Fabrikantj fcOOll l: 5 - iPoulaert et all 
jl978h ; 6 - iNaii et al.l dl998l) ; 7 - IStancil et all ( 11999) 



S. C. 0. Glover & D. W. Savin 

Table A7. Chemical processes: collisional detachment and collisional dissociation (CD) 



No. 



Reaction 



Rate coefficient (cm 3 s 1 ) 



Notes 



Ref. 



CD1 



CD2 



H -+ H + H- 



fe CD1 = exp[-1.801849334 x 10 1 
+ 2.36085220 x 10° lnT c 

- 2.82744300 x lO-^lnT,,) 2 
+ 1.62331664 x 10- 2 (lnT e ) 3 

- 3.36501203 x 10- 2 (lnT o ) 4 
+ 1.17832978 x 10- 2 (lnT e ) 5 

- 1.65619470 x 10- 3 (lnT e ) 6 
+ 1.06827520 x 10~ 4 (lnT e ) 7 

- 2.63128581 x 10~ 6 (lnT c ) 8 ] 
fc C D2 = 2.5634 x io-9r o 178186 

= exp[-2.0372609 x 10 1 
+ 1.13944933 x 10° lnT e 

- 1.4210135 x lO-^lnTe) 2 
+ 8.4644554 x 10- 3 (lnT c ) 3 

- 1.4327641 x 10- 3 (lnT c ) 4 
+ 2.0122503 x 10" 4 (lnT c ) 5 
+ 8.6639632 x 10~ 5 (lnT e ) 6 

- 2.5850097 x 10- 5 (lnT c ) 7 
+ 2.4555012 x 10- 6 (lnT c ) 8 

- 8.0683825 x 10- 8 (lnT e ) 9 ] 



CD3 


H~ 


+ D^H + D + c~ 


fcCD3 


= &CD2 




CD4 


H~ 


+ He — > H + He + c - 


fcCD4 


= 4.1 X 10" 17 T 2 exp (- 


19870 ^ 
T ) 


CD5 


D" 


+ e~ — ► D + e~ +e~ 


fcCD5 


= fccDl 




CD6 


D" 


+ H^D + H + e" 


fcCD6 


= fcCD2 




CD7 


D" 


+ D^D + D + c- 


fcCD7 


= fccD2 




CD8 


D" 


+ He — > D + He + e - 


fcCD8 


= 1.5 X KT 17 T 2 exp (- 


19870 ^ 


CD9 


H 2 + H^H + H + H 


fcCD9 


= 6.67 x 10- 12 T 5 exp 


[ a + 63 ™ 3 



CD10 H 2 + H 2 ^H + H + H 2 
CD11 H 2 +He ^ H + H + He 



= k TB i/K 

, _ 5.996xl0" 30 T 41881 

KCD10 - n nxK 7 f;ivin-fi T )B.fi^l ex P 



CD12 H 2 +e" 



H + H- 



CD13 HD + H^H + D + H 

CD14 HD + H 2 -> H + D + H 2 

CD15 HD + He -> H + D + He 

CD16 HD + e- -> H + D + e" 

CD17 D 2 +H^D + D + H 

CD18 D 2 + H 2 -* D + D + H 2 

CD19 D 2 + He -> D + D + He 

CD20 D 2 + e" ^ D + D + e- 



fccmi 

fccD12 

&CD13 
^CD14 
&CD15 : 
fcCD16 

&CD17 
fcCD18 
fcCD19 
fcCD20 



(1.0+6.761x10 
kTB2/K 

dex [-27.029 + 3.801 log T — 
6.6 x 10- 10 T°- 115 exp(- 52000) 
4.49xl0- 9 T°- 11 exp(-i2^ 8 ) 
1.91 x 10- 9 T - 136 exp 

&CD9 
&CD10 
feCDll 

5.09 x 10- 9 T 0128 exp( 
1.04 X 10~ 9 T - 218 exp ( 

&CD10 
fecDll 

8.24 X 10- 9 T°- 126 exp 
2.75 x 10- 9 T°- 163 exp 



t . 

53070.7 



T 

5333P. 



CD21 


LiH+ - 


ID -> 


Li+ + H + D 


fcCD21 


= 1.0 x 10" 


-Q 

" exp 


CD22 


LiH+ - 


ID -+ 


Li + H+ + D 


&CD22 


= 1.0 x 10" 


* exp 


CD23 


LiH+ - 


ID -> 


Li + H + D+ 


fcCD23 


= 1.0 x 10" 


9 

* exp 


CD24 


LiD+ - 


1D^ 


Li+ + D + D 


&CD24 


= 1.0 x 10" 


9 

exp 


CD25 


LiD+ - 


FD^ 


Li + D+ + D 


kCT>25 


= 1.0 x 10" 


9 

■' exp 


CD26 


LiH+- 


r-H 2 - 


-> Li+ + H 2 + H 2 


kcT>26 


= 1.0 x 10" 


9 exp 



1400 

T ) 
97500 ^ 
T I 
97500 \ 
T ) 
1400 \ 
T ) 
97500 'i 
T ) 
3000 \ 
T ) 



T e < 0.1 eV 
T c > 0.1 eV 



v=0 
LTE 
v = 
LTE 
v = 
LTE 
v = 
LTE 

See oHiXol 
See oHiXol 
See oHiXol 

v = 
LTE 



v = 
LTE 



2 
3 
2 
2 
2 
4 
5 
6 
7 
6 
8 
9 

10 
10 
2 
2 
2 

11 
11 

2 
2 
2 
10 
10 
12 
12 
12 
12 
12 
13 



Notes: T is the gas temperature in K and T e is the gas temperature in eV. K is the equilibrium 
constant relating reactions TBI and CD9, and reactions TB2 and CD10; its value is given in 
33X51 

References: 1 — iJanev et al. 

I dl987T) ; 2 — Assumed same as corresponding H reaction; 3 — 
Huq et al l dl982f): 4 - S ame as corresponding H reaction, but scaled by D reduced mass; 5 - 
Mac Low &: ShulU lll986|); 6 — determined f rom t hree- body rate coefficie nt by detailed balance 
(see fl3Tt~7ll: 7 -iMartin. Keogh. fc Mandvl l|l998h : 8 - iDove et all (|l9S7h : 9 - determined from 
the IWa lkauslras"~fclKaufman ni975l) rate coefficient for reaction TB3 by detailed balance; 10 - 
iTrevisan fc Tennyson! 1 12002311 ; 11 - iTrevisan fc Tennyson! d2002bh ; 13 - estimate - see also &13.1.11I 



Is H% cooling important? 53 



Table A8. Chemical processes: mutual neutralization (MN) 



No. 



Reaction 



Rate coefficient (cm 3 s 



Ref. 



MN1 

MN2 

MN3 

MN4 

MN5 

MN6 

MN7 

MN8 

MN9 

MN10 

MN11 

MN12 

MN13 

MN14 

MN15 

MN16 

MN17 

MN18 

MN19 

MN20 

MN21 

MN22 

MN23 

MN24 

MN25 

MN26 

MN27 

MN28 

MN29 

MN30 

MN31 

MN32 

MN33 

MN34 

MN35 

MN36 

MN37 

MN38 

MN39 

MN40 

MN41 

MN42 

MN43 

MN44 



H+ + H" 
D+ + H- 
H+ + D- 

D+ + D" 



H 2 + - 
H+- 



■H- 
-D- 
H+ + D- 
HD+ + H" 
HD+ + H" 
HD+ + D" 
HD+ + D" 



H 



D 



H+ 



H 2 D+ 
H 2 D+ 
H 2 D+ 
H 2 D+ 
H 2 D+ 
H 2 D+ 
H 2 D+ 
HD+ - 
HD+ - 
HD+ - 
HD+ - 
HD+- 
HD+ - 
HD+ - 



f H" 
f H" 
f H" 
f D" 
f D" 
f D" 
f D" 
H~ 
H~ 
H" 
H _ 
D~ 
D~ 
D~ 



D 3 + " 

D i- 

He+ 
He+ 
Li+- 



H" 
H~ 
D" 
D" 
f H" 
f D" 
- H" 



Li+ + D" 
Li" + H+ 
Li" + D+ 



H + H 
D + H 
D + H 
D + D 
H 2 +H 
H + H + 
H 2 + D 
H + H + D 
-> HD + H 
-» D + H + H 
-» HD + D 
-^D + H + D 
D 2 + H 
D + D + H 
ID 
D + D 
f H + H 

f H 2 
f H + D 
f HD 

H 2 + H + D 



D 2 
■ D4 
H 2 
H 2 
H 2 
H 2 



-> H 2 - 
-* HD 
-+H 2 - 
-* H 2 - 
-> HD 
-> HD 
H 2 + 
HD 4 
HD 4 
D 2 + 
HD 4 
HD 4 



-HD 
f H + H 
-D + D 
-D 2 

f H + D 
f HD 
D 2 

H + D 
HD 
H + H 
D + D 
D 2 



-^D 2 +H + D 
HD + D 2 
D 2 + H + D 
D 2 + D + D 
D 2 + D 2 

^Hc + H 

4 He + D 

• Li + H 

• Li + D 
•Li + H 
> Li + D 



^MNl = 
fcMN2 = 
^MN3 = 
^MN4 = 
&MN5 = 
%N6 = 
&MN7 = 
^MN8 = 
fcMN9 = 
&MN10 
^MNll 
fcMNl 2 
&MN13 
fcMN14 
^MN15 
kMN16 
^MN17 
fcMN18 
^MN19 
fcMN 2 
^MN21 
&MN22 
^MN 2 3 
^MN24 
*MN25 
fcMN 2 6 
^MN27 
?£MN28 
fclVlN29 
&MN30 
^MN31 
^MN32 
^MN33 
^MN34 
^MN35 
fcMN36 
^MN37 
fcMN38 
^MN39 
^MN40 
fcMN41 
^MN42 
^MN43 
^MN44 



2.4 x lCT 6 ! 1 - 1 / 2 (l.O + 5.0 x 10~ 5 T) 

1.1 X fc M Nl 
1.1 X fc MN1 

1.3 x A;mni 



1.4 x 10 
1.7 x 10 
1.7 x 10" 7 t!: 

1.5 x 10~ 7 T 3 ~ 



-0.5 
L 3 

-7T-0.5 
3 

-7rp-0.5 
"-0.S 



7rp-0.5 

3 

7 T -0.5 
3 

7^-0.5 
3 

7rp-0.5 
3 

7 T -0.5 

J 3 
7 T -0.5 



-in 



1.5 x 10 
1.9 x 10 
1.9 x 10 
1.5 x 10 

1.5 x 10 
2.0 x 10 

2.0 x 10 
2.3 x 10 

2.3 x 10 
2.9 x 10 
2.9 x 10 

1.6 x 10 
1.6 x 10 
1.6 x 10 
1.5 x 10 
1.5 x 10 
1.5 x 10 

1.5 x 10 
1.2 x 10 
1.2 x 10 
1.2 x 10 

1.2 x 10 

2.1 x 10 
2.1 x 10 
2.1 x 10 

2.4 x 10 
2.4 x 10 

3.3 x 10 
3.3 x 10 
2.32 x 10- 7 T. 
3.03 x 10~ 7 T. 
2.93 x 10 

2.06 x 10~ 7 T. 
1.8 x 10 



0.5 
-0.5 



0.5 



T T -0.5 
3 

7ti-0.5 

J 3 
7 T -0.5 
3 

7 T -0.5 

J 3 
7T-0.5 
3 

7 T -0.5 
3 

7 T -0.5 

J 3 
7 T -0.5 

1 3 
7 T -0.5 

J 3 
7 T -0.5 
3 

7 T -0.5 
3 

7^-0.5 

J 3 
7^-0.5 
3 

7ti-0.5 

J 3 
7^-0.5 
3 

7t,-0.5 
3 

7^-0.5 

J 3 
7rp-0.5 
3 

7rp-0.5 



-in 



0.5 



Tt-0.6 



°- S2 expC T \ 

3 CX P I 22400 J 

3 exp V 22400 ) 
0.477 / T 



7 T - 
^3 
7T-0 



eX P V 23200 / 
T \ 



■ 2.06 x 10 



3 eX P V 18300 ) 
7 T -0.477 / T 

± 3 ex fV 23200 

-7 T -0.5 / T \ 
± 3 LX P \ i 8 300 , 



Notes: T is the gas temperature in K, and T3 = T/300 K. Some of the mutual neutralization 
reactions listed her e also includ e disso ciation or transfer in the process. 

References: 1 — ICroft et al, I jl999h: 2 - S ame as cor responding H reaction , but scaled 
by D reduced mass ; 3 - iDalgarno fc Leppl jl987h ; 4 - iDalearno fc McDowell dl956h ; 5 - 
iLe Teuff et al. (2000); 6 - As 2, with the additional assumption of equally probable outcomes; 7 
- IPeart fc Havtonl l|l994h 
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Table A9. Chemical processes: three-body association (TB) 



No. 


Reaction 






Rate coefficient (cm 6 s 


- 1 ) 




Ref. 


TBI 


H + 


H4-H 


^H 2 +H 


See §3.1.7| 








TB2 


H + 


H4H 2 -» 


H 2 + H 2 


See §3.1.7| 








TB3 


H + 


H4-Hc -> 


H 2 +Hc 


&TB3 = 


= 6.9 x 10" 32 T- 


0.4 




1 


TB4 


H + 


D + H 


^HD + H 


See §3.1.7| 








TB5 


H + 


D+H; 




HD + H 2 


See §3.1.7| 








TB6 


H + 


D 4- He — > 


HD + Hc 


&TB6 - 


= 6.9 x 10" 32 T- 


0.4 




2 


TB7 


D + 


D + H 


^D 2 +H 


See §3.1.7| 








TBS 


D + 


D + H 2 ^ 


D 2 + H 2 


See §3.1.7] 








TB9 


D + 


D4-Hc 


D 2 + He 


&TB9 = 


= 6.9 x 10" 32 T- 


0.4 




2 


TB10 


H+ 


4-H + H -> 


H++H 


^TBIO 


= 1.203 x 1CT 29 


T~ 


1.041 


3 


TB11 


D+ 


+ H + H^ 


HD+ + H 


&TB11 


= 1.203 x 10 -29 


T~ 


1.041 


2 


TB12 


H+ + D + H — > 


HD+ + H 


&TB12 


= 1.203 x 10~ 29 


T~ 


1.041 


2 


TB13 


D+ 


+ D4H^ 


D++H 


^TB13 


= 1.203 x 10~ 29 


T~ 


1.041 


2 


TB14 


H+ 


+ H 2 + 


H - 


^H+ + H 


KTB14 


1 n \y in — 28 
= l.U X 1U 






4 


TB15 


H+ 


+ H 2 + 


H 2 


— H+ + H 2 


fcTB15 


= 5.4 x 10~ 29 






5 


TB16 


H+ 


+ H 2 + 


He 


-> H+ + He 


^TB16 


= 1.07 x 10" 28 






5 


TB17 


D+ 


4-H 2 + 


H - 


-> H 2 D+ + H 




= 1.0 x 10" 28 






4 


TB18 


D+ 


+ H 2 + 


H 2 


-> H 2 D+ + H 2 


fcTB18 


= 5.4 x 10~ 29 






2 


TB19 


D+ 


4-H 2 + 


He 


-> H 2 D+ + He 


&TB19 


= 1.07 x 10~ 28 






2 


TB20 


H+ 


4-HD4 


H 


-+ H 2 D+ +H 


^TB20 


= 1.0 x 10~ 28 






4 


TB21 


H+ 


4-HD4 


H 2 


-> H 2 D+ + H 2 


&TB21 


= 5.4 x 10" 29 






2 


TB22 


H+ 


4-HD4 


He 


-+ H 2 D+ +He 


^TB22 


= 1.07 x 10~ 28 






2 


TB23 


D+ 


4-HD4 


-H 


-> HD+ + H 


&TB23 


= 1.0 x 10" 28 






4 


TB24 


D+ 


4-HD4 


•H a 


HD+ + H 2 


^TB24 


= 5.4 x 10" 29 






2 


TB25 


D+ 


4-HD4 


-He 


-» HD+ + He 


^TB25 


= 1.07 x 10~ 28 






2 


TB26 


H+ 


4-D 2 + 


H - 


■+ HD+ + H 


&TB26 


= 1.0 x 10" 28 






4 


TB27 


H+ 


4-D 2 + 


H 2 


-> HD+ + H 2 


^TB27 


= 5.4 x 10~ 29 






2 


TB28 


H+ 


+ D 2 + 


He 


-> HD+ + He 


&TB28 


= 1.07 x 10~ 28 






2 


TB29 


D+ 


+ D 2 + 


H - 


^D+ + H 


&TB29 


= 1.0 x 10" 28 






i 


TB30 


D+ 


+ D 2 + 


H 2 


-> D+ + H 2 


^TB30 


= 5.4 x 10~ 29 






2 


TB31 


D+ 


+ D 2 + 


He 


-> D+ + He 


&TB31 


= 1.07 x 10~ 28 






2 


TB32 


Li 4 


-H + H 


-+ LiH + H 


^TB32 


= 2.5 x 10" 29 T 


-1 




6 


TB33 


Li 4 


H + H; 


2 — * 


LiH + H 2 


&TB33 


= 4.1 x 10" 30 T 


-1 




6 


TB34 


Li 4 


-D + H 


-> LiD + H 


^TB34 


= 2.5 x 10~ 29 T 


-1 




2 


TB35 


Li 4 


- D + H 2 -> 


LiD + H 2 


&TB35 


= 4.1 x 10~ 30 T 


-1 




2 



Notes: T is the ga s temperature in K. 

Reference s: 1 — [ Wa lkauskas & Kaufman (1975J); 2 — Sa me as corresponding H re action; 
3 - iKrstic. Janev. fc Schultj ll2003r i; 4 - Estimate; 5 - iGerlich fc Horning! dl992t) ; 6 - 
iMizusawa. Omukai. fc Nishil i hooa i 
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Table A10. Chemical processes: isotopic exchange (IX) 



No. 


Reaction 




Rate coefficient (cm 3 s 


- 1 ) 


1X1 


H+ + D - 


* HD+ +H 


fc Ixl = 1.07 x 10- 9 T« 


° 62 -p(-23loo 


1X2 


H+ + D - 


* HD + H+ 


fc IX2 = 1.0 x 10~ 9 




1X3 


HD+ + H 


+D 


fc IX3 = 1.0 x 10~ 9 exp 


(-¥) 


1X4 


HD+ + H 


-> H 2 + D+ 


fc IX4 = 1.0 x 10~ 9 




1X5 


HD+ + D 


D+ +H 


fe IXS = 1.0 x 10~ 9 




1X6 


HD+ + D 


^D 2 +H+ 


fc IX6 = 1.0 x 10~ 9 




1X7 


D++H- 


* HD+ +D 


fc IX7 = 1.0 x 10~ 9 exp 




1X8 


D++H- 


* HD + D+ 


fc IX8 = 1.0 x 10" 9 




1X9 


H 2 + D+ 


-> HD + H+ 


fci X9 = 4.17xl0- 10 + 


8.46 x 10" 10 lo 








- 1.37 x 10~ 10 


(logT) 2 


1X10 


H 2 + D+ 


-► HD+ + H 


feixio = [1.04 x 10- 9 - 


f 9.52 x 10- 9 (- 



Notes 



1X11 
1X12 
1X13 
1X14 



1X21 
1X22 
1X23 
1X24 
1X25 
1X26 
1X27 
1X28 
1X29 
1X30 
1X31 
1X32 
1X33 
1X34 
1X35 



HD + H+ 
HD + H+ 
HD + D+ 
HD + D+ 



■ H 2 + D+ 
H+ +D 
D 2 +H+ 



D^ 



H 



1X15 D 2 + H+ -> HD + D+ 
1X16 D 2 + H+ -► HD+ + D 

1X17 H 2 + D -> HD + H 



1X18 HD + H -> H 2 + D 

1X19 HD + D -> D 2 + H 
1X20 D 2 + H ->• HD + D 



H 



3 

H 2 D 



H+ +D 



D -» H 2 D 
+ H - 

H 2 D+ 
HD+ + H -> H 2 D+ 



D 



HD^ 



D D 3 
I -> HD^ 



H 

- D 



H+ + HD^H 2 D++H 2 



H+ + D 2 
H+ + D 2 
H 2 D+ 
H 2 D+ 



H 2 - 
HD 



H 2 D+ 
HD+4 
H+ 



Hj 



f- HD 
H 2 
^HD 
+ D 2 



H 2 D+ + HD -► HD+ + H 2 
H 2 D+ +D 2 -> HD+ +HD 
H 2 D+ +D 2 -> D+ +H 2 



&IX12 
fclX13 
^1X14 



&IX15 
^1X16 

&IX17 



^1X19 
&IX20 



HD^ 



H 2 -» 



D 2 



felX21 
&IX22 
&IX23 
^1X24 
fclX25 
&IX26 
fclX27 
fclX28 
&IX29 
&IX30 
^1X31 
^1X32 
&IX33 
&IX34 
&IX35 



1.81xl0-»(^)^ exp(-2fHl) 



1.1 X 10- 9 exp (-4P 



lOOOO / 

21000 i 



T 
21600 
T 



cxp 



1.0 X 10" 9 exp 

1.0 x 10~ 9 

[3.54 x 10~ 9 + 7.50 x 10- 10 
-2.92xlO- 1 «( T?5 2 n 

2.1 x 10- 9 cxp(-i|!-) 
[5. 18 x 10- 11 + 3.05 x 10~ 9 ( T 

_ 5 .42xl0- 10 (^ 2 ' 

dex [-56.4737 + 5.88886 log T 
+ 7.19692(log T) 2 
+ 2.25069(logT) 3 
- 2.16903(log T) 4 
+ 0.317887(log T) 5 1 



(-21100) 



exp | 



10000 / 

20100 1 



5.25 x 10 



-10 


exp 


' 5207 


s T 


-11 


exp 


1 4430 


s T 


-11 


exp 


1 4430 


s T 


-11 


exp 


' 3220 


v T 



173900 ~> 



dex [-86.1558 + 4.53978 logT 
+ 33.5707(log T) 2 
- 13.0449(logT) 3 
+ 1.22017(log T) 4 
+ 0.0482453(log T) 5 ] 

2.67 x 10- 1 " cxp (-^ 

1.0 x 10~ 9 

1.0 x 10- 9 cxp (-^?) 
1.0 x 10~ 9 

1.0 x 10- 9 cxp (-^?) 
1.0 x 10~ 9 



T < 2000 K 



T > 2000 K 
T ^ 200 K 
T > 200 K 

T < 2200 K 



T > 2200 K 



(-¥) 



1.0 X 10~ 9 exp 
3.5 x 10~ 10 
3.5 x lO-HTj" 019 
9.64 x 10- 10 T 3 - 024 



= 1.4 x 10" 



'exp(-^) 

-0.19. 



1.75 x lO" 1 ^" 019 exp (-if 3 -) 
2.6 x 10- 10 
8.5 x 10- 10 
8.5 x 10~ 10 

2.0 x 10- 10 exp(-^P) 



Ref. 

1 

2 
3 
2 
4 
2 
4 
2 
5 



9 
10 
4 
4 
4 
4 
4 
11 
11 
4 
12 
4 
4 
4 
13 
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Table A10 - continued 



No. 



Reaction 



Rate coefficient (cm 3 s 1 ) 



Notes Ref. 



1X36 HD^ 



1X37 HD^ 



1X38 
1X39 
1X40 
1X41 
1X42 



HD 2 
HD^ 



+ H 2 - 

+ HD 

+ HD - 
+ D 2 - 
H 2 

H 2 -+ 
HD -» 



H 2 D+ +HD 
♦ H 2 D+ +D 2 



1.0 x 10- 10 exp (-i^2) 
x [l.O + exp 



^D+4 
H 2 D+ 
HD+4 
HD+ - 



-H 2 
HD 

f D 2 
HD 

-D 2 



&IX37 

fclX38 
fclX39 
^1X40 
felX41 
&IX42 



1.0 x 10- 10 exp(-i2M) 
x [l.O + exp (-S§£)] 
: 2.0 x lO" 10 
: 8.7 x lO" 10 
1.5 x 10~ 9 exp 
1.5 x 10" 9 exp 
3.75 x 10- 10 exp 
x [l.0 + 2.0exp( 



:U2.2 

T 
233.8 

T 



155 \ 
T 7 
50.4 \ 



1X43 


HeH+ 


f D- 


-> HeD+ 


+ H 


&IX43 = 


1.0 x 10 


-9 


1X44 


HeD+ 


f H - 


-* HeH+ 


+ D 


^1X44 = 


8.0 x 10 


- 10 exp( 


1X45 


LiH+ 4 


- D - 


+ LiD+ 4 


H 


^1X45 = 


1.0 x 10 


-9 


1X46 


LiD+ 4 


-H- 


+ LiH+ 4 


D 


^1X46 = 


1.0 x 10 


~ 9 exp (— 



T J 



T ) 



- exp I 



)] 



13 

13 

4 
14 
13 
13 
13 

3 
3 
2 
2 



Notes: T is the gas temperature in K, and T3 = T/300 K 



References: 1 



Linder. Janev fe Boterol lll995[); 2 - estimate: 3 -iDalgarno fc McDowelll lligsell 



scaled as in lStancil et al. I (|l998l) : 4 - IWalmslev? Flower, fc Pineau des Foretsl <2004l) ; 5 Gcrlich 



1982) 



2003) 



1989) 



Wang fc Stancill j2002 l): 7 - Our fits to lMielke et al 



Our fits to cro ss-se ctions from 
- IShavittl fl959h: 9 - [ Millar. Bennett, fc Herbstl jl98fih : 10 - IPineau des Forets e~ 



11 - iMovano fc Collinsl ll2003h: 12 - De rived from forward reaction, using equi l ibrium 
constant from Ramanl al fc Tennvsonl l l2004h: 13 — iFlow er. Pineau des Forets. fc Walmsle vl l|2004l l: 
14 — Derived from inverse reaction in fwalmslev, FloweT!^r?ineau^^s*^br^t?720*0*4T) 
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Table All. Chemical processes: transfer reactions (TR) 



No. Reaction Rate coefficient (cm 3 s 1 ) Ref. 



TR1 


Hj + ri2 — > rig + xl 


feTRl = 


01 \* in — 9"/ 1 ().()42 
= Z.Z4 X 1U vjv-v-*" exp i 


f T 


) 1 


\ 46600 


TR2 


H J + HD — > H J + D 


^TR2 = 


= 1.05 x 10 -9 




2 


TR3 


tt4- i TTT~\ TT T"~\ — 1- i XT 

HJ + HD — > H2D^ + H 


fcTR3 = 


= 1.05 x 10~ 9 




2 


TR A 


TT + j_ pi . XJ„"n+ 1 IT. 

tl.2 -+- L)2 — > T J--* 


fc TR4 - 


= 1.05 x 10- 9 




Q 
O 


1 JaO 


TT+ i T\_, . TTTT + I "P" 

xl 2 ~r *-*2 — > XILJ2 ~r n 


&TR5 = 


= 1.05 x 10~ 9 




q 



TR a 


TTP.+ 1 T1J , TJ-\- 1 T\ 

ILL) -\- IT — ► rig T" J--* 


feTR6 = 


= 0.5 x fc T Ri 




1 


TR 7 
±1X1 


XJT)+ 1 TTT i. Tr"r,Ft+ _U TT 

JTJT -p it 2 — * JTo U ~r Jt 


^TR7 = 


= 0.5 x fc TR1 




1 
1 


1 xvo 


TTF1+ T- TTT) — *. MoT)+ 4- TT 

n U -p JT1T r JT9 JT -p IT 


^TR8 - 


= 1.05 x 10~ 9 




Q 
O 




TTH+ 4- TTT) TTFl"'" 4- TT 

JTiT ~p JT.JT * JTiT 2 T JT 


^TR9 = 


= 1.05 x 10~ 9 




Q 
O 


TR 1 n 

1 TxlU 


TTTT+ 1 TT„ . TT"T\+ 1 T~\ 
TTTJ 1 + JJ2 — *■ J ^L J -- , 2 


KTR10 


= 1.05 x 10" 9 




Q 
O 


TR "1 1 


TTTT + 1 . r»+ 1 TT 

JTJT -p JT 2 — ' JTi^ -p JT 


^TRll 


= 1.05 x 10~ 9 




Q 
O 


T P "1 O 


r\+ 1 TT , TTT TT + 1 TT 

T JT2 — TT2TJ T JT 


h 

«TR12 


= 1.05 x 10~ 9 




q 



T P "1 ^ 

1 Trio 


TT+ j_ TT ., . TT"n+ 1 TT 

J-*2 T JT2 — ► 111^2 T TT 


fc TR13 


= 1.05 x 10~ 9 




q 



TR 1 A 


n+ _l_ TTFl 1. TTTt+ _i_ n 

JT2 *t JT JT — * JTJT2 ~r *-) 


^TR14 


= 1.05 x 10~ 9 




Q 
O 


TR 1 ^ 

1 JTTO 


pl+ 1 TTTT . T)+ 1 JT 

it 2 ~p nu — >■ jTo ~~r n 


^TRIS 


= 1.05 x 10" 9 







1 JTTO 


n+ _L Tl^ i P> + __L TT 

JTo T TT2 — * TJ^ -p JT 


fcTR16 


= 2.1 x 10~ 9 




Q 
O 


1 Til ( 


TJ+ 1 TT , TJ+ 1 TT 

rig + rl — ► Jr. 2 + ir.2 


KTR17 


- 7.7 x 10- 9 cxp ( 17 ^ () ) 


A 

4 


TR18 


H3" + D — > + HD 


fcTR18 


= 0.5 X fc T R17 




5 


TR19 


TT-T 1 T~\ TTFl-l- 1 TT 

HJ + D —> HD + + H2 


fcTR19 


= 0.5 X fe T R17 




5 


TR20 


TT T~»+ 1 TT , TT+ 1 TT] v 

1I2U + rl — > H2 + tiL) 


KTR20 


= 0.5 X fe T R17 




5 


TR21 


TT T\4- 1 TJ , TJT'l ' 1 TT 

H2JJ + H — ^ rlL) n + H2 


fc TR21 


= 0.5 x fc TR i 7 




5 


TR OO 


TT -pi+ 1 T-\ T4"+ 1 

li2D T + JJ — > rl 2 + L>2 


fc TR22 


= 0.333 x fc TR i7 




r 
O 


TR '3 


X12U T" tj — * riij -f nu 


&TR23 


= 0.333 x fe T Ri7 




O 


TR24 


H 2 D+ + D — > + H 2 


,v 1 R24 


= 0.333 x fe T Ri7 




5 


TR25 


HD+ + H -r H+ + D 2 


^TR25 


= 0.333 x fe T Ri7 




5 


TR26 


HD+ + H -> HD+ + HD 


feTR26 


= 0.333 x fe T Ri7 




5 


TR27 


HD+ + H -r D+ + H 2 


^TR27 


= 0.333 x fe TR i7 




5 


TR28 


HD+ + D -> HD+ + D 2 


^TR28 


= 0.5 X fc T R17 




5 


TR29 


HD+ + D -> D+ + HD 


^TR29 


= 0.5 X fe T R17 




5 


TR30 


T"\ ~\~ 1 TT tttn-I- 1 T™\ 

D^ + H — > HD + + D2 


fcTR30 


= 0.5 X fc T R17 




5 


TR31 


D+ 4- H — > D + 4- HD 

A_J 1^ ]^ ± T " L_s 2 T^ T T i—S 


^TR31 


= 0.5 X fe T R17 




5 


TR32 


Dj + D -> D+ + D 2 

3 2 z 


^TR32 


= ^TR17 




5 


TR33 


He + H+ -> HeH+ + H 


^TR33 


= 3.0x 10- 10 exp(-^I) 


6 


TR34 


He + HD+ -» HcH+ + D 


fcTR34 


= ^TR33 




7 


TR35 


He + HD+ -» HcD+ + H 


^TR35 


= ^TR33 




7 


TR36 


He + D+ -» HcD+ + D 


^TR36 


= ^TR33 




8 


TR37 


HeH+ + H -> H+ + He 


fcxR37 


= 1.04 X 10- 9 T 3 °' 13 exp 1 




) 1 


( 33100, 


TTi 
1 ixoo 


jnejTL ^ 1 ' ' [117 ^ jne 


^TR38 


= 8.5 x 10- 10 T 3 ' 13 exp| 




) 9 


> 33100, 


TR QO 


TT„TT+ 1 TT Tt4 i TT 

xlCxl -|- xl2 — * rig ~r~ tie 


feTR39 


= 1.53 x 10" 9 T 3 °- 24 exp 1 




\ 1 


I 14800 


TR40 


HeH+ + HD — > H 2 D+ + He 


fcTR40 


= 1.20 x 10" 9 T 3 ' 24 exp 1 




) 2 


V 14800 


TR41 


HeH+ + D 2 -» HD+ + He 


^TR41 


= 1.1 X 10- 9 T0- 24 cxp (• 


T \ 
14800 ) 


10 


TR42 


TT T~\ — |— 1 TT TTT\-r- i TT 

HeD T + H — > HD T + He 


^TR42 


= 9.1 X 10" 10 T 3 013 exp( 


T ' 


) 9 


> 33100 1 


TR43 


HeD+ + D -» D+ + He 


^TR43 


= 8.5 X 10" 10 T 3 ' 13 exp| 




) 11 


» 33100 , 


TR44 


HeD+ + H 2 -> H 2 D+ + He 


^TR44 


= 1.24 X 10- 9 T 3 °' 24 exp 1 


( T 


) 2 


\ 14800 


TR45 


HcD+ + HD HD+ + He 


fcTR45 


= 1.2 X 10- 9 T 3 - 24 cxp ( 


T \ 
14800 


10 


TR46 


HeD+ + D 2 -> D+ + He 


^TR46 


= 1.1 X 10- 9 T 3 °- 24 exp ( 


T \ 
14800 ) 


10 


TR47 


LiH+ + H -> Li+ + H 2 


^TR47 


= 3.0 x lO" 10 




12 


TR48 


LiH+ + D -» Li+ + HD 


&TR48 


= 3.0 x 10" 1() 




13 


TR49 


LiD+ + H -> Li+ + HD 


^TR49 


= 3.0 x lO" 10 




14 


TR50 


LiD+ + D -> Li+ + D 2 


fcTR50 


= 3.0 x lO" 10 




14 


TR51 


LiH+ + H -> Li + H+ 


fcTR51 


= 9.0 x 10- 10 cxp ( 66400 ) 


12 


TR52 


LiH+ + D -> Li + HD+ 


^TR52 


= feTR51 




13 


TR53 


LiD+ + H Li + HD+ 


&TR53 


= ^TRSl 




14 


TR54 


LiD+ + D -> Li + D+ 


^TR54 


= &TR51 




14 
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Table All - continued 



.No. 


Reaction 




Rate coefficient (cm s ) 


net. 


TR55 


LiH- 


t-H+ 


-> Li+ + H 2 


^TRSS = 2.0 X 10" 15 


15 


TR56 


LiH- 


hD+ 


-+ Li+ + HD 


fc TR56 = 2.0 x lO" 15 


16 


TR57 


LiD- 


t- H+ 


-* Li+ + HD 


fc TR57 = 2.0 x 10~ 15 


16 


TR58 


LiD- 


t- D+ 


-* Li+ + D 2 


fc TR58 = 2.0 x 10~ 15 


16 


TR59 


LiH- 


hH+ 


-> Li + H+ 


fc T R59 = 1.0 X 10" 9 


12 


TR60 


LiH- 


hD+ 


-> Li + HD+ 


fc T R60 = 1.0 X 10" 9 


16 


TR61 


LiD- 


t- H+ 


-+ Li + HD+ 


fc T R61 = 1-0 X 10" 9 


16 


TR62 


LiD- 


b- D+ 


-> Li + D+ 


fc T R62 = LO X 10" 9 


16 


TR63 


LiH- 


hH - 


Li + H 2 


fc TR63 = 1.55 x io-"T°- 4247 


17 


TR64 


LiH- 


|-D- 


Li + HD 


fc T R 6 4 = 1-2 x io-Ht - 4247 


11 


TR65 


LiD- 


m - 


Li + HD 


fc T R65 = 1-54 x io-"T - 4247 


11 


TR66 


LiD- 




Li + D 2 


fc TR 66 = 1-2 x i -"T°- 4247 


11 



T/300 K. 
Jl995Tl: 2 



IStancil et alJ dl998T): 3 
Siilhu. Miller, fc Tennyson! dl992t) : 5 



Notes: T is the gas te mperature in K, and T3 
References: 1 iLinde r. Jane v. fc Boten 
IWalmslev. Flower, fe Pineau des Foretj ||2004|); 4 

estimate, based onlSidhu. Miller, fc TennvsonTlll992f) : 6 -IBlackl dl978f>: 7 -IStancil et al. 
based on IBlackl dl97ST): 8 — estima te, based on iBlaclq dl978l ): 9 — ILind er. Jan ev"*fc"^oterc 
scaled as in [Stancil et al . (1998) ; 10 — Estimate, based on IStancil et al .1 dl99Sft ; li- 
as corresponding H reaction, but scaled by D reduced m ass; 12 - IStancil et alJ dl996h : 
IStancil et "all lll998l1. b ased on corr esponding H reacti on in IStancil et al. I dl996l) ; 14 



1998), 



1995), 



based on IStancil et al. |_dl996l) ; 15 
17 - iDefazio et alJ <2005h 



Same 
13 - 
estimate, 



iBodo et al.l ll200ll) : 16 — same as corresponding H reaction; 



Is H% cooling important? 



Table A12. Chemical processes: background radiation induced photodctachmcnt, photodissoci- 
ation and photoionization (BP) 



No. Reaction Rate {J 2 \ s 1 ) R- e f- 



BP1 


H- + 


7 -» 


H + e" 


^BPl = 


= 1.36 X 10 


-11 


1 


BP2 


D" + 


7 -» 


D + e~ 


^?BP 2 = 


= 1.36 X 10 


-11 


2 


BP3 


H+ + 


7 -> 


H + H+ 


^BP3 = 


= 4.11 x 10 


-12 


3 


BP4 


HD+ 


+ 7" 


-^H + D+ 


i?BP4 = 


= 2.05 x 10 


-12 


2 


BP5 


HD+ 


+ 7" 


-> D + H+ 


^BP5 = 


= 2.05 x 10 


— 12 


2 


BP6 


D+ + 


7 -* 


D + D+ 


^BP6 = 


= 4.11 x 10 


— 12 


2 


BP7 


H 2 + 


7 -» 


H + H 


RBP7 - 


= 1.3 x 10" 


12 f 

/sh,H 2 


5 


BP8 


HD + 


7 -» 


H + D 


Rbps = 


= 1.45 x 10 


~ 12 /sh,HD 


6 


BP9 


D 2 + 


7 -» 


D + D 


RBP9 = 


= 1.3 x 10" 


12 


7 


BP10 


H+ + 


7 -+ 


H+ +H 


^BPIO 


= 2.4 x 10 


- 16 


8 


BP11 


H+ + 


7 -> 


H 2 +H+ 


Kbph 


= 2.4 x 10 


- 16 


8 


BP12 


H 2 D+ +7 


-> H+ + D 


^?BP1 2 


= 1.2 x 10 


-16 


9 


BP13 


H 2 D 4 


+ 7 


-» H 2 + D+ 


^BP13 


= 1.2 x 10 


-16 


9 


BP 14 


H 2 D+ 


+ 7 


-> HD+ + H 


Rbpu 


= 1.2 x 10 


-16 


9 


BP15 


H 2 D+ 


+ 7 


-» HD + H+ 


^?BP15 


= 1.2 x 10 


-16 


9 


BP16 


HD+ 


+ 7" 


-> HD+ +D 


RBP16 


= 1.2 x 10 


-16 


9 


BP17 


HD+ 


+ 7" 


-> HD + D+ 


KBP17 


= 1.2 x 10 


-16 


9 


BP18 


HD+ 


+ 7 - 


-^D++H 


-KBP18 


= 1.2 x 10" 


-16 


9 


BP19 


HD+ 


+ 7" 


^D 2 +H+ 


^BP19 


= 1.2 x 10 


-16 


9 


BP20 


D+ + 


7 -> 


D+ +D 


fiBP 2 () 


= 2.4 x 10 


-16 


9 


BP21 


4 + 


7 -» 


D 2 +D+ 


^?BP 2 1 


= 2.4 x 10 


-16 


9 


BP22 


HeH+ 


+ 7 


-> He + H+ 


^BP 22 


= 1.0 x 10 


-17 


10 


BP23 


HeD 4 


+ 7 


-> He + D+ 


^BP 2 3 


= 1.0 x 10 


-17 


10 


BP24 


He+- 


^7- 


+ He + He+ 


RBP24, 


= 1.0 x 10 


-12 


11 


BP25 


Li + 7 — ► Li+ + e - 


J?BP 2 5 


= 1.4 x 10 


-12 


12 


BP26 


Li~ + 


7 -» 


Li + e - 


^BP 2 6 


= 1.2 x 10 


-11 


13 


BP27 


LiH+ 


+ 7 


-» Li+ + H 


^BP 2 7 


= 5.0 x 10 


-18 


13 


BP28 


LiH+ 


+ 7 


-» Li + H+ 


^?BP 2 8 


= 9.3 x 10 


-9 


13 


BP29 


LiD+ 


+ 7 


— ► Li+ + D 


RBP29 


= 5.0 x 10 


-18 


2 


BP30 


LiD+ 


+ 7 


-» Li + D+ 


^?BP30 


= 9.3 x 10 


-9 


2 


BP31 


LiH 4 


7 -4 


Li + H 


^BP31 


= 4.4 x 10 


-14 


14 


BP32 


LiD ■+ 


7 -• 


Li + D 


^BP32 


= 4.4 x 10 


-14 


2 



Notes: 7 represents a photon from the external background radiation field. The listed reaction 
rates were computed assuming that this background has the spectrum of a 10 5 K diluted 
black-body, cut-off above hv = 13.6 eV, as described in Ej3] With this spectrum, reactions with 
threshold energies greater than 13.6 eV do not occur and are not listed in the table. / s h,H 2 
and /sh.Hn are the self-shie lding factors for H 2 and HD photodissociation, respectively (see 
e.g., [Glover fe Jappsenl I2007T) . Note that in this paper, we consider only the limiting cases 

/sh,H 2 = /sh.HD = and / s h.Ho = /sh.HD = !• 

References: 1 - Wishart (1979); 2 — assumed same as for corresponding H reaction; 3 - iDunnl 
(|l968h : 4 — total rate a ssumed same as for c orresp ondi ng H reaction, indiv i dual o utcomes assumed 
e nually probabl e ; 5 - iDraine fe Bertoldil dl996l); 6 - lAbgrall fc Rouefj 1 120061); 7 - estimate; 8 
van Dishoeckl <198Sl ) : 9 — estimat e, based on Ivan DishoeckT (|l988l l : 10 -iRoberge fc Dalgarnol 
Jl982h: 11 - IStancill il994h ; 12 - IVerner fc Ferlandl jl996ft ; 13 - iGalli fc Pallal lll998f ); 1~ 
iKirbv fc Dalgarnd <197sF 
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Table A13. Chemical processes: cosmic ray ionization (CR) 



No. Process Rate (Ci/Ch) Reference 



CR1 


H + 


CR. - 


► H+ + e- 


1.0 




CR2 


H 2 - 


f CR. - 


-+ H+ +e" 


2.09 


1 


CR3 


H 2 - 


f- CR. - 


-> H + H+ + e~ 


0.09 


1 


CR4 


H 2 - 


f CR. - 


-» H + H 


3.26 


1 


CR5 


He- 


h CR. - 


-> He+ + e" 


1.09 


1 


CR6 


D + 


CR. - 


» D+ +c- 


1.0 


2 


CR7 


HD 


+ C.R. 


-> HD+ + e- 


2.09 


2 


CR8 


HD 


+ C.R. 


-> H + D+ +e" 


0.04 


2 


CR9 


HD 


+ C.R. 


-> H+ + D + e- 


0.04 


2 


CR10 


HD 


+ C.R. 


^H + D 


3.26 


2 


CR11 


D 2 - 


f CR. - 


-+ T>2 + e_ 


2.09 


2 


CR12 


D 2 - 


f CR. - 


-> D + D+ +e" 


0.09 


2 


CR13 


D 2 - 


t- CR. - 


-^D + D 


3.26 


2 



Notes: CR. represents a cosmic ray. £h, the cosmic ray ionization rate of atomic hydrogen, is 
an adjustable par ameter in our models. 

References: 1 - Walmsle v. Flower. fc Pineau d cs Forfits (2004); 2 - assumed same as corre- 
sponding H process 



Table A14. Chemical processes: cosmic ray induced photodetachment, photodissociation and 
photoionization (CP) 



No. 


Reaction 


°-X,eff,H 2 ( Mb ) 


^X.cff.H (Mb) 


Ref. 


CP1 


H" + 7„ -» H + e- 


5.0 


5.8 


1 


CP2 


D~ + 7 cr -» D + c" 


5.0 


5.8 


2 


CP3 


H+ + 7cr H + H+ 


5.0 


6.6 


3 


CP4 


HD+ + 7„ -> H + D+ 


2.5 


3.3 


2 


CP 5 


HD+ + 7cr -» D + H+ 


2.5 


3.3 


2 


CP6 


D+ + 7cl — > D + D+ 


5.0 


6.6 


2 


CP 7 


Li + 7cr — * Li + + e~ 


1.0 


1.3 


4 


CP8 


Li" + 7 cr — * Li + e~ 


1.0 


1.0 


5 


CP9 


He+ + 7cr — He + He+ 


5.0 


5.0 


6 


CP10 


LiH+ + 7cr -> Li + H+ 


100 


100 


7 


CPU 


LiD+ + 7cr -> Li + D+ 


100 


100 


2 



Notes: 7 cr represents a secondary photon, produced by cosmic-ray induced excitation of H or H 2 , 
as described in i|3,3l The references listed are the sources from which we have taken our photodis- 
sociation or photoionization cross-sections. The emission probabil ities Phq {v) used to calculate 
ox.eff .Ho are rough estimates based on the emission spectra given in lSternberg. Dalgarno. fc LeprJ 

lll987ti and are li kely accurate on ly to within a factor of a few. 

Refer ence s: 1 -IWisharj lll979h: 2 — assumed same as for corresponding H reaction; 3 - iDunnl 
lll968l); 4 -IVerner fc Ferlandl J1996T) ; 5 — order of magnitu de estimate; 6 - esti mate, based on 
Stanci 3 l ll994h ; 7 - rough estimate, based on thermal rate in iGalli fc Palial dl998l) 



